<p class="separator" style="text-align: center; clear: both;"><br></p><p class="separator" style="text-align: center; clear: both;"><br></p><p class="separator" style="text-align: center; clear: both;"><br></p><p class="separator" style="text-align: center; clear: both;"><a imageanchor="1" href="https://lh3.googleusercontent.com/-dyTmC1pMzY0/WjexEXeyflI/AAAAAAAADfA/-imk0ASHQYkR4Oi3byRMlJUFcoRZA1fGQCLcBGAs/s1600/.png" style="margin-left: 1em; margin-right: 1em;"><img src="https://lh3.googleusercontent.com/-dyTmC1pMzY0/WjexEXeyflI/AAAAAAAADfA/-imk0ASHQYkR4Oi3byRMlJUFcoRZA1fGQCLcBGAs/s320/.png" border="0" style="" width="320" height="292"></a><a imageanchor="1" href="https://lh3.googleusercontent.com/-NEGCu6ANRYs/WjexBaiZkSI/AAAAAAAADe8/pnQvzWDI8q45VJ6yYnXKbCMYXoVsgVuMQCLcBGAs/s1600/cp2k11.png" style="margin-left: 1em; margin-right: 1em;"><img src="https://lh3.googleusercontent.com/-NEGCu6ANRYs/WjexBaiZkSI/AAAAAAAADe8/pnQvzWDI8q45VJ6yYnXKbCMYXoVsgVuMQCLcBGAs/s320/cp2k11.png" border="0" width="301" height="320"></a></p><div>Hi CP2K group,</div><div><br></div>I am trying to create a PSF file for ALOOH surface. I have a particular requirement that the psf file should have connections shown as in left figure.<div><br></div><div>I have created the PSF file, but CP2K shows error: </div><div>CP2K requires molecules to be contiguous and we have detected a non *                                                                                                                               <div> * [ABORT]    contiguous one!! In particular a molecule defined from index (   *                                                                                                                               </div><div> *  \___/  1) to (   210) contains other molecules, not connected! </div><div><br></div><div><br></div><div>My PSF file look like this: </div><div><br></div><div><div>     256 !NATOM                                                                                                                                                                                                </div><div>       1 14   01   DIAS Al1  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       2 14   01   DIAS Al2  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       3 03   01   DIAS Al3  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       4 03   01   DIAS Al4  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       5 13   01   DIAS Al5  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       6 13   01   DIAS Al6  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       7 02   01   DIAS Al7  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       8 02   01   DIAS Al8  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       9 12   01   DIAS Al9  YA     0.000000        0.0000           0                                                                                                                                         </div><div>      10 12   01   DIAS Al1A YA     0.000000        0.0000           0                                                                                                                                         </div><div>      11 01   01   DIAS Al1B YA     0.000000        0.0000           0 </div></div><div><br></div><div>Questions:</div><div>1) If it is same molecule, then in PSF file should SEGMENT ID should be defined exactly in a row? for instance</div><div><br></div><div>which format is correct for CP2K. See SEGMENT ID (second column).. Should each molecule be defined in a row or it can be expressed randomly?</div><div><br></div><div><div> 256 !NATOM                                                                                                                                                                                                </div><div>       1 AA   01   DIAS Al1  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       2 BB  01   DIAS Al4  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       3 AA  01   DIAS Al3  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       4 AA 01   DIAS Al2  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       5 BB  01   DIAS Al5  YA     0.000000        0.0000           0      </div></div><div><br></div><div><br></div><div><div> 256 !NATOM                                                                                                                                                                                                </div><div>       1 AA   01   DIAS Al1  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       2 AA  01   DIAS Al2  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       3 AA  01   DIAS Al3  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       4 BB  01   DIAS Al4  YA     0.000000        0.0000           0                                                                                                                                         </div><div>       5 BB  01   DIAS Al5  YA     0.000000        0.0000           0      </div></div><div><br></div><div> </div><div>2) In my current ALO(OH), except for the non-hydroxy oxygen everything are connected. I could create PSF connecting all atoms (easier to create for me compared to the PSF separating the non-hydroxy oxygen). Is there any easier way to define non-hydroxy oxygens as separate molecules, so that I can have 1-3 VdW, EI.</div><div><br></div><div><br></div><div>Please let me know. Attached are my PSF, input and PDB files.</div><p class="separator" style="text-align: center; clear: both;"><br></p><div><br></div><div><br></div><div><br>On Saturday, February 14, 2015 at 2:53:48 PM UTC+1, Luca wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Christian,
<br>I have investigated in your problem,
<br>Maybe the solution is in the definition of molecule 
<br>in the PSF file . I have just reproduced your error using a small
<br>system.
<br>The original PSF created with psfgen for 10 waters and 2 ions (SOD and
<br>CLA) does not work in cp2k as it is. 
<br>
<br> 29 !NATOM
<br>   [..]   
<br>      26 W1   9    TIP3 H1   HT     0.417000        1.0080           0
<br>      27 W1   9    TIP3 H2   HT     0.417000        1.0080           0
<br>      28 O1   10   CLA  CLA  CLA   -1.000000       35.4500           0
<br>      29 O1   11   SOD  SOD  SOD    1.000000       22.9898           0
<br>   [..]
<br>
<br>Because, cp2k read info about the molecule from segment definition the
<br>error arises from the fact that CLA and SOD have same fragment/segment
<br>name, which is "O1".
<br>If you change this name in the PSF  file, it should be work.
<br>
<br>      28 O1   10   CLA  CLA  CLA   -1.000000       35.4500           0
<br>      29 O2   11   SOD  SOD  SOD    1.000000       22.9898           0
<br>    
<br>Good luck,
<br>
<br>Luca    
<br>
<br>
<br>
<br>
<br>> My apologies. I've reported it now.
<br>> 
<br>> 
<br>> On 12 Feb 2015 21:47, "Teodoro Laino" <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">teod...@gmail.com</a>> wrote:
<br>>         Which tone? If you believe that you have a problem with psfgen
<br>>         you should report this to the psfgen mailing list.
<br>>         
<br>>         > On 12 Feb 2015, at 15:07, Christian Jorgensen
<br>>         > <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">chr...@gmail.com</a>> wrote:
<br>>         > 
<br>>         > Sorry I do not understand why the tone on
<br>>         > this board is like this. I'm simply asking for help.
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > On Thu, Feb 12, 2015 at 7:16 AM, Teodoro Laino
<br>>         > <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">teod...@gmail.com</a>> wrote:
<br>>         >         report this to the proper mailing list.
<br>>         >         
<br>>         >         > On 11 Feb 2015, at 23:04, Christian Jorgensen
<br>>         >         > <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">chr...@gmail.com</a>> wrote:
<br>>         >         > 
<br>>         >         > ATOM  49880  SOD SOD I  37      21.151   0.526
<br>>         >         > -30.628  1.00  0.00      ION NA
<br>>         >         > ATOM  49881  SOD SOD I  38      24.667 -35.051
<br>>         >         >  19.114  1.00  0.00      ION NA
<br>>         >         > ATOM  49882  SOD SOD I  39     -12.620   2.426
<br>>         >         >  26.969  1.00  0.00      ION NA
<br>>         >         > ATOM  49883  SOD SOD I  40     -10.705 -33.380
<br>>         >         > -11.158  1.00  0.00      ION NA
<br>>         >         > ATOM  49884  SOD SOD I  41     -19.491  29.350
<br>>         >         > -34.679  1.00  0.00      ION NA
<br>>         >         > ATOM  49885  SOD SOD I  42      40.802 -30.866
<br>>         >         >  10.911  1.00  0.00      ION NA
<br>>         >         > ATOM  49886  SOD SOD I  43     -21.411 -20.389
<br>>         >         >  17.270  1.00  0.00      ION NA
<br>>         >         > ATOM  49887  SOD SOD I  44      13.100  24.287
<br>>         >         >  18.857  1.00  0.00      ION NA
<br>>         >         > ATOM  49888  SOD SOD I  45      33.628  30.716
<br>>         >         >  18.333  1.00  0.00      ION NA
<br>>         >         > ATOM  49889  SOD SOD I  46      14.256 -31.556
<br>>         >         >  -9.097  1.00  0.00      ION NA
<br>>         >         > ATOM  49890  SOD SOD I  47      -1.600  38.825
<br>>         >         >  28.079  1.00  0.00      ION NA
<br>>         >         > ATOM  49891  SOD SOD I  48      26.911  -9.319
<br>>         >         >  19.217  1.00  0.00      ION NA
<br>>         >         > ATOM  49892  CLA CLA I  49     -21.849 -29.233
<br>>         >         >  -4.389  1.00  0.00      ION CL
<br>>         >         > ATOM  49893  CLA CLA I  50     -27.778  19.207
<br>>         >         >  -8.907  1.00  0.00      ION CL
<br>>         >         > ATOM  49894  CLA CLA I  51      -1.036 -18.525
<br>>         >         > -33.557  1.00  0.00      ION CL
<br>>         >         > ATOM  49895  CLA CLA I  52     -37.955  23.751
<br>>         >         >  -9.529  1.00  0.00      ION CL
<br>>         >         > ATOM  49896  CLA CLA I  53     -15.636 -34.005
<br>>         >         >  12.873  1.00  0.00      ION CL
<br>>         >         > ATOM  49897  CLA CLA I  54      40.122  21.296
<br>>         >         > -26.349  1.00  0.00      ION CL
<br>>         >         > ATOM  49898  CLA CLA I  55      32.593  -9.101
<br>>         >         >  17.836  1.00  0.00      ION CL
<br>>         >         > ATOM  49899  CLA CLA I  56       8.168 -24.215
<br>>         >         > 8.246  1.00  0.00      ION CL
<br>>         >         > ATOM  49900  CLA CLA I  57      -4.530 -32.691
<br>>         >         >  -6.031  1.00  0.00      ION CL
<br>>         >         > ATOM  49901  CLA CLA I  58     -14.998  37.422
<br>>         >         > 2.603  1.00  0.00      ION CL
<br>>         >         > ATOM  49902  CLA CLA I  59     -22.559 -24.866
<br>>         >         > -37.160  1.00  0.00      ION CL
<br>>         >         > ATOM  49903  CLA CLA I  60      22.442 -17.653
<br>>         >         > -37.184  1.00  0.00      ION CL
<br>>         >         > ATOM  49904  CLA CLA I  61     -31.657   0.756
<br>>         >         > -22.411  1.00  0.00      ION CL
<br>>         >         > ATOM  49905  CLA CLA I  62     -25.707   4.569
<br>>         >         >  26.200  1.00  0.00      ION CL
<br>>         >         > ATOM  49906  CLA CLA I  63       9.816  21.844
<br>>         >         > -29.523  1.00  0.00      ION CL
<br>>         >         > ATOM  49907  CLA CLA I  64      20.348  15.655
<br>>         >         >  31.373  1.00  0.00      ION CL
<br>>         >         > ATOM  49908  CLA CLA I  65     -33.471  28.927
<br>>         >         >  25.816  1.00  0.00      ION CL
<br>>         >         > ATOM  49909  CLA CLA I  66     -23.076 -36.320
<br>>         >         > -14.581  1.00  0.00      ION CL
<br>>         >         > ATOM  49910  CLA CLA I  67      11.228  23.055
<br>>         >         >  26.442  1.00  0.00      ION CL
<br>>         >         > ATOM  49911  CLA CLA I  68     -31.967 -25.220
<br>>         >         >  12.122  1.00  0.00      ION CL
<br>>         >         > ATOM  49912  CLA CLA I  69       7.122 -32.163
<br>>         >         > -23.902  1.00  0.00      ION CL
<br>>         >         > ATOM  49913  CLA CLA I  70     -30.542 -17.921
<br>>         >         >  15.502  1.00  0.00      ION CL
<br>>         >         > ATOM  49914  CLA CLA I  71     -35.163 -32.044
<br>>         >         >  11.857  1.00  0.00      ION CL
<br>>         >         > ATOM  49915  CLA CLA I  72      -1.340 -18.251
<br>>         >         > -21.790  1.00  0.00      ION CL
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > In relation to the error
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         >        49892           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         >        49893           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         >        49894           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         >        49895           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         >        49896           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         >        49897           0
<br>>         >         >  Two molecules have been defined as identical
<br>>         >         > molecules but atoms mismatch charges!!
<br>>         >         > "
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > How can this be? I've generated this with the
<br>>         >         > standard PSFGEN / Ionize plugin written
<br>>         >         > for VMD. 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > On Wed, Feb 11, 2015 at 9:14 PM, Christian
<br>>         >         > Jorgensen <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">chr...@gmail.com</a>> wrote:
<br>>         >         >         Hello Dr Laino
<br>>         >         >         
<br>>         >         >         This is a problem with the ions. The atom
<br>>         >         >         indices it is warning about are all
<br>>         >         >         chlorine added to achieve electrical
<br>>         >         >         neutrality ( in VMD with Psfgen for NAMD)
<br>>         >         >         
<br>>         >         >         I do not understand because I've used the
<br>>         >         >         standard procedure from Psfgen
<br>>         >         >         
<br>>         >         >         On 11 Feb 2015 20:57, "Teodoro Laino"
<br>>         >         >         <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">teod...@gmail.com</a>> wrote:
<br>>         >         >                 double check your PSF or PDB.. the
<br>>         >         >                 error is self-esplicative : there
<br>>         >         >                 must be two or more molecules that
<br>>         >         >                 have the same name but have
<br>>         >         >                 different atom sequence..
<br>>         >         >                 
<br>>         >         >                 
<br>>         >         >                 the fact your PSF or PDB may work
<br>>         >         >                 with other codes means only that
<br>>         >         >                 such codes are not performing
<br>>         >         >                 tight checks.
<br>>         >         >                 
<br>>         >         >                 
<br>>         >         >                 Teo
<br>>         >         >                 
<br>>         >         >                 > On 11 Feb 2015, at 18:51, cjor
<br>>         >         >                 > <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">chr...@gmail.com</a>> wrote:
<br>>         >         >                 > 
<br>>         >         >                 > Hello, 
<br>>         >         >                 > I'm trying to run a QMMM
<br>>         >         >                 > simulation with CP2K based on a
<br>>         >         >                 > CHARMM FF MD simulation.
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > I first tried to run MM to see
<br>>         >         >                 > it works, but I get the
<br>>         >         >                 > following error:
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > "
<br>>         >         >                 >  CP2K| version string:
<br>>         >         >                 > 
<br>>         >         >                 >  CP2K version 2.3
<br>>         >         >                 >  CP2K| source code revision
<br>>         >         >                 > number:
<br>>         >         >                 >                12343
<br>>         >         >                 >  CP2K| is freely available from
<br>>         >         >                 > 
<br>>         >         >                 >   <a href="http://www.cp2k.org/" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.cp2k.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEfybep2Ix3JOTL9zjRCIb7jDt62g';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.cp2k.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEfybep2Ix3JOTL9zjRCIb7jDt62g';return true;">http://www.cp2k.org/</a>
<br>>         >         >                 >  CP2K| Program compiled at
<br>>         >         >                 >                      Thu Feb 28
<br>>         >         >                 > 16:12:09 GMT 2013
<br>>         >         >                 >  CP2K| Program compiled on
<br>>         >         >                 > 
<br>>         >         >                 >            popoca
<br>>         >         >                 >  CP2K| Program compiled for
<br>>         >         >                 > 
<br>>         >         >                 >  Linux-x86-64-gfortran
<br>>         >         >                 >  CP2K| Input file name
<br>>         >         >                 > 
<br>>         >         >                 >            mm.inp
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 >  GLOBAL| Force Environment
<br>>         >         >                 > number
<br>>         >         >                 >                      1
<br>>         >         >                 >  GLOBAL| Basis set file name
<br>>         >         >                 > BASIS_SET
<br>>         >         >                 >  GLOBAL| Geminal file name
<br>>         >         >                 > BASIS_GEMINAL
<br>>         >         >                 >  GLOBAL| Potential file name
<br>>         >         >                 > POTENTIAL
<br>>         >         >                 >  GLOBAL| MM Potential file name
<br>>         >         >                 > MM_POTENTIAL
<br>>         >         >                 >  GLOBAL| Coordinate file name
<br>>         >         >                 > 
<br>>         >         >                 >      ionized.pdb
<br>>         >         >                 >  GLOBAL| Method name
<br>>         >         >                 > 
<br>>         >         >                 >              CP2K
<br>>         >         >                 >  GLOBAL| Project name
<br>>         >         >                 > 
<br>>         >         >                 >          phd2_md
<br>>         >         >                 >  GLOBAL| Preferred FFT library
<br>>         >         >                 > FFTSG
<br>>         >         >                 >  GLOBAL| Run type
<br>>         >         >                 > MD
<br>>         >         >                 >  GLOBAL| All-to-all
<br>>         >         >                 > communication in single
<br>>         >         >                 > precision
<br>>         >         >                 >    F
<br>>         >         >                 >  GLOBAL| FFTs using library
<br>>         >         >                 > dependent lengths
<br>>         >         >                 >                    F
<br>>         >         >                 >  GLOBAL| Global print level
<br>>         >         >                 > 
<br>>         >         >                 >              LOW
<br>>         >         >                 >  GLOBAL| Total number of message
<br>>         >         >                 > passing processes
<br>>         >         >                 > 1
<br>>         >         >                 >  GLOBAL| Number of threads for
<br>>         >         >                 > this process
<br>>         >         >                 >                  1
<br>>         >         >                 >  GLOBAL| This output is from
<br>>         >         >                 > process
<br>>         >         >                 > 0
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 >  MEMORY| system memory details
<br>>         >         >                 > [Kb]
<br>>         >         >                 >  MEMORY|
<br>>         >         >                 >  rank 0           min
<br>>         >         >                 > max       average
<br>>         >         >                 >  MEMORY| MemTotal
<br>>         >         >                 >  132089732     132089732
<br>>         >         >                 > 132089732     132089732
<br>>         >         >                 >  MEMORY| MemFree
<br>>         >         >                 > 2537384       2537384
<br>>         >         >                 > 2537384       2537384
<br>>         >         >                 >  MEMORY| Buffers
<br>>         >         >                 >  407408        407408
<br>>         >         >                 >  407408        407408
<br>>         >         >                 >  MEMORY| Cached
<br>>         >         >                 > 34948436      34948436
<br>>         >         >                 >  34948436      34948436
<br>>         >         >                 >  MEMORY| Slab
<br>>         >         >                 > 732060        732060
<br>>         >         >                 >  732060        732060
<br>>         >         >                 >  MEMORY| SReclaimable
<br>>         >         >                 > 581768        581768
<br>>         >         >                 >  581768        581768
<br>>         >         >                 >  MEMORY| MemLikelyFree
<br>>         >         >                 >  38474996      38474996
<br>>         >         >                 >  38474996      38474996
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 >        49892           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 >        49893           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 >        49894           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 >        49895           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 >        49896           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 >        49897           0
<br>>         >         >                 >  Two molecules have been defined
<br>>         >         >                 > as identical molecules but atoms
<br>>         >         >                 > mismatch charges!!
<br>>         >         >                 > "
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > I'm calling a PSF and PDF file
<br>>         >         >                 > that I know works. The error
<br>>         >         >                 > lies with the parsing of the 
<br>>         >         >                 > ion section (SOD and CLA).
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > My MM input is:
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > &GLOBAL
<br>>         >         >                 >   PRINT_LEVEL LOW
<br>>         >         >                 >   PREFERRED_FFT_LIBRARY  FFTSG
<br>>         >         >                 >   PROJECT phd2_md
<br>>         >         >                 >   RUN_TYPE MD
<br>>         >         >                 > &END GLOBAL
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > &MOTION
<br>>         >         >                 >   &MD
<br>>         >         >                 >     ENSEMBLE   NPT_I
<br>>         >         >                 >     STEPS      6000
<br>>         >         >                 >     TIMESTEP   0.48
<br>>         >         >                 >     TEMPERATURE 298.0
<br>>         >         >                 >     &THERMOSTAT
<br>>         >         >                 >       TYPE NOSE
<br>>         >         >                 >       REGION MASSIVE
<br>>         >         >                 >       &NOSE
<br>>         >         >                 >         TIMECON   [wavenumber_t]
<br>>         >         >                 > 1000
<br>>         >         >                 >       &END NOSE
<br>>         >         >                 >     &END THERMOSTAT
<br>>         >         >                 >     &PRINT
<br>>         >         >                 >       &ENERGY
<br>>         >         >                 >         FILENAME =phd2_md.ener
<br>>         >         >                 >         &EACH
<br>>         >         >                 >           MD 1
<br>>         >         >                 >         &END EACH
<br>>         >         >                 >       &END ENERGY
<br>>         >         >                 >       &PROGRAM_RUN_INFO
<br>>         >         >                 >         &EACH
<br>>         >         >                 >           MD 1
<br>>         >         >                 >         &END EACH
<br>>         >         >                 >       &END PROGRAM_RUN_INFO
<br>>         >         >                 >     &END PRINT
<br>>         >         >                 >   &END MD
<br>>         >         >                 >   &PRINT
<br>>         >         >                 >     &TRAJECTORY
<br>>         >         >                 >       FILENAME =<a href="http://phd2_md.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fphd2_md.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGMn0T7gYW78mQ5S6rphrP-xrJEWQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fphd2_md.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGMn0T7gYW78mQ5S6rphrP-xrJEWQ';return true;">phd2_md.xyz</a>
<br>>         >         >                 >         &EACH
<br>>         >         >                 >           MD 1
<br>>         >         >                 >         &END EACH
<br>>         >         >                 >     &END TRAJECTORY
<br>>         >         >                 >     &RESTART
<br>>         >         >                 >       FILENAME =phd2_md.restart
<br>>         >         >                 >       BACKUP_COPIES 1
<br>>         >         >                 >       &EACH
<br>>         >         >                 >         MD 10
<br>>         >         >                 >       &END EACH
<br>>         >         >                 >     &END RESTART
<br>>         >         >                 >     &RESTART_HISTORY OFF
<br>>         >         >                 >     &END RESTART_HISTORY
<br>>         >         >                 >   &END PRINT
<br>>         >         >                 > &END MOTION
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > &FORCE_EVAL
<br>>         >         >                 >   METHOD FIST
<br>>         >         >                 >                 ! Using
<br>>         >         >                 > Molecular Mechanics
<br>>         >         >                 >   &MM
<br>>         >         >                 >     &FORCEFIELD
<br>>         >         >                 >       parm_file_name
<br>>         >         >                 > par_all27_prot_lipid.inp
<br>>         >         >                 >       parmtype CHM
<br>>         >         >                 >       ei_scale14 1.0
<br>>         >         >                 >       vdw_scale14 1.0
<br>>         >         >                 >       &SPLINE
<br>>         >         >                 >         emax_spline 1.0
<br>>         >         >                 >         rcut_nb 12
<br>>         >         >                 >       &END SPLINE
<br>>         >         >                 >     &END FORCEFIELD
<br>>         >         >                 >     &POISSON
<br>>         >         >                 >       &EWALD
<br>>         >         >                 >         ewald_type SPME
<br>>         >         >                 >         alpha 0.44
<br>>         >         >                 >         gmax 81
<br>>         >         >                 >       &END EWALD
<br>>         >         >                 >     &END POISSON
<br>>         >         >                 >   &END MM
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 >   &SUBSYS
<br>>         >         >                 >     &CELL
<br>>         >         >                 >       abc 80 80 80
<br>>         >         >                 >       periodic xyz
<br>>         >         >                 >     &END CELL
<br>>         >         >                 >     &KIND H
<br>>         >         >                 >       BASIS_SET DZVP-GTH-BLYP
<br>>         >         >                 >       POTENTIAL GTH-BLYP-q1
<br>>         >         >                 >     &END KIND
<br>>         >         >                 >     &KIND O
<br>>         >         >                 >       BASIS_SET DZVP-GTH-BLYP
<br>>         >         >                 >       POTENTIAL GTH-BLYP-q6
<br>>         >         >                 >     &END KIND
<br>>         >         >                 >     &KIND N
<br>>         >         >                 >       BASIS_SET DZVP-GTH-BLYP
<br>>         >         >                 >       POTENTIAL GTH-BLYP-q5
<br>>         >         >                 >     &END KIND
<br>>         >         >                 >     &KIND C
<br>>         >         >                 >       BASIS_SET DZVP-GTH-BLYP
<br>>         >         >                 >       POTENTIAL GTH-BLYP-q4
<br>>         >         >                 >     &END KIND
<br>>         >         >                 >     &KIND Fe2
<br>>         >         >                 >       BASIS_SET
<br>>         >         >                 > DZVP-MOLOPT-SR-GTH
<br>>         >         >                 >       POTENTIAL
<br>>         >         >                 > DZVP-MOLOPT-SR-GTH
<br>>         >         >                 >     &END KIND
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 >     &TOPOLOGY
<br>>         >         >                 >       CONNECTIVITY    UPSF
<br>>         >         >                 >       COORDINATE      PDB
<br>>         >         >                 >       COORD_FILE_NAME
<br>>         >         >                 > ionized.pdb
<br>>         >         >                 >       CONN_FILE_NAME
<br>>         >         >                 >  ionized.psf
<br>>         >         >                 >     &END TOPOLOGY
<br>>         >         >                 >   &END SUBSYS
<br>>         >         >                 > &END FORCE_EVAL
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > 
<br>>         >         >                 > -- 
<br>>         >         >                 > You received this message
<br>>         >         >                 > because you are subscribed to
<br>>         >         >                 > the Google Groups "cp2k" group.
<br>>         >         >                 > To unsubscribe from this group
<br>>         >         >                 > and stop receiving emails from
<br>>         >         >                 > it, send an email to cp2k
<br>>         >         >                 > +<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">unsu...@googlegroups.com</a>.
<br>>         >         >                 > To post to this group, send
<br>>         >         >                 > email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>>         >         >                 > Visit this group
<br>>         >         >                 > at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.
<br>>         >         >                 > For more options,
<br>>         >         >                 > visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>         >         >                 > 
<br>>         >         >                 
<br>>         >         >                 
<br>>         >         >                 
<br>>         >         >                 
<br>>         >         >                 -- 
<br>>         >         >                 You received this message because
<br>>         >         >                 you are subscribed to a topic in
<br>>         >         >                 the Google Groups "cp2k" group.
<br>>         >         >                 To unsubscribe from this topic,
<br>>         >         >                 visit <a href="https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe';return true;" onclick="this.href='https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe';return true;">https://groups.google.com/d/<wbr>topic/cp2k/kkfMLpVx2Xo/<wbr>unsubscribe</a>.
<br>>         >         >                 To unsubscribe from this group and
<br>>         >         >                 all its topics, send an email
<br>>         >         >                 to cp2k
<br>>         >         >                 +<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">unsu...@googlegroups.com</a>.
<br>>         >         >                 To post to this group, send email
<br>>         >         >                 to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>>         >         >                 Visit this group
<br>>         >         >                 at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.
<br>>         >         >                 For more options,
<br>>         >         >                 visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > 
<br>>         >         > -- 
<br>>         >         > _______
<br>>         >         > 
<br>>         >         > 
<br>>         >         > Christian Jørgensen
<br>>         >         > 
<br>>         >         > 
<br>>         >         > -- 
<br>>         >         > You received this message because you are
<br>>         >         > subscribed to the Google Groups "cp2k" group.
<br>>         >         > To unsubscribe from this group and stop receiving
<br>>         >         > emails from it, send an email to cp2k
<br>>         >         > +<a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">unsu...@googlegroups.com</a>.
<br>>         >         > To post to this group, send email
<br>>         >         > to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>>         >         > Visit this group
<br>>         >         > at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.
<br>>         >         > For more options,
<br>>         >         > visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>         >         
<br>>         >         
<br>>         >         
<br>>         >         
<br>>         >         -- 
<br>>         >         You received this message because you are subscribed
<br>>         >         to a topic in the Google Groups "cp2k" group.
<br>>         >         To unsubscribe from this topic, visit
<br>>         >         <a href="https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe';return true;" onclick="this.href='https://groups.google.com/d/topic/cp2k/kkfMLpVx2Xo/unsubscribe';return true;">https://groups.google.com/d/<wbr>topic/cp2k/kkfMLpVx2Xo/<wbr>unsubscribe</a>.
<br>>         >         To unsubscribe from this group and all its topics,
<br>>         >         send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br>>         >         To post to this group, send email to
<br>>         >         <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>>         >         Visit this group at
<br>>         >         <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.
<br>>         >         For more options, visit
<br>>         >         <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>         >         
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > 
<br>>         > -- 
<br>>         > _______
<br>>         > 
<br>>         > 
<br>>         > Christian Jørgensen
<br>>         > 
<br>>         > 
<br>>         > -- 
<br>>         > You received this message because you are subscribed to the
<br>>         > Google Groups "cp2k" group.
<br>>         > To unsubscribe from this group and stop receiving emails
<br>>         > from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br>>         > To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="gRPU4urLpFYJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>>         > Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.
<br>>         > For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>>         > 
<br>>         
<br>>         
<br>> 
<br>> 
<br>
<br></blockquote></div></div>