Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Sunita
Hello Venkat,

Ya, I am sure that I use 1st.cpt to run next 40 ns job
 Find below the output of trjcat

Summary of files and start times used:

  FileStart time   Time step
-
complex_prod1_old.xtc0.000 ps   50.000 ps
complex_prod1.xtc0.000 ps   50.000 ps WARNING: same
Start time as previous
 
Thanks

--
View this message in context: 
http://gromacs.5086.x6.nabble.com/Fwd-Concatenating-two-trajectories-tp5015956p5015962.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Venkat Reddy
So, that is the problem. Your second trajectory also has the same start
time. Some how you didn't extend your second run from where your 1st run
has finished.

Out of curiosity, why is your time step 50ps? Its too big for a simulation.
Generally it should be 2fs.


On Mon, Apr 21, 2014 at 12:30 PM, Sunita  wrote:

> Hello Venkat,
>
> Ya, I am sure that I use 1st.cpt to run next 40 ns job
>  Find below the output of trjcat
>
> Summary of files and start times used:
>
>   FileStart time   Time step
> -
> complex_prod1_old.xtc0.000 ps   50.000 ps
> complex_prod1.xtc0.000 ps   50.000 ps WARNING: same
> Start time as previous
>
> Thanks
>
> --
> View this message in context:
> http://gromacs.5086.x6.nabble.com/Fwd-Concatenating-two-trajectories-tp5015956p5015962.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>



-- 
With Best Wishes
Venkat Reddy Chirasani
PhD student
Laboratory of Computational Biophysics
Department of Biotechnology
IIT Madras
Chennai
INDIA-600036
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] System coordinates not getting freezed

2014-04-21 Thread Mark Abraham
You're just trying to measure a potential energy, so do it the right way
and don't get involved wiht making mistakes trying to freeze atoms. See
http://www.gromacs.org/Documentation/How-tos/Single-Point_Energy

Mark


On Mon, Apr 21, 2014 at 7:38 AM, sukriti002  wrote:

> Hi everyone,
>
> My system contains one copper ion and one water molecule with no solvent. I
> considered complete system as one group and freezed the system. But the
> coordinates are not getting fixed. Also when i took the same mdp file with
> 2
> copper ions instead of one cu and one water, the system was getting frozen.
> Can anyone please let me know what is the problem.
>
> Regards
> Sukriti
>
> --
> View this message in context:
> http://gromacs.5086.x6.nabble.com/System-coordinates-not-getting-freezed-tp5015955.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Concatenating two trajectories

2014-04-21 Thread Monoj Mon Kalita
Dear Sunita

Are you sure you have added two trajectory files, one from 0 to 40 and the
other from 40 onward ! Because

Here







*Summary of files and start times used:  FileStart
time   Time
step-
complex_prod1_old.xtc0.000 ps   50.000 ps
complex_prod1.xtc0.000 ps   50.000 ps WARNING: sameStart time
as previous*


it is saying that you are joining the same two files with same starting
time. However, do you have .trr files, GROMACS generally creates .trr files
for every successful run. You can try to add those .trr files and later
convert them into .xtc files. And you just need to follow like this

trjcat -f 1.trr 2.trr -o all.trr

That is it.

In another case, you said that you are getting strange output in RMSD. Did
you check, is there any broken molecule or not ! sometime the broken
molecule, due to periodicity, leads to strange RMSD outputs.

Thanks

Monoj
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Position Restraints with NVT and NPT

2014-04-21 Thread Mark Abraham
The restraints just push on the atoms. If other stuff pushes harder, then
the restraints don't work. More relevant is whether you think your starting
position models something with accuracy down to a resolution of 0.01nm, or
some aspect of the starting position is sensitive to such a difference.
(Answer: if it came from an experiment or a simulation, no)

Mark


On Mon, Apr 21, 2014 at 8:41 AM, zanemarkson  wrote:

> Thank you Justin!
>
> I understand usually for regular equilibration, restraints should be
> enough.
> But I am still wondering to what extent position restraint can restrain my
> atoms. In fact, I ran a relatively longer NVT using position restraints and
> I saw that both the absolute coordinates of my "restrained" atoms and the
> distances between my "restrained" atoms changed by ~0.01nm and I think that
> is not negligible any more. I don't know whether this is normal or a
> problem
> that I should somehow fix. Any thoughts? Thank you!
>
> --
> View this message in context:
> http://gromacs.5086.x6.nabble.com/Position-Restraints-with-NVT-and-NPT-tp5015953p5015961.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Mark Abraham
If you want to share a file, please use a file sharing service and share
the link.

Mark


On Mon, Apr 21, 2014 at 7:32 AM, sunita gupta  wrote:

> Hello Everyone,
>
>
> I have two trajectories of 40 ns each, the first one is from 0 to 40 ns and
> next one is extended to next 40 ns taking .cpt file of the first
>  I concatenating the two using two commands:
> 1. trjcat -f 1st.xtc 2nd.xtc -o combine1.xtc -cat
> 2. trjcat -f 1st.xtc 2nd.xtc -o combine2.xtc -settime (selecting 0.000 ps
> for 1st and 4 ps for 2nd)
>
> But when I plotted the rmsd using combine1.xtc and combine2.xtc I am
> getting very strange plots
>
> Unfortunately I am not able to post the plots as its saying msg is too big
> to be posted and bouncing back...
> Can anyone suggest me whats wrong with this method.
>
>
>
>
>
> Best Regards
> --
> SUNITA GUPTA
> Senior Research Fellow
> Bioinformatics Centre
> Jawaharlal Nehru University
> New Delhi- 110067
> Email- sunita@gmail.com
>
>
>
>
>
> --
> --
> SUNITA GUPTA
> Senior Research Fellow
> Bioinformatics Centre
> Jawaharlal Nehru University
> New Delhi- 110067
> Email- sunita@gmail.com
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Tsjerk Wassenaar
Hi Venkat,

