[gmx-users] How to use "tabulate" to add new potential function

2012-02-21 Thread smith bill
I want to add new potential function into Gromacs, and I know there is a way to do this by using "tabulate". But I am not sure about how to use tabulate,and how many new potential function could be added? -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/

Re: [gmx-users] Re: Placing ions in the specified positions

2012-02-21 Thread James Starlight
Dear all :) I have one extra question about Genion. I want to neitralise my system with the qtot - 1.550002e by placing some Na and Cl ions under the physiological concetrantion 100 mmol/l I've performed genion -s ions.tpr -o b2ar_ions.gro -p topol.top -pname NA -nname CL -conc 0.1 -neutral bu

[gmx-users] modifications required for the ffnonbonded.itp while adding a new residue

2012-02-21 Thread ramesh cheerla
Dear Gromacs users, i am adding new residue to the existing charmm27 force field of the gromacs and i am planning to use buckingham potential for the non bonded interactions. for that i made following changes to the ffnonbonded.itp file [ atomtypes ] ;name at.num m

Re: [gmx-users] About the removal of the center of mass

2012-02-21 Thread xiaojing gong
Many thanks, can you tell me the reason? becuase I see in the old manual, there said "if there is fixed part in the system, then you can not use the removal of the center of the mass." XJ 2012/2/21, Yang Ye : > yes. you need to and can still put System for the removal group. > >> >> On Feb 21, 201

RE: [gmx-users] interface

2012-02-21 Thread Dallas Warren
What does it look like when you visualise the coordinate file? That is how you can answer that question, as we certainly can't. And as you will soon find out, you cannot have a single interface in a simulation box, you have to have two. Catch ya, Dr. Dallas Warren Medicinal Chemistry and Drug

RE: [gmx-users] air-water interface

2012-02-21 Thread Dallas Warren
Search the literature and see what others have done to do exactly this. Catch ya, Dr. Dallas Warren Medicinal Chemistry and Drug Action Monash Institute of Pharmaceutical Sciences, Monash University 381 Royal Parade, Parkville VIC 3010 dallas.war...@monash.edu +61 3 9903 9304

[gmx-users] Distance Restraints on Protein - possible at all?

2012-02-21 Thread NG HUI WEN
From: NG HUI WEN Sent: Sunday, February 19, 2012 1:19 PM To: gmx-users@gromacs.org Subject: Distance Restraints on Protein - possible at all? Dear gmxusers, I have been trying to apply distance restraints on my protein but have been unsuccessful thus far. I have consulted the user forum

Re: [gmx-users] Pressure coupling and membrane-type simulations

2012-02-21 Thread Peter C. Lai
On 2012-02-21 03:41:07PM -0500, Justin A. Lemkul wrote: > > > Andrew DeYoung wrote: > > Hi, > > > > I am interested in doing a membrane-type simulation, in which I have > > all-atom membrane "walls" parallel to xy plane, at z = -z_0 and z = +z_0 > > (where z_0 is a constant). I would like to ru

Re: [gmx-users] COMPASS force field revisited Bond Dihedral and Angle Dihedral cross terms

2012-02-21 Thread Broadbent, Richard
On 21/02/2012 19:02, "David van der Spoel" wrote: > On 2/21/12 6:17 PM, Richard Broadbent wrote: >> Dear All, >> >> I am considering conducting a simulation of a polymeric system in >> gromacs. I would like to use the COMPASS forcefield as it has a complete >> parameter set for my molecule. >

Re: [gmx-users] interface

2012-02-21 Thread Mark Abraham
On 22/02/2012 3:11 AM, mohammad agha wrote: Dear Gromacs Specialists, I made a box consists of water with box lengths: 6nm * 6nm * 6nm , then I equilibrated it with NPT ensemble, box size increased to 6.66176, then I kept the /x-/ and /y-/dimensions fixed, and double the system size in /z

Re: [gmx-users] Pressure coupling and membrane-type simulations

2012-02-21 Thread Justin A. Lemkul
Andrew DeYoung wrote: Hi, I am interested in doing a membrane-type simulation, in which I have all-atom membrane "walls" parallel to xy plane, at z = -z_0 and z = +z_0 (where z_0 is a constant). I would like to run an NPT simulation at 1 atm. What type of pressure coupling should I use? Is

Re: [gmx-users] several questions about g_hbond -contact

2012-02-21 Thread Justin A. Lemkul
Anna Marabotti wrote: Dear users, I'm using g_hbond with option -contact in order to find the contacts between my protein and a series of ligands, in a radius of 0.5 nm. I made several calculations but I'm quite uncertain about the results, therefore I'm asking you some questions. The comm

[gmx-users] several questions about g_hbond -contact

2012-02-21 Thread Anna Marabotti
Dear users, I'm using g_hbond with option -contact in order to find the contacts between my protein and a series of ligands, in a radius of 0.5 nm. I made several calculations but I'm quite uncertain about the results, therefore I'm asking you some questions. The command I used is: g_hbond -f my_

Re: [gmx-users] COMPASS force field revisited Bond Dihedral and Angle Dihedral cross terms

