<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">Hi 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">The attached revised input converged after a few steps using SCCS with a dielectric constant of 80. I changed the method and adapted the rho_min and rho_max values which define the surface of the cavity and thus they are system dependent to some extent. If rho_max is chosen too large for instance, the surface might be too close to the molecule. It is often helpful to plot the cavity, e.g. with VMD, using a cube file dump of the <a href="https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/SCCS/DIELECTRIC_FUNCTION.html"> dielectric function</a>.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#1F497D">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Matthias</span><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"><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> Freitag, 15. Mai 2020 17:47<br> <b>To:</b> cp2k <...@googlegroups.com><br> <b>Subject:</b> Re: [CP2K:13314] 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">Hi Matthias,<o:p></o:p></p> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">Thank you again for all your help. I went a bit beyond your suggestion and changed EPS_SCF to 1E-12 and changed cell size to 10 x 10 x 10 A3. To answer your questions:<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">1.) I was able to run GEO_OPT without SCCS and the simulation completed. <o:p></o:p></p> </div> <div> <p class="MsoNormal">2.) I was able to converge and complete a SCCS GEO_OPT simulation with a dielectric constant of 1. <o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">I then took the dieletric constant back up to 80 (with no other changes), and the convergence fails - in fact the convergence Hartree's are trending greater than 1. <o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">I'm going to do a small ramp up in dielectric constant: 2, 5, 10, etc. To see what happens. But I should be able to use 80+ as the dielectric constant without issues.<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">What do you think? Am I missing something?<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">Again input, output, and XYZ attached. <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 SCCS<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"> SAVE_MEM TRUE<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-12<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 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> CUTOFF [Ry] 600<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<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] 0.0<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"> GAMMA 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_SCCS 1.0E-12<o:p></o:p></p> </div> <div> <p class="MsoNormal">! EPS_SCF 1.0E-10<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.008<o:p></o:p></p> </div> <div> <p class="MsoNormal">! RHO_MIN 0.00015<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-12<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &DIAGONALIZATION ON<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">! &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_ALL<o:p></o:p></p> </div> <div> <p class="MsoNormal">! ENERGY_GAP 0.002<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &END OT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &OUTER_SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_SCF 1e-12<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MAX_SCF 30<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END OUTER_SCF<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 off<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 off<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 off<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"> TYPE MINIMIZATION<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"> &CONSTRAINT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &FIXED_ATOMS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> COMPONENTS_TO_FIX XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal"> LIST 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END FIXED_ATOMS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END CONSTRAINT<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> <p class="MsoNormal"><o:p> </o:p></p> </div> <p class="MsoNormal"><br> On Wednesday, May 13, 2020 at 8:08:47 PM UTC-4, Matthias Krack wrote:<o:p></o:p></p> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi Dev<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">The settings in the regression tests are often physically not meaningful, since these test runs have to be completed after a few seconds.<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">SCCS%EPS_SCF<span style="color:#1F497D"> </span>delays the start of the inner SCCS convergence loop only for values greater zero. Otherwise that value is not considered and the inner SCCS cycle is activated immediately. <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">SCCS%EPS_SCCS, the convergence threshold of the inner SCCS loop for the polarization charge should be smaller (tighter) than 1.0E-6, at least 1.0E-8. As soon as the inner SCCS loop is activated, the timings for each SCF iteration step will increase significantly. Also the implicit solvent comes not for free. You can reduce the box size to save computer resources, because your system is small. Does the setup for your system CAl4 work/converge properly for a geo_opt run in gas phase, i.e. without SCCS? If that is the case then I would start with a smaller dielectric constant than 80 for the SCCS and see how the system evolves. The print level medium provides a detailed SCCS output.<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">HTH<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Matthias<o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><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"> <a href="javascript:" target="_blank">c...@googlegroups.com</a> <<a href="javascript:" target="_blank">c...@googlegroups.com</a>> <b>On Behalf Of </b>Dev Rana<br> <b>Sent:</b> Donnerstag, 14. Mai 2020 00:29<br> <b>To:</b> cp2k <<a href="javascript:" target="_blank">c...@googlegroups.com</a>><br> <b>Subject:</b> Re: [CP2K:13298] SCCS Solvation Model Implementation - Molecule Radii Explodes</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I take back what I previously said. The 0th geo opt step finishes quickly with SCF steps taking 20-25s each. The 2st geo opt step does not finish and each SCF step takes 5+ minutes. Here's my input file and output for that. I'm also running this on GPU. <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">What do you advise? I've attached the input output files and paste input below. The convergence values keep increases, rather than decreases. And Total Energy keep lowering. <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&GLOBAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> PROJECT SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> RUN_TYPE geo_opt<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> PRINT_LEVEL low<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> SAVE_MEM TRUE<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&END GLOBAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&FORCE_EVAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> METHOD Quickstep<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> # STRESS_TENSOR analytical<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &DFT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> BASIS_SET_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/BASIS_MOLOPT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> POTENTIAL_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/POTENTIAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &QS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> METHOD GPW<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EPS_DEFAULT 1.0E-12<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EXTRAPOLATION ASPC<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END QS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &MGRID<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> NGRIDS 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> CUTOFF [Ry] 600<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END MGRID<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ALPHA [N*m^-1] 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> BETA [GPa] 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DELTA_RHO 2.0E-5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DERIVATIVE_METHOD cd5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DIELECTRIC_CONSTANT 80<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GAMMA [mN/m] 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EPS_SCCS 1.0E-6<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EPS_SCF 1.0E-7<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> MAX_ITER 100<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! METHOD Andreussi<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> METHOD Fattebert-Gygi<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> MIXING 0.6<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! &ANDREUSSI<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! RHO_MAX 0.001<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! RHO_MIN 0.0001<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! &END ANDREUSSI<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &FATTEBERT-GYGI<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> BETA 1.3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> RHO_ZERO 0.00078<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END FATTEBERT-GYGI<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> SCF_GUESS ATOMIC<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> MAX_SCF 300<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EPS_SCF 1.0E-7<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &DIAGONALIZATION ON<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ALGORITHM STANDARD<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END DIAGONALIZATION<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! &OT on<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! MINIMIZER DIIS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! PRECONDITIONER FULL_ALL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! ENERGY_GAP 0.002<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! &END OT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &OUTER_SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> EPS_SCF 1e-7<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> MAX_SCF 30<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END OUTER_SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &XC<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &XC_FUNCTIONAL PBE<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END XC_FUNCTIONAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END XC<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &PDOS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> APPEND TRUE<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END PDOS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &SCCS ON<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &DENSITY_GRADIENT off<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END DENSITY_GRADIENT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &DIELECTRIC_FUNCTION off<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END DIELECTRIC_FUNCTION<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &POLARISATION_POTENTIAL off<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END POLARISATION_POTENTIAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END DFT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &SUBSYS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ABC [angstrom] 30 30 30<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! PERIODIC XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! MULTIPLE_UNIT_CELL 1 1 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &TOPOLOGY<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> COORD_FILE_NAME Al4C.xyz<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> COORD_FILE_FORMAT XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">! MULTIPLE_UNIT_CELL 1 1 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &KIND C<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ELEMENT C<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> BASIS_SET DZVP-MOLOPT-GTH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> POTENTIAL GTH-PBE-q4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END KIND<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &KIND Al<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ELEMENT Al<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> BASIS_SET DZVP-MOLOPT-SR-GTH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> POTENTIAL GTH-PBE-q3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END KIND<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END SUBSYS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&END FORCE_EVAL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&MOTION<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &GEO_OPT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> TYPE MINIMIZATION<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> OPTIMIZER BFGS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> MAX_ITER 100<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END GEO_OPT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &CONSTRAINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &FIXED_ATOMS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> COMPONENTS_TO_FIX XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> LIST 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END FIXED_ATOMS<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END CONSTRAINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &STRUCTURE_DATA<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 1 2<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 1 3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 1 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 1 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 2 3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 2 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 2 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 3 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 3 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> DISTANCE 4 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 2 3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 2 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 2 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 3 2<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 3 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 3 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 4 2<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 4 3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 4 5<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 5 2<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 5 3<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> ANGLE 1 5 4<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END STRUCTURE_DATA<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &TRAJECTORY<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END TRAJECTORY<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &VELOCITIES<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END VELOCITIES<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &FORCES<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END FORCES<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &RESTART<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END RESTART<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> &END PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">&END MOTION<o:p></o:p></p> </div> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p> </div> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">-- <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="javascript:" target="_blank">c...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b0c10452-d528-48da-9f83-07d9b01d5dfb%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/b0c10452-d528-48da-9f83-07d9b01d5dfb%40googlegroups.com</a>.<o:p></o:p></p> </div> </div> </blockquote> </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/cc63d43c-380a-4575-ac41-5b981251d294%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/cc63d43c-380a-4575-ac41-5b981251d294%40googlegroups.com</a>.<o:p></o:p></p> </div> </body> </html>