Re: [gmx-users] Re: Polarizable halide (Iodide)

2011-03-02 Thread David van der Spoel

On 2011-03-02 09.44, Mikhail Stukan wrote:

David,

What should not be working in the current version?
Do you mean that any polarizable model (water for example) can run on a single 
processor only?
I am currently using the model by Roux et al. for NaCl solution in SWM4-NDP 
model and everything seems to be working fine.
I also had some problems with the collapsing of the system which seems to be 
solved by setting explicitly all pair interactions in [nonbond_param] statement 
(I had a feeling that mixing rule did not work for some reason). The results 
are very reasonable.
Does it mean that I am doing something completely wrong?


I have an open issue on this:

http://redmine.gromacs.org/issues/713

maybe you can download the test files (based on the same model) that I 
put there and check them. Let's continue the discussion in the redmine 
forum.



Many thanks in advance,
Mikhail


On 2011-03-01 22.06, Abhijeet Joshi wrote:

Hi all,

 I am trying implement polarizable model for halides
suggested by Roux group.
I am not able to equilibrate such system.


In present gromacs versions this runs on 1 processor (core) only.



System: Sodium Iodide in water
I am using following topology file for iodide.

Anybody has any clue what might be happening?

Question 2: I am having trouble finding implementation details for
thole_polarization.
How would that appear in topology file?



;
; Topology file for iodide polarizable

[ moleculetype ]
; molnamenrexcl
I1

[ atoms ]
; idat typeres nr residu nameat namecg nrcharge
1   II 1  SICI11   3.733085
2   IS1  SISI1-4.733085


[ polarization ]
; See notes above.alpha (nm^3)
121 0.007439



[ exclusions ]
; iatom excluded from interaction with i
12
21

#ifdef POSRES
; Restrain the oxygen...
[ position_restraints ]
; iatom typefxfyfz
11100100100
#endif



Your help will be greatly appreciated

Thanks
~ Abhijeet




--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell&  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se





--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Re: Polarizable halide (Iodide)

2011-03-02 Thread Mikhail Stukan
David,

What should not be working in the current version?
Do you mean that any polarizable model (water for example) can run on a single 
processor only?
I am currently using the model by Roux et al. for NaCl solution in SWM4-NDP 
model and everything seems to be working fine.
I also had some problems with the collapsing of the system which seems to be 
solved by setting explicitly all pair interactions in [nonbond_param] statement 
(I had a feeling that mixing rule did not work for some reason). The results 
are very reasonable.
Does it mean that I am doing something completely wrong?

Many thanks in advance,
Mikhail


On 2011-03-01 22.06, Abhijeet Joshi wrote:
> Hi all,
>
> I am trying implement polarizable model for halides
> suggested by Roux group.
> I am not able to equilibrate such system.

In present gromacs versions this runs on 1 processor (core) only.

>
> System: Sodium Iodide in water
> I am using following topology file for iodide.
>
> Anybody has any clue what might be happening?
>
> Question 2: I am having trouble finding implementation details for
> thole_polarization.
> How would that appear in topology file?
>
>
>
> ;
> ; Topology file for iodide polarizable
>
> [ moleculetype ]
> ; molnamenrexcl
> I1
>
> [ atoms ]
> ; idat typeres nr residu nameat namecg nrcharge
> 1   II 1  SICI11   3.733085
> 2   IS1  SISI1-4.733085
>
>
> [ polarization ]
> ; See notes above.alpha (nm^3)
> 121 0.007439
>
>
>
> [ exclusions ]
> ; iatom excluded from interaction with i
> 12
> 21
>
> #ifdef POSRES
> ; Restrain the oxygen...
> [ position_restraints ]
> ; iatom typefxfyfz
> 11100100100
> #endif
>
>
>
> Your help will be greatly appreciated
>
> Thanks
> ~ Abhijeet
>


--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists