[gmx-users] distorted structure

2011-03-30 Thread Olga Ivchenko
Dear gromacs users, I want to simulate small molecule - choline which is not neutral, it has uncompensated positive charge. When I transform mol2 file to pdb structure looks distorted (not correct) in VMD. I added with genion to the system a negative ions, but the pdb structure is still distorted.

Re: [gmx-users] distorted structure

2011-03-30 Thread Mark Abraham
On 30/03/2011 7:50 PM, Olga Ivchenko wrote: Dear gromacs users, I want to simulate small molecule - choline which is not neutral, it has uncompensated positive charge. When I transform mol2 file to pdb structure looks distorted (not correct) in VMD. I added with genion to the system a negativ

Re: [gmx-users] Strange behavior from g_sas

2011-03-30 Thread Mark Abraham
On 30/03/2011 5:58 PM, Terry wrote: Hi, all, I have four trajectories which I want to analyze Surface on. So I run for i in 1 2 3 4; do echo -e " non-Water \n Protein " | g_sas -f $i.xtc -s $i.tpr -o $i.area.xvg -n index.ndx -q $i.pdb ; done But the strange thing is, g_sas analyzed the first

[gmx-users] How to merge 2 top-files

2011-03-30 Thread Ghassen Hassani
Dear Friends, Right now i am working on the docking of a small ligand to a receptor. in order to begin md simulations in gromacs, i need to create some files.. I have already docked my ligand to the receptor and i have those files:    - ligand_receptor.gro    - posre.itp (for only the receptor

Re: [gmx-users] How to merge 2 top-files

2011-03-30 Thread Mark Abraham
On 30/03/2011 8:14 PM, Ghassen Hassani wrote: Dear Friends, Right now i am working on the docking of a small ligand to a receptor. in order to begin md simulations in gromacs, i need to create some files.. I have already docked my ligand to the receptor and i have those files: - ligand_rec

Re: [gmx-users] Strange behavior from g_sas

2011-03-30 Thread Terry
Hi Mark, The output is like this: WARNING: masses and atomic (Van der Waals) radii will be determined based on residue and atom names. These numbers can deviate from the correct mass and radius of the atom type. 2565 out of 4931 atoms were classified as hydrophobic Back Off! I

Re: [gmx-users] Error in GROMACS installation

2011-03-30 Thread Monisha Hajra
Hi Justin, I checked carefully. The make command also gave the same error as pasted below. How should I deal this error? at' follow .libs/libxdrf.o:libxdrf.c:(.text+0x31e3): undefined reference to `_xdr_int' .libs/libxdrf.o:libxdrf.c:(.text+0x340b): undefined reference to `_xdr_int' .libs/libxdrf

Re: [gmx-users] Error in GROMACS installation

2011-03-30 Thread Justin A. Lemkul
Monisha Hajra wrote: Hi Justin, I checked carefully. The make command also gave the same error as pasted below. How should I deal this error? For Windows systems, Gromacs needs to use internal XDR definitions (rather than those that are available system-wide on most Unix/Linux machines).

Re: [gmx-users] Strange behavior from g_sas

2011-03-30 Thread Mark Abraham
On 30/03/2011 9:21 PM, Terry wrote: Hi Mark, The output is like this: WARNING: masses and atomic (Van der Waals) radii will be determined based on residue and atom names. These numbers can deviate from the correct mass and radius of the atom type. 2565 out of 4931 atoms were

[gmx-users] index file

2011-03-30 Thread Elisabeth
Dear gromacs users, I am planning to produce C-C rdf for my system and I have difficulty making the index file with a group of all C in the system, [Carbon] . .If there are thousands of atoms in the system do I have to go through the top file and note down C atom numbers and use those numbers to m

Re: [gmx-users] index file

2011-03-30 Thread Justin A. Lemkul
Elisabeth wrote: Dear gromacs users, I am planning to produce C-C rdf for my system and I have difficulty making the index file with a group of all C in the system, [Carbon] . .If there are thousands of atoms in the system do I have to go through the top file and note down C atom numbers an

[gmx-users] Thermal expansion coeff.

2011-03-30 Thread Elisabeth
Dear all, I asked this question a few days ago but did not receive a solid answer. Has anyone tried to produce the properties below? Can you comment on accuracy of these properties using g_energy? Thermal expansion coeff. (NPT): Enthalpy, Vol, Temp Isothermal compressibility: Vo

[gmx-users] gromacs/CPMD. Is LAM mandatory?

2011-03-30 Thread Elena Formoso
Hi I am trying to do a QM/MM calculation with gromacs/CPMD. When I try to use a parallel version of mdrun and CPMD I get a segmentation fault. I have seen in the examples that LAM is used in runcpmd. It is mandatory for parallel runs? Thanks *** Elena ETH Zürich an

[gmx-users] g_dipole: dipole moment autocorrelation function

2011-03-30 Thread Nilesh Dhumal
Hello, I am trying to calculate the dipole moment autocorrelation function for my system (ionic liquid). I am using gromacs 4.0.7 version. I run the simulation for 4 ns. I run the following command to calculate the dipole moment autocorrelation function g_dipoles -f water.trr -s water.tpr -co

[gmx-users] g_energy -ravg option

2011-03-30 Thread simon sham
Hi, I need to calculate a running average on pressure, and have problem in using the -ravg option in g_energy command: I have tried the followings: 1. g_energy -f file.edr -ravg ravg.xvg -o pressure.xvg 2. g_energy -f file.edr -ravg -o pressure.xvg I could get the fluctuations in pressure.xvg, b

Re: [gmx-users] g_dipole: dipole moment autocorrelation function

2011-03-30 Thread David van der Spoel
On 2011-03-30 18.54, Nilesh Dhumal wrote: Hello, I am trying to calculate the dipole moment autocorrelation function for my system (ionic liquid). I am using gromacs 4.0.7 version. I run the simulation for 4 ns. I run the following command to calculate the dipole moment autocorrelation funct

Re: [gmx-users] g_energy -ravg option

2011-03-30 Thread Justin A. Lemkul
simon sham wrote: Hi, I need to calculate a running average on pressure, and have problem in using the -ravg option in g_energy command: I have tried the followings: 1. g_energy -f file.edr -ravg ravg.xvg -o pressure.xvg 2. g_energy -f file.edr -ravg -o pressure.xvg I could get the fluctuat

Re: [gmx-users] g_dipole: dipole moment autocorrelation function

2011-03-30 Thread Nilesh Dhumal
Thanks. How can I take average. How much long I should run the simulation. Nilesh On Wed, March 30, 2011 12:59 pm, David van der Spoel wrote: > On 2011-03-30 18.54, Nilesh Dhumal wrote: > >> Hello, >> >> >> I am trying to calculate the dipole moment autocorrelation >> function for my system (ion

[gmx-users] Splitted DMPC bilayer

2011-03-30 Thread Dr. Ramón Garduño-Juárez
Dear all, Dear Justin, This time I want to ask the gurus about this problem I encountered in the Equilibration step of my system made of 3 individual (small) protein chains in a solvated DMPC bilayer, no ions present since the protein system is neutral... Following the tutorial I started wit

Re: [gmx-users] g_dipole: dipole moment autocorrelation function

2011-03-30 Thread David van der Spoel
On 2011-03-30 20.16, Nilesh Dhumal wrote: Thanks. How can I take average. summing up and dividing by the number of sims. How much long I should run the simulation. until the average converges. Nilesh On Wed, March 30, 2011 12:59 pm, David van der Spoel wrote: On 2011-03-30 18.54, Nilesh D

Re: [gmx-users] Splitted DMPC bilayer

2011-03-30 Thread Justin A. Lemkul
Dr. Ramón Garduño-Juárez wrote: Dear all, Dear Justin, This time I want to ask the gurus about this problem I encountered in the Equilibration step of my system made of 3 individual (small) protein chains in a solvated DMPC bilayer, no ions present since the protein system is neutral... F

[gmx-users] Heat of vap

2011-03-30 Thread Elisabeth
Dear all, I intend to obtain vaporization heat per volume for a *pure alkane system*. Here is the steps I am taking. Please correct me. 1- Obtain total energy of system (kinetic+potential) and divide by number of molecules to obtain energy per mol of molecules. g_energy -f *.edr -nmol XXX 2- Obta

Re: [gmx-users] Heat of vap

2011-03-30 Thread David van der Spoel
On 2011-03-30 21.22, Elisabeth wrote: Dear all, I intend to obtain vaporization heat per volume for a /pure alkane system/. Here is the steps I am taking. Please correct me. 1- Obtain total energy of system (kinetic+potential) and divide by number of molecules to obtain energy per mol of molec

Re: [gmx-users] Heat of vap

2011-03-30 Thread Justin A. Lemkul
Elisabeth wrote: Dear all, I intend to obtain vaporization heat per volume for a /pure alkane system/. Here is the steps I am taking. Please correct me. 1- Obtain total energy of system (kinetic+potential) and divide by number of molecules to obtain energy per mol of molecules. g_energy -

[gmx-users] Help on umbrella sampling

2011-03-30 Thread raghuvir
Dear Users and Developers I am trying to execute md_pull.mdp {given below} on a small molecule embedded at the core of lipid bilayer. I have tried to modify the md_pull.mdp from the Tutorial by Justin Lemkul on Umbrella sampling. I encounter LINCS warning (1000) after about 100-110 ps of

Re: [gmx-users] Help on umbrella sampling

2011-03-30 Thread Justin A. Lemkul
raghu...@bcpindia.org wrote: Dear Users and Developers I am trying to execute md_pull.mdp {given below} on a small molecule embedded at the core of lipid bilayer. I have tried to modify the md_pull.mdp from the Tutorial by Justin Lemkul on Umbrella sampling. I encounter LINCS warning (10

[gmx-users] the total charge of system is not an integer

2011-03-30 Thread ahmet yıldırım
Dear users, Before energy minimization step , I performed the preprosessing step using grompp . However, there are two note that : *NOTE 1 [file topol.top, line 52]:* System has non-zero total charge: -1.50e+01 *NOTE 2 [file topol.top]:* The largest charge group contains 11 atoms. Sinc

Re: [gmx-users] the total charge of system is not an integer

2011-03-30 Thread Tsjerk Wassenaar
Hi Ahmet, As suggested, it's better to break up your molecule into smaller charge groups. Note that charge groups don't need to have zero charge, nor integer charge. In your case, I'd suggest two COH groups for EDO, which will have zero net charge each, and for TRS I'd take the COH groups as separ

Re: [gmx-users] the total charge of system is not an integer

2011-03-30 Thread Mark Abraham
On 31/03/2011 5:18 PM, ahmet y?ld?r?m wrote: Dear users, Before energy minimization step , I performed the preprosessing step using grompp . However, there are two note that : *_NOTE 1 [file topol.top, line 52]:_* System has non-zero total charge: -1.50e+01 This is an integer. See ht

[gmx-users] Re: gmx-users Digest, Vol 83, Issue 217

2011-03-30 Thread Raghuvir R S Pissurlenkar
, Regards, -- next part -- An HTML attachment was scrubbed... URL: http://lists.gromacs.org/pipermail/gmx-users/attachments/20110330/47fd8ce2/attachment-0001.html -- Message: 2 Date: Wed, 30 Mar 2011 21:31:22 +0200 From: David van der Spoel S