<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>See, but that's the funny part. I absolutely DO NOT want any
topologies created and I did not think there was any need for it
in this type of simulations. Other than using a PDB file with
coordinates as input as part TOPOLOGY section, I am not asking
(not explicitly, anyway) for any connectivity tables. should I
just convert to xyz and dump the whole TOPOLOGY section? All
suggestions are very welcome.<br>
</p>
<p>Thanks,</p>
<p>Alex<br>
</p>
<p>p.s. Complete file below:<br>
</p>
&GLOBAL<br>
 PROJECT opt<br>
 PRINT_LEVEL MEDIUM <br>
 RUN_TYPE GEO_OPT<br>
&END GLOBAL<br>
 &MOTION<br>
 &GEO_OPT<br>
   MAX_ITER 400<br>
   OPTIMIZER BFGS <br>
 &END<br>
&END MOTION <br>
&FORCE_EVAL<br>
 METHOD QS<br>
 STRESS_TENSOR ANALYTICAL<br>
 &DFT<br>
#Â Â Â CHARGE 1Â Â Â <br>
   BASIS_SET_FILE_NAME BASIS_SET2<br>
   POTENTIAL_FILE_NAME GTH_POTENTIALS2<br>
   &MGRID<br>
     NGRIDS 5<br>
     CUTOFF 500<br>
   &END MGRID<br>
   &QS<br>
     METHOD GPW<br>
   &END QS<br>
   &SCF<br>
     SCF_GUESS RESTART <br>
     EPS_SCF 5.0E-6<br>
     MAX_SCF 150<br>
     ADDED_MOS 100 <br>
     &MIXING<br>
         METHOD BROYDEN_MIXING<br>
         ALPHA   0.1<br>
         NBROYDEN  8<br>
     &END MIXING<br>
   &END SCF<br>
    &XC<br>
       &XC_FUNCTIONAL PBE<br>
       &END XC_FUNCTIONAL<br>
       &VDW_POTENTIAL<br>
         POTENTIAL_TYPE PAIR_POTENTIAL<br>
         &PAIR_POTENTIAL<br>
           TYPE DFTD3<br>
           PARAMETER_FILE_NAME dftd3.dat<br>
          REFERENCE_FUNCTIONAL PBE<br>
          CALCULATE_C9_TERM T<br>
         &END PAIR_POTENTIAL<br>
       &END VDW_POTENTIAL<br>
     &END XC<br>
 &END DFT<br>
 &SUBSYS<br>
   &CELL<br>
    ABC   22.087 21.859 50.000<br>
   ANGLES   90.000 90.000 90.000<br>
       PERIODIC XYZ <br>
       MULTIPLE_UNIT_CELL 1 1 1 <br>
   &END<br>
   &TOPOLOGY<br>
      &CENTER_COORDINATES TRUE<br>
      &END<br>
      COORD_FILE_NAME pore3_potassium4.pdb<br>
      COORD_FILE_FORMAT PDB<br>
      MULTIPLE_UNIT_CELL 1 1 1<br>
   &END<br>
   &KIND S<br>
     BASIS_SET DZVP-MOLOPT-SR-GTH<br>
     POTENTIAL GTH-PBE-q6<br>
   &END KIND<br>
   &KIND MO<br>
     BASIS_SET DZVP-MOLOPT-SR-GTH<br>
     POTENTIAL GTH-PBE-q14<br>
   &END KIND<br>
   &KIND O<br>
     BASIS_SET DZVP-MOLOPT-SR-GTH<br>
     POTENTIAL GTH-PBE-q6<br>
   &END KIND<br>
   &KIND K <br>
     BASIS_SET DZVP-MOLOPT-SR-GTH<br>
     POTENTIAL GTH-PBE-q9<br>
   &END KIND<br>