2012-02-21 Thread David van der Spoel
On 2/21/12 6:17 PM, Richard Broadbent wrote: Dear All, I am considering conducting a simulation of a polymeric system in gromacs. I would like to use the COMPASS forcefield as it has a complete parameter set for my molecule. I believe the majority of the implementation is simple (though long an

[gmx-users] COMPASS force field revisited Bond Dihedral and Angle Dihedral cross terms

2012-02-21 Thread Richard Broadbent
Dear All, I am considering conducting a simulation of a polymeric system in gromacs. I would like to use the COMPASS forcefield as it has a complete parameter set for my molecule. I believe the majority of the implementation is simple (though long and fiddly). However, it has Bond-Dihedral and An

[gmx-users] Pressure coupling and membrane-type simulations

2012-02-21 Thread Andrew DeYoung
Hi, I am interested in doing a membrane-type simulation, in which I have all-atom membrane "walls" parallel to xy plane, at z = -z_0 and z = +z_0 (where z_0 is a constant). I would like to run an NPT simulation at 1 atm. What type of pressure coupling should I use? Isotropic pressure coupling

[gmx-users] Regarding Cluster g_cluster

2012-02-21 Thread Singam Karthick
Dear all, >I have done the clustering analysis using g_cluster, I got the 11 clusters i >have used the cutoff of 0.4 nm, now i want to see how well the >structure distributed in each clusters. I have got the middle structure of >each clusters, Here i have given the cutoff of 0.4 nm, due to which

[gmx-users] interface

2012-02-21 Thread mohammad agha
Dear Gromacs Specialists, I made a box consists of water with box lengths:  6nm  * 6nm  * 6nm , then I equilibrated it with NPT ensemble, box size increased to 6.66176, then I kept the x- and y-dimensions fixed, and double the system size in z as following: editconf -f pr1.gro -o newbox1.gro -bo

Re: Sv: Re: [gmx-users] double and single precision

2012-02-21 Thread Edvin Erdtman
Hi GMX-users! If someone else get the same problem. I will present how I could recover the edr files: 1. Calculate the number of bytes for header consisted of in single and double precision, respectively. run some short test runs to calculate this 2. Calculate the number of bytes for each fram

[gmx-users] tpbconv - segmantation fault

2012-02-21 Thread יוכבד
Hi I'm trying to set to zero the charges of my ligand group (MOL) using tpbconv: tpbconv_d -s md.tpr -f md.trr -n -o md_lig_zero.tpr -zeroq then I select the MOL group and I get the following message: Selected 14: 'MOL' Segmentation fault (core dumped) Any explanation to why this may have happ

[gmx-users] Using tpbconv_d -zeroq

2012-02-21 Thread יוכבד
Hi I ran a md simulation for a system with two groups: Surrounding , LIG. I wish to set the charges of the surrounding group to zero and to recalculate the Coulomb interactions between the groups. For this purpose I used tpbconv with the -zeroq option: tpbconv_d -s md_prev.tpr -f md_prev.trr -n

Re: [gmx-users] About the removal of the center of mass

2012-02-21 Thread Yang Ye
yes. you need to and can still put System for the removal group. > > On Feb 21, 2012 5:51 PM, "xiaojing gong" wrote: >> >> Dear all, >> I used the GMX 4.5.5 to run a simulation with a fixed plane and water. >> I wonder to know if I fix the plane, shall I use the command of >> the remov

Re: [gmx-users] how to extract the x y z

2012-02-21 Thread Justin A. Lemkul
lina wrote: Hi, I don't know how to extract the position like the x y z with t or r with t. g_traj prints coordinates from a trajectory. -Justin -- Justin A. Lemkul Ph.D. Candidate ICTAS Doctoral Scholar MILES-IGERT Trainee Department of Biochemist

Re: [gmx-users] Umbrella_pull_simulation

2012-02-21 Thread Justin A. Lemkul
shahid nayeem wrote: Initially I used g_wham -if pullf-files.dat -it tpr-files.dat -unit Kcal -histo. But now as suggested by you I added -b 1000 -e 1 leaving 1ns for equilibriation. The new profile.xvg is attached. How can I further improve it. I don't know. You've never even state

[gmx-users] how to extract the x y z

2012-02-21 Thread lina
Hi, I don't know how to extract the position like the x y z with t or r with t. Thanks for any suggestions, Best regards, -- 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_L

Re: [gmx-users] Correcting PBC effects

2012-02-21 Thread Tsjerk Wassenaar
Hi Ariel, If it works, it's good :) It really only affects the visualization... It doesn't change the system. Cheers, Tsjerk On Tue, Feb 21, 2012 at 10:33 AM, Siew Wen Leong wrote: > Dear all, > > Have anyone tried using CHARMM36 equilibrated POPE bilayer by Dr.Klauda > http://terpconnect.umd.

[gmx-users] About the removal of the center of mass

2012-02-21 Thread xiaojing gong
Dear all, I used the GMX 4.5.5 to run a simulation with a fixed plane and water. I wonder to know if I fix the plane, shall I use the command of the removal of the center of mass for water? I will appreciate it if you will share your opion. Many thanks XJ -- gmx-users mailing list

[gmx-users] Correcting PBC effects

2012-02-21 Thread Siew Wen Leong
Dear all, Have anyone tried using CHARMM36 equilibrated POPE bilayer by Dr.Klauda http://terpconnect.umd.edu/~jbklauda/research/download.html? I am planning to use the membrane patch for a transmembrane protein simulation. At first visualization, the membrane patch consists of a bilayer with a si

[gmx-users] air-water interface

2012-02-21 Thread mohammad agha
Dear Gromacs Specialists, I want to make a system consists of one box of solvent (for example water) that there is 100 surfactant molecules in it and water is in contact with air in surface as same as we do in laboratory in beaker. But I don't know how should I define this air-water interface f