That time step is the time step in the output trajectory. Yours said 2.000
ps, not fs.

Cheers,

Tsjerk
On Apr 21, 2014 9:35 AM, "Venkat Reddy"  wrote:

> So, that is the problem. Your second trajectory also has the same start
> time. Some how you didn't extend your second run from where your 1st run
> has finished.
>
> Out of curiosity, why is your time step 50ps? Its too big for a simulation.
> Generally it should be 2fs.
>
>
> On Mon, Apr 21, 2014 at 12:30 PM, Sunita  wrote:
>
> > Hello Venkat,
> >
> > Ya, I am sure that I use 1st.cpt to run next 40 ns job
> >  Find below the output of trjcat
> >
> > Summary of files and start times used:
> >
> >   FileStart time   Time step
> > -
> > complex_prod1_old.xtc0.000 ps   50.000 ps
> > complex_prod1.xtc0.000 ps   50.000 ps WARNING: same
> > Start time as previous
> >
> > Thanks
> >
> > --
> > View this message in context:
> >
> http://gromacs.5086.x6.nabble.com/Fwd-Concatenating-two-trajectories-tp5015956p5015962.html
> > Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> > --
> > Gromacs Users mailing list
> >
> > * Please search the archive at
> > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> > posting!
> >
> > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> >
> > * For (un)subscribe requests visit
> > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> > send a mail to gmx-users-requ...@gromacs.org.
> >
>
>
>
> --
> With Best Wishes
> Venkat Reddy Chirasani
> PhD student
> Laboratory of Computational Biophysics
> Department of Biotechnology
> IIT Madras
> Chennai
> INDIA-600036
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Venkat Reddy
Hi Tsjerk,
Sorry for the confusion. I thought its time step which is actually
output writing frequency. Thank you sir for the correction.


On Mon, Apr 21, 2014 at 2:39 PM, Tsjerk Wassenaar  wrote:

> Hi Venkat,
>
> That time step is the time step in the output trajectory. Yours said 2.000
> ps, not fs.
>
> Cheers,
>
> Tsjerk
> On Apr 21, 2014 9:35 AM, "Venkat Reddy"  wrote:
>
> > So, that is the problem. Your second trajectory also has the same start
> > time. Some how you didn't extend your second run from where your 1st run
> > has finished.
> >
> > Out of curiosity, why is your time step 50ps? Its too big for a
> simulation.
> > Generally it should be 2fs.
> >
> >
> > On Mon, Apr 21, 2014 at 12:30 PM, Sunita  wrote:
> >
> > > Hello Venkat,
> > >
> > > Ya, I am sure that I use 1st.cpt to run next 40 ns job
> > >  Find below the output of trjcat
> > >
> > > Summary of files and start times used:
> > >
> > >   FileStart time   Time step
> > > -
> > > complex_prod1_old.xtc0.000 ps   50.000 ps
> > > complex_prod1.xtc0.000 ps   50.000 ps WARNING: same
> > > Start time as previous
> > >
> > > Thanks
> > >
> > > --
> > > View this message in context:
> > >
> >
> http://gromacs.5086.x6.nabble.com/Fwd-Concatenating-two-trajectories-tp5015956p5015962.html
> > > Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> > > --
> > > Gromacs Users mailing list
> > >
> > > * Please search the archive at
> > > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> > > posting!
> > >
> > > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> > >
> > > * For (un)subscribe requests visit
> > > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> > > send a mail to gmx-users-requ...@gromacs.org.
> > >
> >
> >
> >
> > --
> > With Best Wishes
> > Venkat Reddy Chirasani
> > PhD student
> > Laboratory of Computational Biophysics
> > Department of Biotechnology
> > IIT Madras
> > Chennai
> > INDIA-600036
> > --
> > Gromacs Users mailing list
> >
> > * Please search the archive at
> > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> > posting!
> >
> > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> >
> > * For (un)subscribe requests visit
> > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> > send a mail to gmx-users-requ...@gromacs.org.
> >
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>



-- 
With Best Wishes
Venkat Reddy Chirasani
PhD student
Laboratory of Computational Biophysics
Department of Biotechnology
IIT Madras
Chennai
INDIA-600036
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] all-angles

2014-04-21 Thread Maria Astón Serrano
Dear Gromacs users

I would like to know whether you have ever run a simulation with
constraints in "all-angles", or if you know about somebody who did it.
When I try to run a protein with this option, it crashes, and only
artificial systems like linear chains or modified small molecules seem
to work with all-angles.

Best regards,

Maria
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] how to couple a organic molecule in water

2014-04-21 Thread xiao
Dear gromacs users,


I am doing a simulation of an organic molecule in water solution. The organic 
molecule has 137 atoms (including hydrogen atoms). I have no idea on how to do 
temperature coupling. I made two temperature groups, one is water, and the 
other is the organic molecule. I do not  know whether it is reasonable. Another 
way is to  couple the organic molecule with the water molecules.  Which way is 
better?
Any information on that is appreciated.


Best regards
Fugui
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] g_potential help