&END SUBSYS<br>
&END<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">On 11/23/2017 2:33 AM, Gabriele wrote:<br>
</div>
<blockquote type="cite"
cite="mid:34cc39e3-80ff-4d...@googlegroups.com">
<div dir="ltr">Hi Alex,
<div><br>
</div>
<div>It looks like the error comes from the SUBSYS section
in <a class="moz-txt-link-freetext" href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/">https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/</a>
and not from the DFT section. If you are using just DFT there
is no need to generate the connectivity, which is where you
get the warning. We can have a closer look if you upload the
whole input file. </div>
<div><br>
</div>
<div>Cheers,</div>
<div>Gabriele</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div><br>
On Thursday, November 23, 2017 at 3:02:52 AM UTC+1, Alex
wrote:
<blockquote class="gmail_quote" style="margin:
0;margin-left: 0.8ex;border-left: 1px #ccc
solid;padding-left: 1ex;">
<div dir="ltr">Hi all,
<div><br>
</div>
<div>I have a system that contains a molybdenum
disulfide (MoS2) sample and a charged (+1) potassium
ion. The MoS2 structure itself gets energy-minimized
without any issues, but in combination with this ion I
get a nonstop stream of GENERATE warnings like this:Â </div>
<div><br>
</div>
<div>
<div><br>
</div>
<div>Â GENERATE|Â Preliminary Number of Bonds
generated:Â Â Â Â Â Â Â Â Â Â Â Â Â Â 360</div>
<div><br>
</div>
<div>Â GENERATE| WARNING in connectivity generation!</div>
<div>Â GENERATE| Two molecules/residues named (MOL1)
have different number of atoms.</div>
<div>Â GENERATE| Molecule starting at position (13) has
Nr. <0> of atoms.</div>
<div>Â GENERATE| while the other same molecules have
Nr. <162> of atoms!</div>
<div>Â GENERATE| Increasing bondparm_factor by 1.05..
An error was found in the generated</div>
<div>Â GENERATE| connectivity. Retry...</div>
<div>Â GENERATE| Present value of BONDPARM_FACTOR (Â
 Infinity ).</div>
</div>
<div><br>
</div>
<div>When I replace the potassium ion with e.g. a
proton, everything is back to normal.</div>
<div><br>
</div>
<div>Here's the DFT section of my input:</div>
<div><br>
</div>
<div>
<div>Â &DFT</div>
<div>Â Â CHARGE 1</div>
<div>Â Â BASIS_SET_FILE_NAME BASIS_SET2</div>
<div>Â Â POTENTIAL_FILE_NAME GTH_POTENTIALS2</div>
<div>Â Â &MGRID</div>
<div>Â Â Â NGRIDS 5</div>
<div>Â Â Â CUTOFF 500</div>
<div>Â Â &END MGRID</div>
<div>Â Â &QS</div>
<div>Â Â Â METHOD GPW</div>
<div>Â Â &END QS</div>
<div>Â Â &SCF</div>
<div>Â Â Â SCF_GUESS RESTART</div>
<div>Â Â Â EPS_SCF 5.0E-6</div>
<div>Â Â Â MAX_SCF 150</div>
<div>Â Â Â ADDED_MOSÂ 100</div>
<div>Â Â Â &MIXING</div>
<div>Â Â Â Â Â METHOD BROYDEN_MIXING</div>
<div>Â Â Â Â Â ALPHAÂ Â 0.1</div>
<div>Â Â Â Â Â NBROYDENÂ Â 8</div>
<div>Â Â Â &END MIXING</div>
<div>Â Â &END SCF</div>
<div>Â Â Â &XC</div>
<div>Â Â Â Â &XC_FUNCTIONALÂ PBE</div>
<div>Â Â Â Â &END XC_FUNCTIONAL</div>
<div>Â Â Â Â &VDW_POTENTIAL</div>
<div>Â Â Â Â Â POTENTIAL_TYPEÂ PAIR_POTENTIAL</div>
<div>Â Â Â Â Â &PAIR_POTENTIAL</div>
<div>Â Â Â Â Â Â TYPEÂ DFTD3</div>
<div>Â Â Â Â Â Â PARAMETER_FILE_NAME dftd3.dat</div>
<div>Â Â Â Â Â Â REFERENCE_FUNCTIONAL PBE</div>
<div>Â Â Â Â Â Â CALCULATE_C9_TERMÂ T</div>
<div>Â Â Â Â Â &END PAIR_POTENTIAL</div>
<div>Â Â Â Â &END VDW_POTENTIAL</div>
<div>Â Â Â &END XC</div>
<div>Â &END DFT</div>
</div>
<div><br>
</div>
<div>Any suggestions? If you need further info, please
let me know. </div>
<div><br>
</div>
<div>Thank you,</div>
<div><br>
</div>
<div>Alex</div>
</div>
</blockquote>
</div>
</div>
</div>
-- <br>
You received this message because you are subscribed to a topic in
the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a
href="https://groups.google.com/d/topic/cp2k/Nsy-57sQOJ4/unsubscribe"
moz-do-not-send="true">https://groups.google.com/d/topic/cp2k/Nsy-57sQOJ4/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email
to <a href="mailto:cp2k+uns...@googlegroups.com"
moz-do-not-send="true">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a
href="mailto:cp...@googlegroups.com" moz-do-not-send="true">cp...@googlegroups.com</a>.<br>
Visit this group at <a
href="https://groups.google.com/group/cp2k"
moz-do-not-send="true">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a
href="https://groups.google.com/d/optout" moz-do-not-send="true">https://groups.google.com/d/optout</a>.<br>
</blockquote>
<br>
</body>
</html>