[gmx-users] Re: Bilayer COM removal issue: Large VCM

2013-11-12 Thread rajat desikan
Hi All, Any suggestions? Thanks, On Mon, Nov 11, 2013 at 12:38 AM, rajat desikan wrote: > Hi All, > I am experiencing a few problems in membrane simulations wrt COM removal. > I downloaded a 400 ns pre-equilibrated Slipid-DMPC membrane with all the > accompanying files. I then carried out the f

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Dear Justin and Piggot, Thanks for the suggestions. Actually, I have constructed a CG lipid vesicle by placing lipids in random conformation in a simulation box. My lipid system is heterogeneous, i.e., it has different types of lipids (POPC,CHOLESTEROL,PPC...etc). Some of them occupy the outer laye

[gmx-users] help about ibi

2013-11-12 Thread guozhicheng222
Hi: When I am running the ibi procedure, I get the following error message: A coordinate in file conf.gro does not contain a '.' Additionally, I check the coordinate file of confout.gro in step_001. It showed that 'nan' symbol appeared in confout

Re: [gmx-users] Re: Change in the position of structural Zinc and calcium ions during MD

2013-11-12 Thread Justin Lemkul
On 11/12/13 1:47 PM, Rama wrote: Hi Justin, Below I pasted .mdp file and topology. In .log file I could see energy term for position restraints. .mdp file--- title = NPT Equilibration define = -DPOSRES ; position restraints for protein ; Run paramet

Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)

2013-11-12 Thread Dwey Kauffman
Hi Mark and Szilard Thanks for your both suggestions. They are very helpful. > > Neither run had a PP-PME work distribution suitable for the hardware it > was > running on (and fixing that for each run requires opposite changes). > Adding > a GPU and hoping to see scaling requires that there

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Thomas Piggot
Hi Venkat, Can you make it a bit clearer what you actually want? If it is the diffusion of the lipids along the curved surface of the vesicle, rather than simply the overall 3D diffusion, this is not trivial to calculate as I don't believe g_msd will do this for you. This property has been st

[gmx-users] Re: Change in the position of structural Zinc and calcium ions during MD

2013-11-12 Thread Rama
Hi Justin, Below I pasted .mdp file and topology. In .log file I could see energy term for position restraints. .mdp file--- title = NPT Equilibration define = -DPOSRES ; position restraints for protein ; Run parameters integrator = md

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Justin Lemkul
On 11/12/13 1:33 PM, Venkat Reddy wrote: Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? If you just want the overall diffusion constant, that's what g_msd does without any additional options. -Justin On Tue,

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Thanks Justin. So, I have to calculate diffusion coefficient three times (x,y,z) and finally add-up together to get in 3D??? On Tue, Nov 12, 2013 at 11:25 PM, Justin Lemkul wrote: > > > On 11/12/13 12:30 PM, Venkat Reddy wrote: > >> Dear Sir, Thanks for the quick reply. >> So, I have to declare

Re: [gmx-users] Change in the positon of structural Zinc and calcium ions during MD

2013-11-12 Thread Justin Lemkul
On 11/12/13 12:35 PM, Rama wrote: Hi Gromacs users, I'm doing protein-Bilayer MD simulations. Enzyme contains structural zinc and calcium ions during Energy minimization, NVT and NPT stage, ions are changing there position even though I applied position restraints for the atoms and ions. Any

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Justin Lemkul
On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? The options for the -type flag are x, y, or z. You said you wanted the diffusion coefficient in each spatial dimension. That is precisely what this option will d

[gmx-users] Change in the positon of structural Zinc and calcium ions during MD

2013-11-12 Thread Rama
Hi Gromacs users, I'm doing protein-Bilayer MD simulations. Enzyme contains structural zinc and calcium ions during Energy minimization, NVT and NPT stage, ions are changing there position even though I applied position restraints for the atoms and ions. Anyone could help me out. Thanks -- Vi

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Dear Sir, Thanks for the quick reply. So, I have to declare -type no flag. Isn't it?? and I have recently gone through Justin's membrane protein tutorial, where he has calculated diffusion coefficient for lipids in a membrane by creating an index group for a particular atom. So, here also shall I d