2014-04-21 Thread rajat desikan
Hi All,

I have a membrane-protein-solvent system, and I want the potential across
the membrane in the direction of the membrane normal. I have 2 questions in
that regard.

1) We know that epsilon_r = 80 for the water, ~7 for the membrane-water
interface and ~2 for the hydrophobic membrane core. g_potential returns the
result for epsilon_r = 1. How can I incorporate this into the result?

2) Default ng =1? Is the number of groups referring to the number of charge
groups defined in the topology?

Thank you for your time.

-- 
Rajat Desikan (Ph.D Scholar)
Prof. K. Ganapathy Ayappa's Lab (no 13),
Dept. of Chemical Engineering,
Indian Institute of Science, Bangalore
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Cannot write trajectory frame; maybe you are out of disk space

2014-04-21 Thread Andrew Bostick
Dear gromacs users

I am doing md simulation on the a protein with 190 residues.

I did minimization and equilibration without problem and error.

But in production run step (last md simulation with 1500 steps),
in step 650, I encountered with following error:

File input/output error:
Cannot write trajectory frame; maybe you are out of disk space.

my mdp file is as follows:
-
title= opls Protein MD
; Run parameters
integrator= md; leap-frog integrator
nsteps= 1500; 0.002 * 50 = 1000 ps, 30 ns
dt = 0.002; 2 fs
; Output control
nstxout= 3000 ; save coordinates every 2 ps
nstvout= 3000 ; save velocities every 2 ps
nstxtcout= 3000 ; xtc compressed trajectory output
every 2 ps
nstenergy= 3000 ; save energies every 2 ps
nstlog= 3000 ; update log file every 2 ps
energygrps  = Protein Ion
; Bond parameters
continuation= yes; Restarting after NPT
constraint_algorithm = lincs; holonomic constraints
constraints= all-bonds; all bonds (even heavy atom-H
bonds) constrained
lincs_iter= 1; accuracy of LINCS
lincs_order= 4; also related to accuracy
; Neighborsearching
ns_type= grid; search neighboring grid cells
nstlist= 5; 10 fs
rlist= 1.0 ; short-range neighborlist cutoff
(in nm)
rcoulomb= 1.0 ; short-range electrostatic cutoff
(in nm)
rvdw= 1.0 ; short-range van der Waals cutoff
(in nm)
; Electrostatics
coulombtype= PME; Particle Mesh Ewald for
long-range electrostatics
pme_order= 4; cubic interpolation
fourierspacing= 0.16; grid spacing for FFT
; Temperature coupling is on
tcoupl= V-rescale; modified Berendsen thermostat
tc-grps= Protein Non-Protein; two coupling groups - more
accurate
tau_t= 0.10.1; time constant, in ps
ref_t= 300 300; reference temperature, one for
each group, in K
; Pressure coupling is on
pcoupl= Parrinello-Rahman; Pressure coupling on in NPT
pcoupltype= isotropic; uniform scaling of box vectors
tau_p= 2.0; time constant, in ps
ref_p= 1.0; reference pressure, in bar
compressibility = 4.5e-5; isothermal compressibility of water,
bar^-1
; Periodic boundary conditions
pbc= xyz; 3-D PBC
; Dispersion correction
DispCorr= EnerPres; account for cut-off vdW scheme
; Velocity generation
gen_vel= no; Velocity generation is off
-

If I use df -h command to know disk space, I obtained following information:

Filesystem  Size  Used Avail Use% Mounted on
/dev/sda3   220G   14G  196G   7% /
none4.0K 0  4.0K   0% /sys/fs/cgroup
udev3.0G  4.0K  3.0G   1% /dev
tmpfs   607M  1.3M  606M   1% /run
none5.0M 0  5.0M   0% /run/lock
none3.0G  152K  3.0G   1% /run/shm
none100M   64K  100M   1% /run/user
/dev/sda6   245G   61G  184G  25% /media/pdfco/my dear
/dev/sda2   220G   85G  136G  39% /media/pdfco/Windows
/dev/sda5   245G   43G  202G  18% /media/pdfco/nazif
--

Is my computer system suitable for this md simulation?

How to solve this error?

Any help will highly appreciated.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Cannot write trajectory frame; maybe you are out of disk space

2014-04-21 Thread Mark Abraham
If you can rule out file permissions issues (because you wrote files
earlier), and actually running out of disk, usually that would suggest a
flaky network file system.

Mark


On Mon, Apr 21, 2014 at 4:09 PM, Andrew Bostick
wrote:

