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

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

[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 kolnkem...@gmail.com 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

[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:

[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

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

[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

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

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 venkat...@gmail.com 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

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 jalem...@vt.edu wrote: On 11/11/13 12:08 PM, Williams Ernesto Miranda Delgado wrote: Hello If I did the MD simulation using PME and neutralized

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

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 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] 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 jalem...@vt.edu wrote: On

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

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 vvcha...@gmail.comwrote: MSD is 3D

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

[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 You

[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

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,

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

[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

[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

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.

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] 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 jalem...@vt.edu 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*

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 jalem...@vt.edu wrote: On 11/12/13 8:55 AM, Venkat Reddy

[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] 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

[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.gro -p

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,

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

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

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

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.

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 jalem...@vt.edu wrote: On 11/12/13 12:30 PM, Venkat Reddy wrote: Dear Sir, Thanks for the quick reply. So, I have to

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,

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

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 be

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

[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

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

[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 rajatdesi...@gmail.comwrote: 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