<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>