> Dear gromacs users
>
> I am doing md simulation on the a protein with 190 residues.
>
> I did minimization and equilibration without problem and error.
>
> But in production run step (last md simulation with 1500 steps),
> in step 650, I encountered with following error:
>
> File input/output error:
> Cannot write trajectory frame; maybe you are out of disk space.
>
> my mdp file is as follows:
>
> -
> title= opls Protein MD
> ; Run parameters
> integrator= md; leap-frog integrator
> nsteps= 1500; 0.002 * 50 = 1000 ps, 30 ns
> dt = 0.002; 2 fs
> ; Output control
> nstxout= 3000 ; save coordinates every 2 ps
> nstvout= 3000 ; save velocities every 2 ps
> nstxtcout= 3000 ; xtc compressed trajectory output
> every 2 ps
> nstenergy= 3000 ; save energies every 2 ps
> nstlog= 3000 ; update log file every 2 ps
> energygrps  = Protein Ion
> ; Bond parameters
> continuation= yes; Restarting after NPT
> constraint_algorithm = lincs; holonomic constraints
> constraints= all-bonds; all bonds (even heavy atom-H
> bonds) constrained
> lincs_iter= 1; accuracy of LINCS
> lincs_order= 4; also related to accuracy
> ; Neighborsearching
> ns_type= grid; search neighboring grid cells
> nstlist= 5; 10 fs
> rlist= 1.0 ; short-range neighborlist cutoff
> (in nm)
> rcoulomb= 1.0 ; short-range electrostatic cutoff
> (in nm)
> rvdw= 1.0 ; short-range van der Waals cutoff
> (in nm)
> ; Electrostatics
> coulombtype= PME; Particle Mesh Ewald for
> long-range electrostatics
> pme_order= 4; cubic interpolation
> fourierspacing= 0.16; grid spacing for FFT
> ; Temperature coupling is on
> tcoupl= V-rescale; modified Berendsen thermostat
> tc-grps= Protein Non-Protein; two coupling groups - more
> accurate
> tau_t= 0.10.1; time constant, in ps
> ref_t= 300 300; reference temperature, one for
> each group, in K
> ; Pressure coupling is on
> pcoupl= Parrinello-Rahman; Pressure coupling on in NPT
> pcoupltype= isotropic; uniform scaling of box vectors
> tau_p= 2.0; time constant, in ps
> ref_p= 1.0; reference pressure, in bar
> compressibility = 4.5e-5; isothermal compressibility of water,
> bar^-1
> ; Periodic boundary conditions
> pbc= xyz; 3-D PBC
> ; Dispersion correction
> DispCorr= EnerPres; account for cut-off vdW scheme
> ; Velocity generation
> gen_vel= no; Velocity generation is off
>
> -
>
> If I use df -h command to know disk space, I obtained following
> information:
>
> Filesystem  Size  Used Avail Use% Mounted on
> /dev/sda3   220G   14G  196G   7% /
> none4.0K 0  4.0K   0% /sys/fs/cgroup
> udev3.0G  4.0K  3.0G   1% /dev
> tmpfs   607M  1.3M  606M   1% /run
> none5.0M 0  5.0M   0% /run/lock
> none3.0G  152K  3.0G   1% /run/shm
> none100M   64K  100M   1% /run/user
> /dev/sda6   245G   61G  184G  25% /media/pdfco/my dear
> /dev/sda2   220G   85G  136G  39% /media/pdfco/Windows
> /dev/sda5   245G   43G  202G  18% /media/pdfco/nazif
>
> --
>
> Is my computer system suitable for this md simulation?
>
> How to solve this error?
>
> Any help will highly appreciated.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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


[gmx-users] Brownian Dynamics

2014-04-21 Thread Yutian Yang
Dear all, 

I noticed that Gromacs has the option to do Brownian Dynamics, is there anyone 
who has done that or where can I find any tutorials that can teach me to do it? 

Thanks. 
Sincerely
Yutian (Shirley) Yang
Biomedical and Chemical Engineering
329 Link Hall Syracuse University 
Syracuse NY 13244
yyan...@syr.edu



-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] H-bond analysis

2014-04-21 Thread Erik Marklund
Hi,

I've never used it myself, hut hbplus allegedly uses a more complex hb 
criterion. http://www.jhu.edu/fleming/Modeling/HBPLUS.html

Kind regards,
Erik

On 17 Apr 2014, at 10:38, michael.b  wrote:

> 
> hi,
> 
> I try to analyze H-bonds in an MD trajectory.
> when i look at the definition of an H-bond used by g_hbond it appears
> that only three atoms (acceptor, donor, hydrogen, ADH) are used
> in the analysis performed by this program ... however, i would assume that
> in fact 4 (or up to 6) atoms should be considered because the existence and
> the strength of an H-bond not only depend on the A-D distance and
> the ADH angle  but also on the orientation of the lone pairs of the
> acceptor, and the latter are actually determined by the geometry of
> the atom(s), call them X, the acceptor is covalently bound to, aren't they?
> 
> So, in short, my question is: is anybody aware of a tool that
> can analyze the H-bonds in MD trajectories considering not only
> A-D distance and ADH angle but also the position of X (basically
> the dihedral DHAX) ?
> 
> thanks for any hints
> cheers,
> michael
> 
> 
> 
> 
> --
> View this message in context: 
> http://gromacs.5086.x6.nabble.com/H-bond-analysis-tp5015899.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> -- 
> Gromacs Users mailing list
> 
> * Please search the archive at 
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> 
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> 
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
> mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] negative HB lifetime

2014-04-21 Thread Erik Marklund
Hi,

Yes. Or use larger analysis groups if possible.

Kind regards,
Erik

On 17 Apr 2014, at 10:58, Djamila  wrote:

> Thanks Erik;
> You're alright, The problem is in the ACF; it falls dramatically to zero.
> In this case, what's the solution? should I extend my simulation to more
> time?
> Thanks
> 
> Djamila
> 
> --
> View this message in context: 
> http://gromacs.5086.x6.nabble.com/negative-HB-lifetime-tp5015830p5015900.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> -- 
> Gromacs Users mailing list
> 
> * Please search the archive at 
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> 
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> 
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
> mail to gmx-users-requ...@gromacs.org.

-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] How to fix or constrain Center of Mass

2014-04-21 Thread Tom
Dear Gromacs Users,

Is there a way to constrain or fix Center of Mass
of lipid bilayer or proteins in a MD run?

