Respectable David/Justin/Mark/Berk
In my study of diffusion of oxygen in water, I solved the segmentation fault 
problem as below
1. At first I ran for 16 ps beyond which segmentation fault occured in the 
previous run. I generated the .gro file of that run
2. I then ran for 100 ps using the .gro file of step 1. Also generated the .gro 
file from this step. I constantly checked the pressure content. It really 
decreased from 4000's to 1400s. No segmentation fault occurred.
3 I than ran again for 500 ps using the same technique.
4. When I ran for 1 ns, there was no segmentation fault as well as the pressure 
was about 1 bar which I needed. The temperature was 298 K as I needed.
 
Now the run upto 100 ns went without any segmentation fault.
Is this the way I should solve the problem?
 
Now when I saw the msd curve, it was nearly straight line and the calculation 
yielded the diffusion coefficient of oxygen as 40000 x 10(-9), where as the 
experimental value is nearly 2 x 10(-9).
Actually, as the time increases it increases, Why is this happening?
We are doubtful about the force field parameters we have supplied. 
Specially, the nonbonded parameters.
 
I have tried to define bare oxygen in .atp file 
Is not that the charge of bare oxygen is zero?
 
Could you tell me what is the value for nonbonded parameters for oxygen 
molecule and spc water molecule? 
OB--------OW
OB---------HW
 
In the nb.itp the nonbonded parameter (Linard-Jonnes) between oxygen and water 
hydrogen is given as zero.
I am stuck at this step and I need to solve it anyway for my Msc. Thesis. 
Please help me.
 
 


      
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php

Reply via email to