R: [gmx-users] implicit water and a layer of explicit water molecule

2011-04-04 Thread battis...@libero.it
Dear Mark,

about point 2, yes I need to have 
a uniform distribution of a defined numberof water molecule (eg. 100 water 
molecule ) into my box.
Is it possible with genbox?

After, I'll have to make the md simulation for my
system in implicit solvent
(I'll have  protein + 100 molecule SOL + implicit solvent)

So my next problem is to set the parameter into mdp file, for this mixed type 
of kind of water.

Thanks for your reply

Anna




 Dear all,


 I got some question about the implicit solvent.

 1)  In gromacs, is it possible simulate a protein in a layer  of 
 explicit water, and put this system (protein + SOL) into a big box and 
 make the MD simulation with implicit solvent?

 How I have to set the md.mdp parameter  (; IMPLICIT SOLVENT ALGORITHM) 
 in this case?

 Do you know some tutorial about this method?


It can't work as simply as that, because the waters on the edge will fly 
off into the implicit solvent region. People have tried various things - 
check out the literature.

 2) I'd like to put into my box a definied number of explicit number of 
 molecule of water eg. 100, so

 I used

 genbox -cp conf.gro -cs -maxsol 100  -p topol.top -o out.gro



 but the water molecule is not in the random position in the box, but 
 is in clustered conformation.

 Is it possible tell to genbox to put in a random way, in all space the 
 defined number of molecule?


What do you actually want - a uniform gas of a given density?

 3) For that system (100 explicit solvent molecule + implicit solvent) 
 I generated the topol.tpr using the following set up into the mdp file:


 ; IMPLICIT SOLVENT ALGORITHM
 implicit_solvent = GBSA
 ;GENERALIZED BORN ELECTROSTATICS;
 ;Algorithm for calculating Born radii
 gb_algorithm = Still


 grompp do not give problem at all,  but mdrun give problem:

 segmantation fault or

 Norm of force =nan

 I think that the problem is the use of explicit water molecule and the 
 implcit water together.


Maybe. We haven't got enough information to know.

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

Re: R: [gmx-users] implicit water and a layer of explicit water molecule

2011-04-04 Thread Mark Abraham

On 4/04/2011 7:12 PM, battis...@libero.it wrote:

Dear Mark,

about point 2, yes I need to have
a uniform distribution of a defined number
of water molecule (eg. 100 water molecule ) into my box.
Is it possible with genbox?


Yes, but not by starting with a uniform distribution of a 
condensed-phase density. You need one of the right density to start with.


A better approach is to decide how large a box of what density you want. 
Work out how much volume that gives to each molecule. Take a single 
molecule and put it in a box of that size with editconf. Move the 
molecule a bit off-center. Then use genconf -rot to replicate that box 
into a large one. Then equilibrate that thoroughly to get rid of the 
residual ordering. If, later on, you want a different size, genbox with 
the box you've equilibrated here will be a good approach.



After, I'll have to make the md simulation for my
system in implicit solvent
(I'll have  protein + 100 molecule SOL + implicit solvent)

So my next problem is to set the parameter into mdp file, for this mixed type 
of kind of water.


You have more problems than that. The force fields probably don't have 
implicit solvation parameters for water atom types. You'll need to 
source them somehow. And like I told you last time, your solvent 
molecules are not going to stay happily around your solvent like you 
hope. I'm going to stop repeating myself :-)


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


R: [gmx-users] implicit water and a layer of explicit water molecule

2011-04-04 Thread battis...@libero.it
Thank you very much for your suggestions!
Anna
-- 
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