Thanks!

Thomas
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Cannot write trajectory frame; maybe you are out of disk space

2014-04-21 Thread Andrew Bostick
Dear Mark

Thanks for your reply.


I am beginner in linux and gromacs.

How to rule out file permissions issues? Which command is true for doing
this?

Please guide me to solve this problem.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] all-angles

2014-04-21 Thread Justin Lemkul



On 4/21/14, 6:21 AM, Maria Astón Serrano wrote:

Dear Gromacs users

I would like to know whether you have ever run a simulation with
constraints in "all-angles", or if you know about somebody who did it.
When I try to run a protein with this option, it crashes, and only
artificial systems like linear chains or modified small molecules seem
to work with all-angles.



None of the force fields provided with Gromacs were parametrized with fully 
rigid angles, so using that option probably just breaks the model physics.  I 
would expect such an option to only be useful for very small molecules that can 
be sensibly represented as rigid.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Cannot write trajectory frame; maybe you are out of disk space

2014-04-21 Thread Mark Abraham
Did earlier frames get written? Use ls -l, gmxcheck etc. Google about file
permissions - unlike on (say) an old Windows machine, on a compute cluster
you will not be able to write files anywhere but where you're allowed to.
Guessing here, because you haven't supplied any context.

Mark


On Mon, Apr 21, 2014 at 7:27 PM, Andrew Bostick
wrote:

> Dear Mark
>
> Thanks for your reply.
>
>
> I am beginner in linux and gromacs.
>
> How to rule out file permissions issues? Which command is true for doing
> this?
>
> Please guide me to solve this problem.
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] how to couple a organic molecule in water

2014-04-21 Thread Justin Lemkul



On 4/21/14, 7:53 AM, xiao wrote:

Dear gromacs users,


I am doing a simulation of an organic molecule in water solution. The organic 
molecule has 137 atoms (including hydrogen atoms). I have no idea on how to do 
temperature coupling. I made two temperature groups, one is water, and the 
other is the organic molecule. I do not  know whether it is reasonable. Another 
way is to  couple the organic molecule with the water molecules.  Which way is 
better?
Any information on that is appreciated.



For a small molecule in water, a single thermostat to the whole system seems to 
be the most sensible approach to me.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] g_potential help

2014-04-21 Thread Justin Lemkul



On 4/21/14, 9:58 AM, rajat desikan wrote:

Hi All,

I have a membrane-protein-solvent system, and I want the potential across
the membrane in the direction of the membrane normal. I have 2 questions in
that regard.

1) We know that epsilon_r = 80 for the water, ~7 for the membrane-water
interface and ~2 for the hydrophobic membrane core. g_potential returns the
result for epsilon_r = 1. How can I incorporate this into the result?



I'd venture a guess that the value of epsilon_r in the output is the same as the 
interpretation of epsilon_r in the .mdp file - it is the relative dielectric 
permittivity.  A value of 1 means don't scale anything; see previous discussions.



2) Default ng =1? Is the number of groups referring to the number of charge
groups defined in the topology?



No, that means the number of groups you want to analyze, just like g_rdf and 
other tools sometimes give the option of doing multiple analyses at the same time.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] How to fix or constrain Center of Mass

2014-04-21 Thread Justin Lemkul



On 4/21/14, 12:23 PM, Tom wrote:

Dear Gromacs Users,

Is there a way to constrain or fix Center of Mass
of lipid bilayer or proteins in a MD run?



To fix the position of the COM, use nstcomm = 1 (but your performance may suffer 
badly), otherwise removal of COM drift periodically is handled by any value of 
nstcomm > 0.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] [gmx-developers] Regarding Amber topology and co-ordinate file

2014-04-21 Thread Justin Lemkul


Please ask usage questions on the gmx-users mailing list.  The development list 
is for coding and planning discussions; this post is not appropriate here.  I am 
CC'ing this message over to gmx-users.  Please post anything further there.


On 4/21/14, 10:28 AM, MOHD HOMAIDUR RAHMAN wrote:

Dear Gromacs User

I have ILs system in small box that is running on Amber 12 MD package. Now I
want to use Gromacs MD package for my bigger system (Large no of molecules with
proteins) with the same parameter.

So, it is possible to convert Amber topology format to gromacs topology format.
I have also prep and frcmod file for that molecules. which is use in tleap to
generate topology file and Co-ordinate file to direct use in Amber MD package.



There are various scripts floating around the web that do such conversions; 
amb2gmx.pl comes to mind.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Using a forcefield from NAMD

2014-04-21 Thread Ali Khan
Hi,

I am interested in incorporating a nitroxide spin label into my protein
simulation. I want to use the Charmm36 force field, and to my knowledge
there is no force field for the MTSL spin label in charmm36 for gromacs.
What I do have is a charmm22 force field of the MSTL label from NAMD. Is
there a way to incorporate this NAMD force field into my gromacs
simulations.

Best Wishes,
Ali Khan
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] how to couple a organic molecule in water

2014-04-21 Thread xiao
Hi Justin,


Thank you very much for your information


Best regards


Fugui