Re: [gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun

2013-11-12 Thread Justin Lemkul
On 11/12/13 12:14 PM, cjalmeciga wrote: The output of energy minimization was Potential Energy = -1.42173622068236e+06 Maximum force = 9.00312066109319e+02 on atom 148 Norm of force = 2.06087515037187e+01 OK, reasonable enough. How about a description of what the system is, whi

[gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun

2013-11-12 Thread cjalmeciga
The output of energy minimization was Potential Energy = -1.42173622068236e+06 Maximum force = 9.00312066109319e+02 on atom 148 Norm of force = 2.06087515037187e+01 Thanks Javier Justin Lemkul wrote > On 11/12/13 10:58 AM, cjalmeciga wrote: >> I run >> >> grompp -f nvt.mdp -c em.gr

Re: [gmx-users] Re: Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread Justin Lemkul
On 11/12/13 12:07 PM, arunjones wrote: Thank you Sir, initially i was running on 60 threads, now i changed it to 100. simulation is running with out any error, but i found a note in the log file as fallows #nodes mismatch, current program: 100 checkpoint file: 60 #PME-nodes mi

[gmx-users] Re: Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread arunjones
Thank you Sir, initially i was running on 60 threads, now i changed it to 100. simulation is running with out any error, but i found a note in the log file as fallows #nodes mismatch, current program: 100 checkpoint file: 60 #PME-nodes mismatch, current program: -1 checkpoint

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Justin Lemkul
On 11/12/13 11:25 AM, Venkat Reddy wrote: Then, how to mention the direction for spherical particles Sir? Read g_msd -h again, paying specific attention to the -type flag. -Justin On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul wrote: On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Then, how to mention the direction for spherical particles Sir? On Tue, Nov 12, 2013 at 7:28 PM, Justin Lemkul wrote: > > > On 11/12/13 8:55 AM, Venkat Reddy wrote: > >> Thank you sir for the prompt reply. >> *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* >> >> Here I

Re: [gmx-users] Re: Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread Justin Lemkul
On 11/12/13 11:10 AM, arunjones wrote: Hello Sir, Thanks for the reply. now is it fine if i use 100 threads in my restart. is there any impact on the over all simulation? Only if that is the number of threads originally used in the run. If not, there will be a mismatch between the DD grid

Re: [gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun

2013-11-12 Thread Justin Lemkul
On 11/12/13 10:58 AM, cjalmeciga wrote: I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. The fact that grompp completes indicates there is nothing syntactically wrong with the input files. Whethe

[gmx-users] Re: Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread arunjones
Hello Sir, Thanks for the reply. now is it fine if i use 100 threads in my restart. is there any impact on the over all simulation? -- View this message in context: http://gromacs.5086.x6.nabble.com/Restarting-a-simulation-after-replacing-an-empty-md-trr-file-tp5012443p5012459.html Sent from the

[gmx-users] Re: segmentation fault on gromacs 4.5.5 after mdrun

2013-11-12 Thread cjalmeciga
I run grompp -f nvt.mdp -c em.gro -p topol.top -n index.ndx -o nvt.tpr and everything looks fine. I check the nvt.tpr, and temperature is ok. the real problem is with the mdrun function. could be a problem of the software? Thanks Javier Justin Lemkul wrote > On 11/11/13 11:24 AM, Carlos

Re: [gmx-users] Invalid order for directive defaults

2013-11-12 Thread Justin Lemkul
On 11/12/13 10:07 AM, Atila Petrosian wrote: Dear all My system contains protein + cnt + water molecules. I have summarized what I did below: --- 1) By pdb2gmx and charmm27 force field, I obtained pr.top for prot

Re: mdrun on 8-core AMD + GTX TITAN (was: Re: [gmx-users] Re: Gromacs-4.6 on two Titans GPUs)

2013-11-12 Thread Szilárd Páll
As Mark said, please share the *entire* log file. Among other important things, the result of PP-PME tuning is not included above. However, I suspect that in this case scaling is strongly affected or by the small size of the system you are simulating. -- Szilárd On Sun, Nov 10, 2013 at 5:28 AM,

[gmx-users] Invalid order for directive defaults

2013-11-12 Thread Atila Petrosian
Dear all My system contains protein + cnt + water molecules. I have summarized what I did below: --- 1) By pdb2gmx and charmm27 force field, I obtained pr.top for protein then I converted it into pr.itp by deleting

[gmx-users] Re: ok, thank you

2013-11-12 Thread Williams Ernesto Miranda Delgado
> Send gmx-users mailing list submissions to > gmx-users@gromacs.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.gromacs.org/mailman/listinfo/gmx-users > or, via email, send a message with subject or body 'help' to > gmx-users-requ...@gromacs.org >

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Justin Lemkul
On 11/12/13 8:55 AM, Venkat Reddy wrote: Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? No. The system is a sphere, so what use

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Thank you sir for the prompt reply. *g_msd -f traj.xtc -s topol.tpr -n msd.ndx -lateral z -o msd.xvg -tu ns* Here I am giving -lateral z (like for membrane simulations). Is it fine for spherical systems also? On Tue, Nov 12, 2013 at 5:32 PM, Dr. Vitaly Chaban wrote: > MSD is 3D by default. > > >

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread Justin Lemkul
On 11/12/13 8:33 AM, bharat gupta wrote: Thanks justin for your replies. I understood the g_analyze related data. I tired g_analyze to dump the structures as you said. But, I didn't find any switch that can be used to dump the structure in pdb format. Because that's not the function of g_ana

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread bharat gupta
Thanks justin for your replies. I understood the g_analyze related data. I tired g_analyze to dump the structures as you said. But, I didn't find any switch that can be used to dump the structure in pdb format. On Tue, Nov 12, 2013 at 10:15 PM, Justin Lemkul wrote: > > > On 11/11/13 6:56 PM, bh

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread Justin Lemkul
On 11/11/13 6:56 PM, bharat gupta wrote: Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ?

Re: [gmx-users] Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread Justin Lemkul
On 11/12/13 1:15 AM, arun kumar wrote: Dear Gromacs users, I am running a 50ns simulation of a protein having nearly 700 residues on 60 threads (Gromacs 4.6.3). At one point i got a disk space problem, so i have deleted the md.trr file and created an empty md.trr file. when i tried to restart

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread Justin Lemkul
On 11/11/13 5:39 PM, bharat gupta wrote: Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.

Re: [gmx-users] Re: Reaction field zero and ions

2013-11-12 Thread Dr. Vitaly Chaban
There are no problems to have ions while using Reaction-Field treatment. Dr. Vitaly V. Chaban On Mon, Nov 11, 2013 at 7:06 PM, Justin Lemkul wrote: > > > On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: >> >> Hello >> If I did the MD simulation using PME and neutralized with ions,

Re: [gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Dr. Vitaly Chaban
MSD is 3D by default. Dr. Vitaly V. Chaban On Tue, Nov 12, 2013 at 6:01 AM, Venkat Reddy wrote: > Dear all, > I am simulating a spherical lipid vesicle. I want to calculate the > diffusion coefficient for each lipid component in 3D. How to calculate it > using g_msd (or any other tool like g_v

Re: [gmx-users] Re: Reaction field zero and ions

2013-11-12 Thread Justin Lemkul
On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: Hello If I did the MD simulation using PME and neutralized with ions, and I want to rerun this time with reaction field zero, is there any problem if I keep the ions? This is for LIE calculation. I am using AMBER99SB. Why do you th

[gmx-users] Restarting a simulation after replacing an empty md.trr file

2013-11-12 Thread arun kumar
Dear Gromacs users, I am running a 50ns simulation of a protein having nearly 700 residues on 60 threads (Gromacs 4.6.3). At one point i got a disk space problem, so i have deleted the md.trr file and created an empty md.trr file. when i tried to restart the simulation from check point file on 100

Re: [gmx-users] hydrogen bond calculation problem

2013-11-12 Thread Justin Lemkul
On 11/11/13 12:36 PM, Sushma Yadav wrote: Dear gromacs users, I have NaCl-water system of finite concentration .I want to calculate water-water hydrogen bond in the first solvation shell around ion..So I used g_hbond -shell option.For single ion in water,its easy to calculate..But say,for wate

[gmx-users] Calculating diffusion coefficient in three dimension

2013-11-12 Thread Venkat Reddy
Dear all, I am simulating a spherical lipid vesicle. I want to calculate the diffusion coefficient for each lipid component in 3D. How to calculate it using g_msd (or any other tool like g_velacc)? Thank you for your concern -- With Best Wishes Venkat Reddy Chirasani -- gmx-users mailing list

[gmx-users] Re: installation error under openSuse 12.2

2013-11-12 Thread kolnkempff
Thank you so much Justin. On the one hand, I feel dumb because I could have sworn that I was using a clean build directory. On the other hand, I obviously lost track of what I was doing because your suggestion worked like a charm! Koln -- View this message in context: http://gromacs.5086.x6.na

[gmx-users] installation error under openSuse 12.2

2013-11-12 Thread Christoph Junghans
> Date: Mon, 11 Nov 2013 08:27:05 -0800 (PST) > From: kolnkempff > Subject: [gmx-users] installation error under openSuse 12.2 > To: gmx-users@gromacs.org > Message-ID: <1384187225465-5012430.p...@n6.nabble.com> > Content-Type: text/plain; charset=us-ascii > > Dear gmx-users, > > I am trying to in

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread bharat gupta
Hi, I tried g_select to dump the structure with the interacting water molecules, but I don't know know how to do that. I searched for some threads in the discussion but wasn't able to find anything related to my need. Can you explain how can I do that ? On Tue, Nov 12, 2013 at 7:39 AM, bharat gu

Re: [gmx-users] Re: g_analyze

2013-11-12 Thread bharat gupta
Sorry, I attached the wrong file . Here's the average file generate from one of the files I sent in my last mail. I used the command g_analyze -f hbond_115-water.xvg -av hbond_115-water-avg.xvg. Here's the file obtained from this command :- https://www.dropbox.com/s/sovzk40cudznfjw/hbond_115-water