At 2014-04-22 01:52:09,"Justin Lemkul"  wrote:
>
>
>On 4/21/14, 7:53 AM, xiao wrote:
>> Dear gromacs users,
>>
>>
>> I am doing a simulation of an organic molecule in water solution. The 
>> organic molecule has 137 atoms (including hydrogen atoms). I have no idea on 
>> how to do temperature coupling. I made two temperature groups, one is water, 
>> and the other is the organic molecule. I do not  know whether it is 
>> reasonable. Another way is to  couple the organic molecule with the water 
>> molecules.  Which way is better?
>> Any information on that is appreciated.
>>
>
>For a small molecule in water, a single thermostat to the whole system seems 
>to 
>be the most sensible approach to me.
>
>-Justin
>
>-- 
>==
>
>Justin A. Lemkul, Ph.D.
>Ruth L. Kirschstein NRSA Postdoctoral Fellow
>
>Department of Pharmaceutical Sciences
>School of Pharmacy
>Health Sciences Facility II, Room 601
>University of Maryland, Baltimore
>20 Penn St.
>Baltimore, MD 21201
>
>jalem...@outerbanks.umaryland.edu | (410) 706-7441
>http://mackerell.umaryland.edu/~jalemkul
>
>==
>-- 
>Gromacs Users mailing list
>
>* Please search the archive at 
>http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
>
>* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
>* For (un)subscribe requests visit
>https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
>mail to gmx-users-requ...@gromacs.org.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] g_potential help

2014-04-21 Thread rajat desikan
Hi Justin,

So different parts of my system in Z-direction should be scaled
differently, right? On a coarse level, aqueous part with a scaling of 80,
membrane with a scaling of 2? How do I correct for this? Can I take the
quantitative output for potential difference across a membrane in water
from g_potential literally?

Thanks,

On Monday, April 21, 2014, Justin Lemkul  wrote:

>
>
> On 4/21/14, 9:58 AM, rajat desikan wrote:
>
>> Hi All,
>>
>> I have a membrane-protein-solvent system, and I want the potential across
>> the membrane in the direction of the membrane normal. I have 2 questions
>> in
>> that regard.
>>
>> 1) We know that epsilon_r = 80 for the water, ~7 for the membrane-water
>> interface and ~2 for the hydrophobic membrane core. g_potential returns
>> the
>> result for epsilon_r = 1. How can I incorporate this into the result?
>>
>>
> I'd venture a guess that the value of epsilon_r in the output is the same
> as the interpretation of epsilon_r in the .mdp file - it is the relative
> dielectric permittivity.  A value of 1 means don't scale anything; see
> previous discussions.
>
>  2) Default ng =1? Is the number of groups referring to the number of
>> charge
>> groups defined in the topology?
>>
>>
> No, that means the number of groups you want to analyze, just like g_rdf
> and other tools sometimes give the option of doing multiple analyses at the
> same time.
>
> -Justin
>
> --
> ==
>
> Justin A. Lemkul, Ph.D.
> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>
> Department of Pharmaceutical Sciences
> School of Pharmacy
> Health Sciences Facility II, Room 601
> University of Maryland, Baltimore
> 20 Penn St.
> Baltimore, MD 21201
>
> jalem...@outerbanks.umaryland.edu | (410) 706-7441
> http://mackerell.umaryland.edu/~jalemkul
>
> ==
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/
> Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>


-- 
Rajat Desikan (Ph.D Scholar)
Prof. K. Ganapathy Ayappa's Lab (no 13),
Dept. of Chemical Engineering,
Indian Institute of Science, Bangalore
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Using a forcefield from NAMD

2014-04-21 Thread Justin Lemkul



On 4/21/14, 3:02 PM, Ali Khan wrote:

Hi,

I am interested in incorporating a nitroxide spin label into my protein
simulation. I want to use the Charmm36 force field, and to my knowledge
there is no force field for the MTSL spin label in charmm36 for gromacs.
What I do have is a charmm22 force field of the MSTL label from NAMD. Is
there a way to incorporate this NAMD force field into my gromacs
simulations.



You've got parameters in one format, and you need them in another, so all you 
need to do is convert them.  If you've got a stream file that you're reading 
into NAMD, our Python script should do the job for you.  Look for 
cgenff_charmm2gmx.py:


http://mackerell.umaryland.edu/CHARMM_ff_params.html

-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] g_potential help

2014-04-21 Thread Justin Lemkul



On 4/21/14, 5:11 PM, rajat desikan wrote:

Hi Justin,

So different parts of my system in Z-direction should be scaled
differently, right? On a coarse level, aqueous part with a scaling of 80,
membrane with a scaling of 2? How do I correct for this? Can I take the
quantitative output for potential difference across a membrane in water
from g_potential literally?



You don't need to scale anything.  The force field parameters for every 
biomolecular force field that I know of are based on epsilon_r = 1, i.e. no 
special alteration of vacuum permittivity.  The results, at face value, should 
be sufficient.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] calculation speed of GROMACS and cut off parameters

2014-04-21 Thread mircial
Hi everyone,

I am using GROMACS to simulate bilayers, and I am testing two sets of 
parameters at present. However, I found significant differences of simulation 
speed of these two parameters sets. 
The mainly differences between these parameters are the cut off parameters or 
of the VDW and coulomb interactions, in the first parameter set, a cut off of 
1.0 nm was used for both of the VDW and the coulomb interactions. In the second 
parameter set, a cut off 1.4 was used for the coulomb interactions, whereas a 
switch from 1.4 nm to 1.5 nm was used for the VDW parameters.  Under the first 
parameter set, my system could run about 20 ns per day, whereas for the second 
parameter set, my system only run 7-8 ns per day.  So I am wondering if this is 
only raised because of the cut off parameters? And is this a normal pheonomen?

I am listing the mdp file at the end of this email. Please note the only 
relevarent differences between these two parameter sets are the cut off 
differences. 

With my best regards, 
rxgu

the first parameter set:

title  = NPT Equilibration


integrator = md
nsteps = 2500
dt = 0.002
nstxout= 0  
nstvout= 0  
nstenergy  = 1000  
nstlog = 2000  
nstxtcout  = 2000
xtc-precision  = 1

constraint_algorithm = lincs  
constraints= all-bonds
lincs_iter = 1  
lincs_order= 4 

;cutoff-scheme  = Verlet  
ns_type= grid  
nstlist= 10  
rlist   = 1.0
  
coulombtype= PME
rcoulomb   = 1.0
  

  
vdwtype   = cut-off 
  
rvdw   = 1.0
  
DispCorr   = EnerPres   
  
pme_order  = 4  
  
fourierspacing = 0.12   
  
optimize_fft   = yes
  

  
pbc= xyz
  

  
tcoupl = V-rescale  
  
tc-grps= POPC SOL   
  
tau_t   = 0.5   0.5 
  
ref_t   = 323  323

pcoupl = Berendsen  
pcoupltype = semiisotropic  
tau_p   = 10.0 10.0  
ref_p   = 1.0  1.0  
compressibility = 4.5e-5 4.5e-5  

gen_vel = yes
gen_temp= 323
gen_seed= 82907543

comm-mode   = Linear
comm-grps = system

the second parameter set:

title  = NPT Equilibration


integrator = md
nsteps = 2500
dt = 0.002
nstxout= 0  
nstvout= 0  
nstenergy  = 1000  
nstlog = 2000  
nstxtcout  = 2000
xtc-precision  = 1

constraint_algorithm = lincs  
constraints= all-bonds
lincs_iter = 1  
lincs_order= 4 

;cutoff-scheme  = Verlet  
ns_type= grid  
nstlist= 10  
rlist   = 1.4
rlistlong  = 1.7
  
coulombtype= PME
rcoulomb   = 1.4 
 
vdwtype   = Switch
rvdw   = 1.5
rvdw-switch= 1.4
DispCorr   = EnerPres
pme_order  = 4  
fourierspacing = 0.12  
optimize_fft   = yes

pbc= xyz

tcoupl = V-rescale  
tc-grps   

Re: [gmx-users] calculation speed of GROMACS and cut off parameters

2014-04-21 Thread Justin Lemkul



On 4/21/14, 5:48 PM, mirc...@sjtu.edu.cn wrote:

Hi everyone,

I am using GROMACS to simulate bilayers, and I am testing two sets of 
parameters at present. However, I found significant differences of simulation 
speed of these two parameters sets.
The mainly differences between these parameters are the cut off parameters or 
of the VDW and coulomb interactions, in the first parameter set, a cut off of 
1.0 nm was used for both of the VDW and the coulomb interactions. In the second 
parameter set, a cut off 1.4 was used for the coulomb interactions, whereas a 
switch from 1.4 nm to 1.5 nm was used for the VDW parameters.  Under the first 
parameter set, my system could run about 20 ns per day, whereas for the second 
parameter set, my system only run 7-8 ns per day.  So I am wondering if this is 
only raised because of the cut off parameters? And is this a normal pheonomen?



Yes, it's normal.  The most important thing to know is you need to use the 
parameters for which the force field was designed.  Lipids are extremely 
sensitive to changes in the cutoff.  Varying these values is very hazardous, so 
if you're doing something non-standard, you'd better have a lot of good proof 
that what you're doing isn't going to lead to artifacts.


-Justin


I am listing the mdp file at the end of this email. Please note the only 
relevarent differences between these two parameter sets are the cut off 
differences.

With my best regards,
rxgu

the first parameter set:

title  = NPT Equilibration


integrator = md
nsteps = 2500
dt = 0.002
nstxout= 0
nstvout= 0
nstenergy  = 1000
nstlog = 2000
nstxtcout  = 2000
xtc-precision  = 1

constraint_algorithm = lincs
constraints= all-bonds
lincs_iter = 1
lincs_order= 4

;cutoff-scheme  = Verlet
ns_type= grid
nstlist= 10
rlist   = 1.0

coulombtype= PME
rcoulomb   = 1.0

vdwtype   = cut-off
rvdw   = 1.0
DispCorr   = EnerPres
pme_order  = 4
fourierspacing = 0.12
optimize_fft   = yes

pbc= xyz

tcoupl = V-rescale
tc-grps= POPC SOL
tau_t   = 0.5   0.5
ref_t   = 323  323

pcoupl = Berendsen
pcoupltype = semiisotropic
tau_p   = 10.0 10.0
ref_p   = 1.0  1.0
compressibility = 4.5e-5 4.5e-5

gen_vel = yes
gen_temp= 323
gen_seed= 82907543

comm-mode   = Linear
comm-grps = system

the second parameter set:

title  = NPT Equilibration


integrator = md
nsteps = 2500
dt = 0.002
nstxout= 0
nstvout= 0
nstenergy  = 1000
nstlog = 2000
nstxtcout  = 2000
xtc-precision  = 1

constraint_algorithm = lincs
constraints= all-bonds
lincs_iter = 1
lincs_order= 4

;cutoff-scheme  = Verlet
ns_type= grid
nstlist= 10
rlist   = 1.4
rlistlong  = 1.7

coulombtype= PME
rcoulomb   = 1.4

vdwtype   = Switch
rvdw   = 1.5
rvdw-switch= 1.4
DispCorr   = EnerPres
pme_order  = 4
fourierspacing = 0.12
optimize_fft   = yes

pbc= xyz

tcoupl = V-rescale
tc-grps= POPC SOL
tau_t   = 0.5   0.5
ref_t   = 323  323

pcoupl = Berendsen
pcoupltype = semiisotropic
tau_p   = 10.0 10.0
ref_p   = 1.0  1.0
compressibility = 4.5e-5 4.5e-5

gen_vel = yes
gen_temp= 323
gen_seed= 82907543

comm-mode   = Linear
comm-grps   = System
--



--
==

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] g_potential help

2014-04-21 Thread rajat desikan
Thank you!

On Tuesday, April 22, 2014, Justin Lemkul  wrote:

>
>
> On 4/21/14, 5:11 PM, rajat desikan wrote:
>
>> Hi Justin,
>>
>> So different parts of my system in Z-direction should be scaled
>> differently, right? On a coarse level, aqueous part with a scaling of 80,
>> membrane with a scaling of 2? How do I correct for this? Can I take the
>> quantitative output for potential difference across a membrane in water
>> from g_potential literally?
>>
>>
> You don't need to scale anything.  The force field parameters for every
> biomolecular force field that I know of are based on epsilon_r = 1, i.e. no
> special alteration of vacuum permittivity.  The results, at face value,
> should be sufficient.
>
> -Justin
>
> --
> ==
>
> Justin A. Lemkul, Ph.D.
> Ruth L. Kirschstein NRSA Postdoctoral Fellow
>
> Department of Pharmaceutical Sciences
> School of Pharmacy
> Health Sciences Facility II, Room 601
> University of Maryland, Baltimore
> 20 Penn St.
> Baltimore, MD 21201
>
> jalem...@outerbanks.umaryland.edu | (410) 706-7441
> http://mackerell.umaryland.edu/~jalemkul
>
> ==
> --
> Gromacs Users mailing list
>
> * Please search the archive at http://www.gromacs.org/
> Support/Mailing_Lists/GMX-Users_List before posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-requ...@gromacs.org.
>


-- 
Rajat Desikan (Ph.D Scholar)
Prof. K. Ganapathy Ayappa's Lab (no 13),
Dept. of Chemical Engineering,
Indian Institute of Science, Bangalore
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] Cannot write trajectory frame; maybe you are out of disk space

2014-04-21 Thread Andrew Bostick
Dear Mark

Thanks for your reply.

I googled file permission in linux OS. After that I did following steps:

I created a directory entitled input.

I put pdb and mdp files in the input directory.

I used following command to change permission:

chmod -R 777 input

Now I will md simulation steps (pdb2gmx  editconf ..)

Is my manner true?
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


Re: [gmx-users] Fwd: Concatenating two trajectories

2014-04-21 Thread Sunita
Hello Venkat,

Yeah, I think by default it will give 0 ps as start time for both
trajectories and it should be set by user using -settime flag
Correct me if I am wrong.

Regards

--
View this message in context: 
http://gromacs.5086.x6.nabble.com/Fwd-Concatenating-two-trajectories-tp5015956p5015992.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.


[gmx-users] No default U-B types

2014-04-21 Thread Shima Arasteh
Hi all,

There is an unnatural sidechain in my input file, so I added new lines as 
new-defined(DPP) residue to the .rtp and .hdb files.
Also there were some missed bonds, angles, and dihedrals in generated .top 
file, which manually added in.
When I go through the grompp to do EM, I find a fatal error about the new 
defined residue:

No default U-B types

As I checked this kind of angle, it exists in other residues such as lysine as 
well. This angle is (defined by atom types): 
CT1 CT2 NH3 


I checked the atom names, It seems ok.

Part of the top file is as follow:

 residue   1 DPP rtp DPP  q +2.0
 1    NH3  1    DPP  N  1   -0.3 14.007   ; qtot 
-0.3
 2 HC  1    DPP H1  2   0.33  1.008   ; qtot 
0.03
 3 HC  1    DPP H2  3   0.33  1.008   ; qtot 
0.36
 4 HC  1    DPP H3  4   0.33  1.008   ; qtot 
0.69
 5    CT1  1    DPP CA  5   0.21 12.011   ; qtot 0.9
 6 HB  1    DPP HA  6    0.1  1.008   ; qtot 1
 7    CT2  1    DPP CB  7   0.21 12.011   ; qtot 
1.21
 8 HA  1    DPP    HB1  8   0.05  1.008   ; qtot 
1.26
 9 HA  1    DPP    HB2  9   0.05  1.008   ; qtot 
1.31
    10    NH3  1    DPP NG 10   -0.3 14.007   ; qtot 
1.01
    11 HC  1    DPP    HG1 11   0.33  1.008   ; qtot 
1.34
    12 HC  1    DPP    HG2 12   0.33  1.008   ; qtot 
1.67
    13 HC  1    DPP    HG3 13   0.33  1.008   ; qtot 2
    14  C  1    DPP  C 14   0.51 12.011   ; qtot 
2.51
    15  O  1    DPP  O 15  -0.51 15.999   ; qtot 2

Fatal error:
ERROR 1 [file topol.top, line 853]:
  No default U-B types

line 853:     5 7    10 5 



Would you please let me your suggestions about the reason of this error? How 
can I pass it? Would you please help me?

Thanks in advance for your suggestions.

Sincerely,
Shima 
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

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

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.