[gmx-users] modification of LJ potential

2013-10-18 Thread Zuzana Benkova
Dear GROMACS users,I went through the mailing list but I am still not sure about my problem. I am going to simulate two chains in a tube using GROMACS. The chains are modeled as strings of connected beads and I would like to describe the interaction with the tube aligned with z axis and passing through [0,0,0] by and additional term to LJ potential in which the r=R-sqrt(x^2+y^2); R is the radius of a tube.I want to simulate this system without pbc. Could you be so kind and give me advice which files of the source code should I modify? Thank you in advance.GreetingsZuzana
-- 
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_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] dispersion corrections in an anisotropic system

2011-08-01 Thread Zuzana Benkova
Dear GROMACS users,my system is composed of a repulsive surface with chains grafted on it and other free chains above this immobilized system. I am hesitating whether to use dispersion corrections for the pressure and energy or not. Without these corrections the repulsion of the grafting surface increases which is understandable. This is however, no a problem because under such conditions I can decrease this repulsion by the manipulation of LJ parameters. However, I think that the dispersion corrections are necessary for representation of the interactions between the bulk free chains wetting the anchored chains. Though in the melt they should be screened, but I am varying the densities of the system and not all cases represent melt situation. I have checked the mailing list but I have not found satisfactory answer. I have not found some way how to treat long range interactions of heterogeneous system.My next question is very trivial, but I can not find in Manual and other sources meaning of AllEnerPres.Thank you for your answer in advance.GreetingsZuzana
-- 
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_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] g_hbond and 4.5.2 version

2011-02-16 Thread Zuzana Benkova
Dear Justin, maybe following information helps you. When I evaluated the number of hydrogen bonds between two different groups, particularly, water and frozen atoms of the siloxane surface I got the same results with both version of GROMACS, 4.0.3 and 4.5.2.GreetingsZuzanaOn 02/15/11, Zuzana Benkova   wrote:Hello Justin, I was using cubic boxes. GreetingsZuzanaOn 02/15/11, Erik Marklund   wrote:Justin A. Lemkul skrev 2011-02-09 23.03:>>>Zuzana Benkova wrote:>>Dear GROMACS users,>>>>I have used g_hbond of version 4.5.2 to analyze number of hydrogen bonds in water. I got the average  number per time frame and number of water oxygen atoms equal to 0.839. When I used g_hbond of version 4.0.7 I got 1.677, which is twice the former value. TIP3P model predicts over 3 hydrogen bonds per one water molecule. I am a bit puzzled. If I multiply the digit from version 4.0.7. by 2 I get the expected number. That is why I supposed that the number of 1.677 means per one water oxygen and  per one water molecule means 2x1.677 since two water molecules participate at one hydrogen bond.>>However, I do not know yet if my interpretation is correct and how to interpret the number obtained by version 4.5.2.>>I would appreciate any help. Thank you in advance.>>>>Try pulling the latest stable development version.  This issue was reported in 4.5.1:>>http://lists.gromacs.org/pipermail/gmx-users/2010-October/054905.html>>but not fixed until after 4.5.3 was released:>>http://lists.gromacs.org/pipermail/gmx-users/2010-December/056406.html>>-Justin>>>Greetings>>Zuzana>>>Are people who are reporting this error using a triclinic boxes or cuboid boxes. That information may help my bugfixing.-- ---Erik Marklund, PhD studentDept. of Cell and Molecular Biology, Uppsala University.Husargatan 3, Box 596,    75124 Uppsala, Swedenphone:    +46 18 471 4537    fax: +46 18 511 755er...@xray.bmc.uu.se    http://folding.bmc.uu.se/-- gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

-- 
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_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] g_hbond and 4.5.2 version

2011-02-15 Thread Zuzana Benkova
Hello Justin, I was using cubic boxes. GreetingsZuzanaOn 02/15/11, Erik Marklund   wrote:Justin A. Lemkul skrev 2011-02-09 23.03:>>>Zuzana Benkova wrote:>>Dear GROMACS users,>>>>I have used g_hbond of version 4.5.2 to analyze number of hydrogen bonds in water. I got the average  number per time frame and number of water oxygen atoms equal to 0.839. When I used g_hbond of version 4.0.7 I got 1.677, which is twice the former value. TIP3P model predicts over 3 hydrogen bonds per one water molecule. I am a bit puzzled. If I multiply the digit from version 4.0.7. by 2 I get the expected number. That is why I supposed that the number of 1.677 means per one water oxygen and  per one water molecule means 2x1.677 since two water molecules participate at one hydrogen bond.>>However, I do not know yet if my interpretation is correct and how to interpret the number obtained by version 4.5.2.>>I would appreciate any help. Thank you in advance.>>>>Try pulling the latest stable development version.  This issue was reported in 4.5.1:>>http://lists.gromacs.org/pipermail/gmx-users/2010-October/054905.html>>but not fixed until after 4.5.3 was released:>>http://lists.gromacs.org/pipermail/gmx-users/2010-December/056406.html>>-Justin>>>Greetings>>Zuzana>>>Are people who are reporting this error using a triclinic boxes or cuboid boxes. That information may help my bugfixing.-- ---Erik Marklund, PhD studentDept. of Cell and Molecular Biology, Uppsala University.Husargatan 3, Box 596,    75124 Uppsala, Swedenphone:    +46 18 471 4537    fax: +46 18 511 755er...@xray.bmc.uu.se    http://folding.bmc.uu.se/-- gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
-- 
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_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_hbond and 4.5.2 version

2011-02-09 Thread Zuzana Benkova
Dear GROMACS users,I have used g_hbond of version 4.5.2 to analyze number of hydrogen bonds in water. I got the average  number per time frame and number of water oxygen atoms equal to 0.839. When I used g_hbond of version 4.0.7 I got 1.677, which is twice the former value. TIP3P model predicts over 3 hydrogen bonds per one water molecule. I am a bit puzzled. If I multiply the digit from version 4.0.7. by 2 I get the expected number. That is why I supposed that the number of 1.677 means per one water oxygen and  per one water molecule means 2x1.677 since two water molecules participate at one hydrogen bond. However, I do not know yet if my interpretation is correct and how to interpret the number obtained by version 4.5.2. I would appreciate any help. Thank you in advance.GreetingsZuzana
-- 
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_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] CHARMM FF and cmap

2010-11-08 Thread Zuzana Benkova
Dear GROMACS users,I would like to run some simulations of systems combined of polyethylene oxide brushes and short homooligopeptides (alanine, 12 residues) using CHARMM force field with GROAMCS 4.0.5 version. I have used parameters from the CHARMM force field files from the newest version of GROMACS. I have prepared rtp files, converted the dihedral potential from type 9 to type 3 (Ryckaert-Bellemans) in ffbonded.itp. I can not prepare topol.top file when using pdb2gms -ff ... -f ... I am geting the following error messageProgram pdb2gmx, VERSION 4.0.3Source code file: resall.c, line: 344Fatal error:in .rtp file in residue cmap at line:   -C  N   CA  C   +NCan I run these simulations with including CMAP.itp in 4.0.5 version. I have tried to omit the cmap.itp in the forcefield.itp, and removed [cmap] section in rtp and thus to neglect the CMAP correction, but I still got the same message. I have checked the mailing list but have not found answer to my questions. I would be very thanful for any suggestions and adviceGreetingsZuzana
-- 
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_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] density tip3p

2010-04-01 Thread Zuzana Benkova
Hi Nisha,I did not try OPLSAA since I need TIP3P model with CHARMM. With this FF hydrogens are also considered as interaction sites. This is not very good combination for GROMACS, but it appeared that the differences are not too big and the time required for simulation with H interacttion side is significantly larger. From my experience rdf is not too sensitive to cut-offs. Try to estimate for instance dielectric constant, dipole moment, diffusion coefficient to be entirely sure. But I think you are getting equlibrated system. Good luckZuzana Dňa 04/01/10, nishap.pa...@utoronto.ca napísal:Hi Zuzana, You are right. I did get my density to ~0.986g/cm3, and I actually came across this other paper by Price and Brooks, and initially they get a density of Tip3p using pme at around 0.979g/cm3, but they tweak some of the parameters in oplsaa to get the correct density. I am not sure if you have tried that and if that helps at all, but also they are using a different simulation package.I am getting good rdf values for the water even with slightly low density, so I guess it is okay.ThanksQuoting Zuzana Benkova :>Hi Nisha,>>I have performed a benchmark simulations of TIP3P water and I still>obtained density 0.984g.cm-3. It looks that this model gives such a>density. If you look at JCP 79, 926, 1983 you will see that the>authors got 0.982g.cm-3 with MC. (Table 3) Try to evaluate alsoo>other properties, each model of water has some deficiencies. I tried>large number of cut-offs combinations LJ and QQ tyoes of>interactions, but the density did not gange by more than 0.001g.cm-3.>>>Greetings>>Zuzana>>Dňa 03/31/10, nishap.pa...@utoronto.ca napísal:Hello,>>    So I looked at the manual and I changed my mdp parameters as>follows to equilibrate at a constant pressure.>> RUN CONTROL PARAMETERS =>integrator   = sd>tinit    = 0>dt   = 0.002>nsteps   = 50>nstcomm  = 1>nstxout  = 0>nstvout  = 0>nstfout  = 0>nstlog   = 500>nstenergy    = 500>nstxtcout    = 0>xtc-precision    = 1000>nstlist  = 5>ns_type  = grid>pbc  = xyz>rlist    = 1.3>coulombtype  = pme>rcoulomb = 1.3>epsilon-r    = 1>vdw-type = switch>rvdw-switch  = 1.0>rvdw = 1.2>DispCorr = EnerPres>fourierspacing   = 0.1>fourier_nx   = 0>fourier_ny   = 0>fourier_nz   = 0>pme_order    = 6>ewald_rtol   = 1e-06>epsilon_surface  = 0>optimize_fft = no>tc_grps  = system>tau_t    = 0.1>ref_t    = 298>Pcoupl   = berendsen>tau_p    = 1.0>compressibility  = 4.5e-05>ref_p    = 1.0>constraints  = all-bonds>constraint-algorithm = Lincs>unconstrained-start  = no>lincs-order  = 4>lincs-warnangle  = 30>lincs-iter   = 1>gen_vel  = yes>gen_temp = 298.0>gen_seed = 173529>>I also tried using mdp parameters as mentioned in some of the gromacs>tutorials. I ran this simulation for 1ns and I am still getting>density around 985g/l, which is low for tip3p water (27nm^3, 900>water molecules). I checked some of the publications and the values>they obtained were ~998g/l. I checked for vacuum in my system, but I>didn't see anything out of the ordinary. Any insights?>>Thanks>>Nisha P>>Quoting nishap.pa...@utoronto.ca:>>>Thanks. I will look up the manual again.>>>>Quoting Sander Pronk :>>>>>Hi Nisha,>>>>>>Looking at your .mdp, there are some issues that might lead to the>>>behavior that you describe:>>>First: you should try to look up the published densities for tip3p>>>water at 300K - they might actually be close to what you get.>>>Second: your neighbor list cut-off (rlist) might be too small to>>>fully contain the charge groups (check the manual, section 7.3.11).>>>Third: You haven't enabled long-range mean field correction for the>>>pressure or energy. Expect the pressure to be strongly dependent on>>>cut-off (see the same section).>>>>>>Sander>>>>>>>>>On Mar 30, 2010, at 22:47 , nishap.pa...@utoronto.ca wrote:>>>>>>>I have a box (3x3x3nm) of Tip3p water molecules ~900 and the>>>>density when I create the box u

Re: [gmx-users] density tip3p

2010-04-01 Thread Zuzana Benkova
Hi Nisha,I have performed a benchmark simulations of TIP3P water and I still obtained density 0.984g.cm-3. It looks that this model gives such a density. If you look at JCP 79, 926, 1983 you will see that the authors got 0.982g.cm-3 with MC. (Table 3) Try to evaluate alsoo other properties, each model of water has some deficiencies. I tried large number of cut-offs combinations LJ and QQ tyoes of interactions, but the density did not gange by more than 0.001g.cm-3. GreetingsZuzanaDňa 03/31/10, nishap.pa...@utoronto.ca napísal:Hello,    So I looked at the manual and I changed my mdp parameters as follows to equilibrate at a constant pressure. RUN CONTROL PARAMETERS =integrator   = sdtinit    = 0dt   = 0.002nsteps   = 50nstcomm  = 1nstxout  = 0nstvout  = 0nstfout  = 0nstlog   = 500nstenergy    = 500nstxtcout    = 0xtc-precision    = 1000nstlist  = 5ns_type  = gridpbc  = xyzrlist    = 1.3coulombtype  = pmercoulomb = 1.3epsilon-r    = 1vdw-type = switchrvdw-switch  = 1.0rvdw = 1.2DispCorr = EnerPresfourierspacing   = 0.1fourier_nx   = 0fourier_ny   = 0fourier_nz   = 0pme_order    = 6ewald_rtol   = 1e-06epsilon_surface  = 0optimize_fft = notc_grps  = systemtau_t    = 0.1ref_t    = 298Pcoupl   = berendsentau_p    = 1.0compressibility  = 4.5e-05ref_p    = 1.0constraints  = all-bondsconstraint-algorithm = Lincsunconstrained-start  = nolincs-order  = 4lincs-warnangle  = 30lincs-iter   = 1gen_vel  = yesgen_temp = 298.0gen_seed = 173529I also tried using mdp parameters as mentioned in some of the gromacs tutorials. I ran this simulation for 1ns and I am still getting density around 985g/l, which is low for tip3p water (27nm^3, 900 water molecules). I checked some of the publications and the values they obtained were ~998g/l. I checked for vacuum in my system, but I didn't see anything out of the ordinary. Any insights?ThanksNisha PQuoting nishap.pa...@utoronto.ca:>Thanks. I will look up the manual again.>>Quoting Sander Pronk :>>>Hi Nisha,Looking at your .mdp, there are some issues that might lead to the>>behavior that you describe:>>First: you should try to look up the published densities for tip3p>>water at 300K - they might actually be close to what you get.>>Second: your neighbor list cut-off (rlist) might be too small to>>fully contain the charge groups (check the manual, section 7.3.11).>>Third: You haven't enabled long-range mean field correction for the>>pressure or energy. Expect the pressure to be strongly dependent on>>cut-off (see the same section).Sander>>On Mar 30, 2010, at 22:47 , nishap.pa...@utoronto.ca wrote:>I have a box (3x3x3nm) of Tip3p water molecules ~900 and the>>>density when I create the box using genbox is 997.177g/l. I did>>>energy minimization run and the potential energy did converge>>>smoothly, so I did NPT equilibration run of 100ps and I got the>>>density value of 975g/l. Why does the density decrease after the>>>run? these are the parameters that I used :>>RUN CONTROL PARAMETERS>>>integrator   = md>>>tinit    = 0>>>dt   = 0.002>>>nsteps   = 5>>>nstcomm  = 0>>>nstxout  = 0>>>nstvout  = 0>>>nstfout  = 0>>>nstlog   = 100>>>nstenergy    = 100>>>nstxtcout    = 0>>>xtc_precision    = 1000>>>nstlist  = 5>>>ns-type  = Grid>>>pbc  = xyz>>>rlist    = 1.1>>>coulombtype  = pme>>>rcoulomb = 1.1>>>epsilon-r    = 1>>>vdw-type = switch>>>rvdw-switch  = 0.8>>>rvdw = 1.0>>>Tcoupl   = V-rescale>>>tc-grps  = System>>>tau_t    = 0.1>>>ref_t    = 300>>>Pcoupl   = berendsen>>>tau_p    = 0.5>>>compressibility  = 4e-05>>>ref_p    = 1.0>>>constraints  = all-bonds>>>constraint-algorithm = Lincs>>>unconstrained-start  = no>>>lincs-order  = 4>>>lincs-iter   = 1>>>lincs-warnangle  = 30>>Am I missing something?>>>Thanks>>Nisha Patel>-->>>gmx-users mailing list    gmx-users@gromacs.org>>>http://lists.gromacs.org/mailman/listinfo/gmx-users>>>Please s

Re: [gmx-users] T-coupling groups

2010-03-19 Thread Zuzana Benkova
Dear Srinivs,thank you for response. You are right I have already found the source of problem. I have used PEO to name the chain group and forgot that the same group holded for residues. GreetingsZuzanaDňa 03/19/10, jampani srinivas   napísal:Dear Zuzana,

Recently i have done similar selection, I have got this error only when i miss some atoms or when i have some atoms in both groups. I would suggest you to check your groups carefully or you can make indexes for these selections in VMD and check whether you are missing some atoms from both groups.


Srinivs.
On Thu, Mar 18, 2010 at 4:55 PM, Zuzana Benkova <zuzana.benk...@savba.sk> wrote:

Dear GROMACS users,I have checked the mails in archive which deal with the problem I am facing but did not find the answer what I do wrong. 
I want to compile a PEO chain centered in box of 2900 water molecules. In my md.mdp the relevant section is as followsTcoupl  =  nose-hoovertc-grps =  PEO  SOLtau_t   =  0.1  0.1
ref_t   =  50   50gen_vel =  yesgen_temp    =  50   50I have generated an index file where I have defined the PEO group. The [ SOL ] group contains 129-8828 atoms and the [ PEO ] contains 1-128atoms. This matches with the gro file. I tried to prepare a tpr file using
grompp -f md.mdp -c PEO18_2900TIP3P_opt.gro -n PEO18_2900TIP3P.ndx  -p topol.top -o PEO18_2900TIP3P_50Kand got messageProgram grompp, VERSION 4.0.3Source code file: readir.c, line: 1050Fatal error:
16 atoms are not part of any of the T-Coupling groupsI tried the compilations with 4.0.3 and 4.0.5 versions but the message was the same.  I would be really thankful for an advice which can move me further. Thank you in advance.
Zuzana--gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!Please don't post (un)subscribe requests to the list. Use thewww interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php-- *
J. Srinivasa Rao Ph.D.Post-doctoral Research AssociateC/o Prof. Luis R Cruz CruzComputational Biophysics GroupDepartment of PhysicsDrexel University3141 Chestnut StPhiladelphia, PA 19104, USA.
Ph:  Off:     215-895-1989       Mob:  704-706-4191Web:http://jsrao.web.officelive.com/default.aspx**


-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] T-coupling groups

2010-03-18 Thread Zuzana Benkova
Dear GROMACS users,
I have checked the mails in archive which deal with the problem I am 
facing but did not find the answer what I do wrong. 

I want to compile a PEO chain centered in box of 2900 water molecules. 

In my md.mdp the relevant section is as follows
Tcoupl  =  nose-hoover
tc-grps =  PEO  SOL
tau_t   =  0.1  0.1
ref_t   =  50   50
gen_vel =  yes
gen_temp    =  50   50

I have generated an index file where I have defined the PEO group. The [ SOL ] group contains 129-8828 atoms and the [ PEO ] contains 1-128atoms. This matches with the gro file. I 
tried to prepare a tpr file using

grompp -f md.mdp -c PEO18_2900TIP3P_opt.gro -n PEO18_2900TIP3P.ndx  -p 
topol.top -o PEO18_2900TIP3P_50K

and got message

Program grompp, VERSION 4.0.3
Source code file: readir.c, line: 1050

Fatal error:
16 atoms are not part of any of the T-Coupling groups

I tried the compilations with 4.0.3 and 4.0.5 versions but the message 
was the same.  I would be really thankful for an advice which can move me 
further. Thank you in advance.

Zuzana




-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] topol.top

2010-03-04 Thread Zuzana Benkova
Dear Mark,I have tried the same processing  with version 4.0.5. After that I got the following error messageProgram grompp_s, VERSION 4.0.5Source code file: confio.c, line: 748Fatal error:The spacing of the decimal points in file conf.gro is not consistent for x, y and zMy conf.gro file was constructing according to the required format. I aslo compared it with other gro files, that worked. After replacing replacing of zero .0 by 0.0 the procesing was successful. Thank you for hint.ZuzanaDňa 03/04/10, Mark Abraham   napísal:On 5/03/2010 4:08 AM, Zuzana Benkova wrote:>Dear GROMACS users,>I wanted to do a test calculation on a linear unrealistic molecule PPPBC>without any contributions to energy exept of LJ nonbenoded interactions.>I have prepared a topol.top file which I am attaching. After processing>this file with grompp I obtained a message>>WARNING 1 [file topol.top, line 63]:>   Too few parameters on line (source file confio.c, line 776)>>Line 63 is the last line of topol.top. I followed the manual and I am>not able to fix the problem. Please could you have a look at this file>and possibly to find out what is missing. Thank you in advanceWhat is your GROMACS version?Mark-- gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease search the archive at http://www.gromacs.org/search before posting!Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] topol.top

2010-03-04 Thread Zuzana Benkova
Dear Mark,I have used version 4.0.3.ZuzanaDňa 03/04/10, Mark Abraham   napísal:On 5/03/2010 4:08 AM, Zuzana Benkova wrote:>Dear GROMACS users,>I wanted to do a test calculation on a linear unrealistic molecule PPPBC>without any contributions to energy exept of LJ nonbenoded interactions.>I have prepared a topol.top file which I am attaching. After processing>this file with grompp I obtained a message>>WARNING 1 [file topol.top, line 63]:>   Too few parameters on line (source file confio.c, line 776)>>Line 63 is the last line of topol.top. I followed the manual and I am>not able to fix the problem. Please could you have a look at this file>and possibly to find out what is missing. Thank you in advanceWhat is your GROMACS version?Mark-- gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease search the archive at http://www.gromacs.org/search before posting!Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] topol.top

2010-03-04 Thread Zuzana Benkova
Dear GROMACS users,
I wanted to do a test calculation on a linear unrealistic molecule
PPPBC without any contributions to energy exept of LJ nonbenoded
interactions. I have prepared a topol.top file which I am attaching.
After processing this file with grompp I obtained a message

WARNING 1 [file topol.top, line 63]:
  Too few parameters on line (source file confio.c, line 776)

Line 63 is the last line of topol.top. I followed the manual and I am
not able to fix the problem. Please could you have a look at this file
and possibly to find out what is missing. Thank you in advance
Zuzana



topol.top
Description: Binary data
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] topol.top

2010-03-04 Thread Zuzana Benkova
Dear GROMACS users,I wanted to do a test calculation on a linear unrealistic molecule PPPBC without any contributions to energy exept of LJ nonbenoded interactions. I have prepared a topol.top file which I am sending. After processing this file with grompp I obtained a messageWARNING 1 [file topol.top, line 63]:  Too few parameters on line (source file confio.c, line 776)Line 63 is the last line of topol.top. I followed the manual and I am not able to fix the problem. Please could you have a look at this file and possibly to find out what is missing. Thank you in advanceZuzana[ defaults ]; nbfunc    comb-rule   gen-pairs   fudgeLJ fudgeQQ  1 2   yes 1.0 1.0[ atomtypes ];at. type at.num mass  charge    ptype sigma   epsilonP   1 1  0. A 1   1B   1 1  0. A 2   2C   1 1  0. A 3   1.5[ bondtypes ]; i    j func    b0  kbP  P    1 1  0.P  B    1 1  0.B  C    1 1  0.[ angletypes ];  i    j    k  func   th0   cthP   P    P   1 0 0.P   B    C   1 0 0.P   P    B   1 0 0.[ dihedraltypes ]P   P   P   B 1 0    0    1P   P   B   C 1 0    0    1[ moleculetype ]; Name    nrexclProtein 3[ atoms ];   nr   type  resnr residue  atom   cgnr charge   mass1   P   1  R    P 1 0   12   P   1  R    P 1 0   13   P   1  R    P 1 0   14   B   1  R    B 1 0   15   C   1  R    C 1 0   1[ bonds ];  ai    aj funct    c0    c1    c2    c31    2  12    3  13    4  14    5  1[ angles ];  ai    aj    ak funct    c0    c1    c2    c31  2    3  12  3    4  13  4    5  1[ dihedrals ];  ai    aj    ak    al funct1    2   3    4   12    3   4    5   1[ system ]; NameProtein[ molecules ]; Compound    #molsProtein    1
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] end-to-end distance

2010-02-18 Thread Zuzana Benkova
Dear Justin,thank you for response and I am sorry not to define my problem properly. I have used the first alternative, which you discribed and thus obtained the weird resluts. Before I used this alternative I did what you mentioned as the second alternative if I understood it correctly. I have specified the first and the last atoms of each chain as separate groups. But this makes 64 pairs of groups for 64 chairs and I need to apply 64 times g_dist, please, could you correct me if I am wrong? If I am right than it is better to write a script. I am appending the example of my group definitions [ PEO1_&_CH3 ]1   87  173  259  345  431  517  603  689  775  861  947 1033 1119 12051291 1377 1463 1549 1635 1721 1807 1893 1979 2065 2151 2237 2323 2409 24952581 2667 2753 2839 2925 3011 3097 3183 3269 3355 3441 3527 3613 3699 37853871 3957 4043 4129 4215 4301 4387 4473 4559 4645 4731 4817 4903 4989 50755161 5247 5333 5419[ PEON_&_CH3 ]83  169  255  341  427  513  599  685  771  857  943 1029 1115 1201 12871373 1459 1545 1631 1717 1803 1889 1975 2061 2147 2233 2319 2405 2491 25772663 2749 2835 2921 3007 3093 3179 3265 3351 3437 3523 3609 3695 3781 38673953 4039 4125 4211 4297 4383 4469 4555 4641 4727 4813 4899 4985 5071 51575243 5329 5415 5501which provided weird results[ R1 ]1[ R11 ]83[ R2 ]
87

[ R22]
169till [ R64 ] [ R6464 ] pairs which step-by step gave realistic results.
Thank you in advanceZuzana.Dňa 02/18/10, "Justin A. Lemkul"   napísal:Zuzana Benkova wrote:>Dear GROMACS users,>>I have a system of 64 chains and want to calculate averaged values of end-to-end distance of all chains as well as the averaged radius of gyration from 1ns trajectory. In the case of end-to-end distance I have prepared index file containing one group of the first atoms in ascending order and the second one containing the terminal atoms in ascending order. I have used g_dist but it looks like it combined all pairs between both groups. I have looked at the mailing list but did not find the answer.>g_dist should only be measuring the COM distance between the two selected groups, so I don't know what you mean by "combined all pairs" - perhaps you can provide an example (index groups, relevant output, what you were expecting, and why you think it's wrong)?  If your groups contain all starting residues in one, then all ending residues in another, then yes, you're going to get weird results.  If you have one chain that is residues 1-100, you should have a group for r_1 and a group for r_100; select those for g_dist analysis.  If another chain is residue 101-200, r_101 and r_200 are what you want, etc.>Is it possible to do such analysis using directly GROMACS tools or I need to prepare some scripts with loops running over all chains. Thank you for response in advance.>It certainly should be possible to iteratively use g_dist to measure distances between whatever groups you select.  Fairly straightforward, actually, so if you're not getting what you expect, provide the information I asked for above and hopefully someone should be able to figure out what's going wrong.-Justin>Zuzana Benkova>-- Justin A. LemkulPh.D. CandidateICTAS Doctoral ScholarMILES-IGERT TraineeDepartment of BiochemistryVirginia TechBlacksburg, VAjalemkul[at]vt.edu | (540) 231-9080http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin-- gmx-users mailing list    gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease search the archive at http://www.gromacs.org/search before posting!Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] end-to-end distance

2010-02-18 Thread Zuzana Benkova
Dear GROMACS users, I have a system of 64 chains and want to calculate averaged values of end-to-end distance of all chains as well as the averaged radius of gyration from 1ns trajectory. In the case of end-to-end distance I have prepared index file containing one group of the first atoms in ascending order and the second one containing the terminal atoms in ascending order. I have used g_dist but it looks like it combined all pairs between both groups. I have looked at the mailing list but did not find the answer. Is it possible to do such analysis using directly GROMACS tools or I need to prepare some scripts with loops running over all chains. Thank you for response in advance.Zuzana Benkova
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] user defined potential

2010-02-17 Thread Zuzana Benkova
Dear GROMACS users,if I want to use table defined potential for nonbonded interactions I have to do so for all pairs between different grops or is there a possibility to specify just one kind lets say BH potential by tables and to use wdvtype=cut-off for LJ interactions. It seems to me that this is not possible. Could you be so kind and confirm my opinion that if I want to mix BH and LJ potential I have to define all interactions by tables. Zuzana Benkova
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] size dependence of fluctuations

2009-10-02 Thread Zuzana Benkova
Dear GROMACS users,

most likely I am asking a trivial question but I need to be sure about it
for calculations of heat capacity based on the fluctuations of energy in
canonical ensemble. I am doing simulation of water using TIP4P model. I
compare two systems consisting of 216 and 8*216 = 1728 molecules. The
expected fluctuations for the larger system should be reduced 8^(1/2)=2.83
times with respect to the fluctuations of the smaller system. Actually, this
applies only to the fluctuations in intensive properties such as the
temperature, pressure and density while the fluctuations in extensive
properties such as energies exceed about 2.83 times the corresponding
fluctuations for smaller system. Could you explain me the inverse ratio of
fluctuations? 

I have also one question regarding the heat capacity calculation. If I want
to apply the formula for canonical ensemble (-^2)/(kbT^2) than the
capacity should be size dependent. I am a bit puzzled about it and any
explanation would be very helpful for me. Thanks in advance

 

Zuzana.

 

___
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] No default LJ-14 types

2009-05-12 Thread Zuzana Benkova
Dear GROMACS users,

I want to use GROMOS force field for the polyethylene simulations to compare
the performance of this force field with those generated for alkanes. I have
prepared my own files based on the ffG53a5 files. In itp I used

[ defaults ]

; nbfunccomb-rule   gen-pairs   fudgeLJ fudgeQQ

  1   1   no  1.0 1.0

 

and in rtp  

[ bondedtypes ]

; bonds  angles  dihedrals  impropers

2   2  1  2

 

In my topology file I required nrexcl = 3 in moleculetype section. I
included [ pairs ] where all 1-4 pairs were given. 

When I did the same for polyalanine system in which case I have used the
ffG53a5 directly from top library, the simulation run and I got  the 1-4
contribution to the potential energy in log file. However in the case of
polyethylene after preprocessing using grompp I got error 41 messages for
each of 41 pairs of the following kind

 

ERROR 35 [file C44.itp, line 134]:

  No default LJ-14 types

 

 From manual I expect that if gen-pairs is no, 1-4 interaction are taken
from pairs section and if yes the 1-4 interaction not presented in pairs
section are generated with corresponding fractions as defined in fudgeLJ and
fudgeQQ. 

 

Previously I removed all the pairs section, simulations run but without 1-4
LJ contribution. I have the same problem with OPLS force field, with which
despite the Ryckaert-Bellemans dihedral potential the 1-4 interaction should
be considered. 

 

I would be very thankful for help. Thank you in advance

 

 

Zuzana

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

RE: [gmx-users] UA simulation of PEO melt

2009-05-12 Thread Zuzana Benkova
Hi Patrick,
Thank you very much for your response. I have just downloaded the paper of
which you sent the reference. I have not found other FF parameters except
for the torsional one about the C-C bond. I have also some papers which
provide good results on PEO simulations (JCP, 124, 234901 (2006); JCP, 102,
9725 (1995).) but I am limited with the order of Ryckaert-Bellemans
potential. Unfortunately I do not have any idea how to perform the code
modification to extent the torsional expansion. This was the second part of
my first mail. I hope somebody will provide me with details necessary to
modify the code appropriately and with all the changes associated with that.


Greetings
Zuzana

-Original Message-
From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org]
On Behalf Of patrick fuchs
Sent: Tuesday, May 12, 2009 3:54 PM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] UA simulation of PEO melt

Hi Zuzana,
I'm actually working on the paramaterization of polyethers with the 
GROMOS developers. Currently in G53a?, there's no such torsion O-C-C-O 
in GROMOS and I don't advice to use any combination of existing 
monoethers. I could figure out that they give the wrong free energy of 
hydration. I hope to have the parameters ready before the end of the 
year. Meanwhile, you can use the recently optimized parameters for 
CHARMM (http://dx.doi.org/10.1529/biophysj.108.133025) though they are 
all-atom.
Cheers,

Patrick

Zuzana Benkova a écrit :
> Dear GROMACS users
> I  want to run simulations of polyethyleneoxide melt. In literature I 
> have found united atoms force fields that performed well for this kind 
> of simulation. For dihedral potential I need to extend the 
> Ryckaert-Bellemans potential up to the order of six. I have looked at 
> the mailing list and found some information which give advice to change 
> the code for this kind of potential. But I could not figure out where to 
> do this changes and how to do them. I would need some detail to do it. I 
> do not have experience with something like this.
> 
> My second question also concerns PEO. I tried to use GROMOS force field 
> as well for simulations. I used the parameters from ffG53a5 files. In 
> ffG53a5bon.itp I am not sure which code of dihedrals sould be used for 
> C-C-O-C and O-C-C-O and which one for the bond angles. In the case of 
> bond angle I found type ga_8 as suitable, but it seems me from the 
> description that it sould be the same for all three bonds in PEO. This 
> seems to be not OK, since the bending of C-O-C angle is characterized by 
> larger bending constant.
> 
> I would appreciate any help. Thank you in advance
> 
> Zuzana.
> 
>  
> 
> 
> 
> 
> ___
> gmx-users mailing listgmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-requ...@gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php

-- 
___
 new E-mail address: patrick.fu...@univ-paris-diderot.fr 
Patrick FUCHS
Dynamique des Structures et Interactions des Macromolécules Biologiques
INTS, INSERM UMR-S665, Université Paris Diderot,
6 rue Alexandre Cabanel, 75015 Paris
Tel : +33 (0)1-44-49-30-57 - Fax : +33 (0)1-47-34-74-31
Web Site: http://www.dsimb.inserm.fr/~fuchs

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] UA simulation of PEO melt

2009-05-12 Thread Zuzana Benkova
Dear GROMACS users
I  want to run simulations of polyethyleneoxide melt. In literature I have
found united atoms force fields that performed well for this kind of
simulation. For dihedral potential I need to extend the Ryckaert-Bellemans
potential up to the order of six. I have looked at the mailing list and
found some information which give advice to change the code for this kind of
potential. But I could not figure out where to do this changes and how to do
them. I would need some detail to do it. I do not have experience with
something like this. 

My second question also concerns PEO. I tried to use GROMOS force field as
well for simulations. I used the parameters from ffG53a5 files. In
ffG53a5bon.itp I am not sure which code of dihedrals sould be used for
C-C-O-C and O-C-C-O and which one for the bond angles. In the case of bond
angle I found type ga_8 as suitable, but it seems me from the description
that it sould be the same for all three bonds in PEO. This seems to be not
OK, since the bending of C-O-C angle is characterized by larger bending
constant. 

I would appreciate any help. Thank you in advance

Zuzana.

 

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

RE: [gmx-users] Polyethylene (PE) simulations in Gromacs - please!

2009-03-03 Thread Zuzana Benkova
Dear Chansoo Kim,

It has passed relatively long time since I attempted to use pdb2gmx to build
my topology for polyethylene. I prepared the relevant files according to
those used for amino acids. I used + and - sign to discern the direction of
connectivity. I have encountered dome problems mentioned in the mail but
they have not been answered satisfactorily. Later on I wanted to run just
united atom simulations and I had next problem with pdb2gmx that always
include hydrogen atoms to my topology though I used option which should have
disabled this. I decided to build my topology file without pdb2gmx and I
applied force field that is not present in GMX directories. 

 

  _  

From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org]
On Behalf Of Chansoo Kim
Sent: Tuesday, March 03, 2009 2:50 PM
To: jalem...@vt.edu; Discussion list for GROMACS users
Subject: Re: [gmx-users] Polyethylene (PE) simulations in Gromacs - please!

 

Dear Justin:

 

 

Thank you for your reply.

I was too urgent, so I did not clarify my questions.

Sorry about that.

 

I totally understand how topology files (tdb, hdb, and etc..) work, but I do
not clearly understand chemical reaction, which I have to define in those
files...

 

 

[Situation]

I want to do simulate polyethylene (PE) using Gromacs and OPLS-AA force
field...

 

[Topology files and etc]

  (1) rtp file

[ Eth ]
 [ atoms ]
   C1opls_136-0.1201
   H11   opls_140 0.0601
   H12   opls_140 0.0601
   C2opls_136-0.1202
   H21   opls_140 0.0602
   H22   opls_140 0.0602
 [ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
   C2+C1
 
  (2) hdb file
Eth 2
2   6   H1   C1  C2 +C1
2   6   H2   C2  C1 -C2

 

 
  (3) -c.tdb file
 
[ Eth ]
[ replace ]
C1  opls_13512.011  -0.18
[ add ]
34  H1  C1  C2
 opls_140   1.008   0.06
[ delete ]
H21
H22

 

 
  (4) -n.tdb file
 
[ Eth ]
[ replace ]
C2  opls_13512.011   -0.18
[ add ]
34  H2  C2  C1
 opls_140   1.008   0.06
[ delete ]
H21
H22

 

   (5) input pdb file

 

ATOM  1  C1  Eth 1   1.000   1.540   0.000

ATOM  2  C2  Eth 1   2.456   2.041   0.000
ATOM  3  C1  Eth 2   2.456   3.581   0.000
ATOM  4  C2  Eth 2   3.912   4.083   0.000
ATOM  5  C1  Eth 3   3.912   5.623   0.000
ATOM  6  C2  Eth 3   5.368   6.124   0.000
END

 

[Error]

When I did run the pdb2gmx, I got the following error.

"Atom -C not found in residue PEth1 while adding hydrogens"

 

[Question #1 and #2]

Therefore, what I should define more...?

What does the "+" and "-" mean here?

 

 

[If...]

If I change the hdb file as follows,

 

 
Eth 2
2   6   H1   C1  C2 C1
2   6   H2   C2  C1 C2

 

After doing pdb2gmx, I got only 12 H-atoms (which means my terminal did NOT
work...).

And those positions are all (0, 0, 0).

 

Following is the results,

 

LE Giving Russians Opium May Alter Current Situation

MODEL1

ATOM  1  C1  PEt 1   1.000   0.000   0.000  1.00  0.00

ATOM  2  H11 PEt 1   0.000   0.000   0.000  1.00  0.00

ATOM  3  H12 PEt 1   0.000   0.000   0.000  1.00  0.00

ATOM  4  C2  PEt 1   2.450   6.000   1.000  1.00  0.00

ATOM  5  H21 PEt 1   0.000   0.000   0.000  1.00  0.00

ATOM  6  H22 PEt 1   0.000   0.000   0.000  1.00  0.00

ATOM  7  C1  PEt 2   2.450   6.000   1.000  1.00  0.00

ATOM  8  H11 PEt 2   0.000   0.000   0.000  1.00  0.00

ATOM  9  H12 PEt 2   0.000   0.000   0.000  1.00  0.00

ATOM 10  C2  PEt 2   3.910   2.000   3.000  1.00  0.00

ATOM 11  H21 PEt 2   0.000   0.000   0.000  1.00  0.00

ATOM 12  H22 PEt 2   0.000   0.000   0.000  1.00  0.00

ATOM 13  C1  PEt 3   3.910   2.000   3.000  1.00  0.00

ATOM 14  H11 PEt 3   0.000   0.000   0.000  1.00  0.00

ATOM 15  H12 PEt 3   0.000   0.000   0.000  1.00  0.00

ATOM 16  C2  PEt 3   5.360   8.000   4.000  1.00  0.00

ATOM 17  H21 PEt 3   0.000   0.000   0.000  1.00  0.00

ATOM 18  H22 PEt 3   0.000   0.000   0.000  1.00  0.00

TER

ENDMDL

 

[Question #3 from If... section ]

I guess that I have to change hdb and tdb files.

Would you please give some guide on that...?

 

Thank you for your care!

 

 

Sincerely yours,

 

C Kim

 

On Tue, Mar 3, 2009 at 8:35 AM, Justin A. Lemkul  wrote:



Chansoo Kim wrote:

Dear Dr. Benkova and others:


I am C. Kim and trying to simulation polymer system.

Since I guessed that polyethylene (PE) is a simple system, it could be not
that hard to simulate it in Gromacs.
Trying to do, I have read Dr. Benkova's articles in gmx-users mailing list
to tackle my problems.
Sorry to say, I

RE: [gmx-users] running a simulation without production

2008-11-06 Thread Zuzana Benkova


-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Zuzana Benkova
Sent: Tuesday, November 04, 2008 3:11 PM
To: 'Discussion list for GROMACS users'
Subject: RE: [gmx-users] running a simulation without production

Dear Mark, 
Thank you very much for valuable advice. 
Previously, I run some simulations and I obtained what I had expected. For
instance, instead of the plate I have just used one single chain which was
treated as frozen and the simulations were running. I obtained all output
files containing appropriate information. 
I have GMX version 3.3.3. I have just tried gmxtest.pl test first in sthe
ingle precision. The output I got was

All 16 simple tests PASSED
FAILED. Check files in dec+water
FAILED. Check files in fe_test
2 out of 14 complex tests FAILED
All 63 kernel tests PASSED
All 45 pdb2gmx tests PASSED

I am not very familiar with programming language C++. As I wanted to test
also the double precision I did some changes in the gmxtest.pl file. I am
not sure if I did it correctly. I have installed GMX with _mpi suffix, so I
have changed the sub setup_vars section 

sub setup_vars {
my $mdprefix = "";
if ( $parallel > 0 ) {
$mdprefix = "mpirun -c $parallel ";
}

if ( $double > 0 ) {
$grompp = "grompp_mpi";
$mdrun  = $mdprefix . "mdrun_mpi";
$ref= "reference_d";
}

and after applying command gmxtest.pl -double all I have obtained 

sh: line 1: 30240 Segmentation fault  mdrun_mpi >mdrun.out 2>&1
FAILED. Check files in rb1
1 out of 16 simple tests FAILED
All 14 complex tests PASSED
All 63 kernel tests PASSED
All 45 pdb2gmx tests PASSED

In the case of the single precision test I can not figure out the error
since I do not have any new files produced after the test in corresponding
directories. In the case of double precision test I have also log file in
rb1 directory which has been terminated exactly in the same way as my
attempts to run simulation of plate-chain system. 

Please, could you give me some explanation and hints how to proceed further?
It would be very helpful for me. I have checked the mailing list but did not
find something that could help me. Thank you in advance.

Zuzana


-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Mark Abraham
Sent: Monday, November 03, 2008 9:55 PM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] running a simulation without production

Zuzana Benkova wrote:
> I have tried also running this job without user specified potential, 
> with relaxed plate, with pbc=no and some other attempts, but I always 
> failed in the same way. I would be very thankful for suggestions how to 
> cope with this problem.

Do other simulations run? Try one of the tutorials or gmxtest.pl.

Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


RE: [gmx-users] running a simulation without production

2008-11-05 Thread Zuzana Benkova
 

 

  _  

From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Zuzana Benkova
Sent: Monday, November 03, 2008 5:49 PM
To: gmx-users@gromacs.org
Subject: [gmx-users] running a simulation without production

 

Dear GROMACS users,

 

I want to simulate a system consisting of a plate modeled by atoms (5000
atoms) and a polyethylene chain (80 carbons). The plate is frozen and I am
using tabulated potentials. In this system, the considered interactions are
chain-chain and chain-plate. There are no electrostatic interactions. The
system is placed in a box, with dimensions corresponding to the dimensions
of plate and substantially larger height. My md.mdp file is as follows

 

cpp =  cpp

define  =

 

integrator  =  md

dt  =  0.002

nsteps  =  7500

nstcomm =  1

comm_mode   =  None

nstxout =  10

nstvout =  10

nstfout =  10

nstxtcout   =  10

nstlog  =  10

nstenergy   =  10

energygrps  =  Chain Plate

freezegrps  =  Plate

freezedim   =  Y Y Y

energygrp_excl  =  Plate Plate

energygrp_table =  Chain Chain Chain Plate

nstlist =  1

ns_type =  simple

rlist   =  0.4

pbc =  xyz

vdw-type=  User

rvdw=  0.4

coulombtype =  cut-off

rcoulomb=  0.1

Tcoupl  =  nose-hoover

tc-grps =  System

tau_t   =  0.1

ref_t   =  450

gen_vel =  yes

gen_temp=  450

constraints =  none

pcoupl  =  no

 

After launching this job the process is running but there is no trajectory
to be produced. The log file finishes in the section containing following
information

 

   energygrp_flags[  0]: 2 2

   energygrp_flags[  1]: 2 1

 

I have tried also running this job without user specified potential, with
relaxed plate, with pbc=no and some other attempts, but I always failed in
the same way. I would be very thankful for suggestions how to cope with this
problem. 

 

Greetings

 

Zuzana

 

 

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

RE: [gmx-users] running a simulation without production

2008-11-04 Thread Zuzana Benkova
Dear Mark, 
Thank you very much for valuable advice. 
Previously, I run some simulations and I obtained what I had expected. For
instance, instead of the plate I have just used one single chain which was
treated as frozen and the simulations were running. I obtained all output
files containing appropriate information. 
I have GMX version 3.3.3. I have just tried gmxtest.pl test first in sthe
ingle precision. The output I got was

All 16 simple tests PASSED
FAILED. Check files in dec+water
FAILED. Check files in fe_test
2 out of 14 complex tests FAILED
All 63 kernel tests PASSED
All 45 pdb2gmx tests PASSED

I am not very familiar with programming language C++. As I wanted to test
also the double precision I did some changes in the gmxtest.pl file. I am
not sure if I did it correctly. I have installed GMX with _mpi suffix, so I
have changed the sub setup_vars section 

sub setup_vars {
my $mdprefix = "";
if ( $parallel > 0 ) {
$mdprefix = "mpirun -c $parallel ";
}

if ( $double > 0 ) {
$grompp = "grompp_mpi";
$mdrun  = $mdprefix . "mdrun_mpi";
$ref= "reference_d";
}

and after applying command gmxtest.pl -double all I have obtained 

sh: line 1: 30240 Segmentation fault  mdrun_mpi >mdrun.out 2>&1
FAILED. Check files in rb1
1 out of 16 simple tests FAILED
All 14 complex tests PASSED
All 63 kernel tests PASSED
All 45 pdb2gmx tests PASSED

In the case of the single precision test I can not figure out the error
since I do not have any new files produced after the test in corresponding
directories. In the case of double precision test I have also log file in
rb1 directory which has been terminated exactly in the same way as my
attempts to run simulation of plate-chain system. 

Please, could you give me some explanation and hints how to proceed further?
It would be very helpful for me. I have checked the mailing list but did not
find something that could help me. Thank you in advance.

Zuzana


-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Mark Abraham
Sent: Monday, November 03, 2008 9:55 PM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] running a simulation without production

Zuzana Benkova wrote:
> I have tried also running this job without user specified potential, 
> with relaxed plate, with pbc=no and some other attempts, but I always 
> failed in the same way. I would be very thankful for suggestions how to 
> cope with this problem.

Do other simulations run? Try one of the tutorials or gmxtest.pl.

Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] running a simulation without production

2008-11-03 Thread Zuzana Benkova
Dear GROMACS users,

 

I want to simulate a system consisting of a plate modeled by atoms (5000
atoms) and a polyethylene chain (80 carbons). The plate is frozen and I am
using tabulated potentials. In this system, the considered interactions are
chain-chain and chain-plate. There are no electrostatic interactions. The
system is placed in a box, with dimensions corresponding to the dimensions
of plate and substantially larger height. My md.mdp file is as follows

 

cpp =  cpp

define  =

 

integrator  =  md

dt  =  0.002

nsteps  =  7500

nstcomm =  1

comm_mode   =  None

nstxout =  10

nstvout =  10

nstfout =  10

nstxtcout   =  10

nstlog  =  10

nstenergy   =  10

energygrps  =  Chain Plate

freezegrps  =  Plate

freezedim   =  Y Y Y

energygrp_excl  =  Plate Plate

energygrp_table =  Chain Chain Chain Plate

nstlist =  1

ns_type =  simple

rlist   =  0.4

pbc =  xyz

vdw-type=  User

rvdw=  0.4

coulombtype =  cut-off

rcoulomb=  0.1

Tcoupl  =  nose-hoover

tc-grps =  System

tau_t   =  0.1

ref_t   =  450

gen_vel =  yes

gen_temp=  450

constraints =  none

pcoupl  =  no

 

After launching this job the process is running but there is no trajectory
to be produced. The log file finishes in the section containing following
information

 

   energygrp_flags[  0]: 2 2

   energygrp_flags[  1]: 2 1

 

I have tried also running this job without user specified potential, with
relaxed plate, with pbc=no and some other attempts, but I always failed in
the same way. I would be very thankful for suggestions how to cope with this
problem. 

 

Greetings

 

Zuzana

 

 

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] g_angle

2008-04-09 Thread Zuzana Benkova

Dear GMX users,
I want to analyse th ebond angle. I have used

mk_angndx -s md1.tpr -n angle.ndx

to prepare an index for g_angle.

After
  g_angle -f md1.xtc -s md1.tpr -n angle.ndx

I have obtained

Software inconsistency error:
Not supported in get_stx_coordnum
  Most likely this is a trivial question and I did some stupid mistake
but I would be thankful for help. Thank you in advance
Zuzana

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] angle analysis

2008-04-04 Thread Zuzana Benkova

Dear GMX users,
I want to analyse th ebond angle. I have used

mk_angndx -s md1.tpr -n angle.ndx

to prepare an index for g_angle.

After
 g_angle -f md1.xtc -s md1.tpr -n angle.ndx

I have obtained

Software inconsistency error:
Not supported in get_stx_coordnum
 Most likely this is a trivial question and I did some stupid mistake  
but I would be thankful for help. Thank you in advance

Zuzana
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] problem in hdb file

2008-02-19 Thread Zuzana Benkova

Dear GMX users,
I want to prepare .gro and topology file starting from the .pdb file  
for a polyethylene chain , whose residue is not defined in .rtp, .hdb.  
and .tdb files. I have defined a new residue in ffoplsaa.rtp file as  
follows


[ Eth ]
 [ atoms ]
   C1opls_136-0.1201
   H11   opls_140 0.0601
   H12   opls_140 0.0601
   C2opls_136-0.1202
   H21   opls_140 0.0602
   H22   opls_140 0.0602
 [ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
   C2+C1
 In ffoplsaa.hdb I defined it as

Eth 2
2   6   H1   C1  C2 +C1
2   6   H2   C2  C1 -C2

At this point I have to note that the presence of the last control  
atoms in this file has not influenced the process.

In ffoplsaa-n.tdb file I have added

[ Eth ]
[ replace ]
C2  opls_13512.011   -0.18
[ add ]
34  H2  C2  C1
 opls_140   1.008   0.06
[ delete ]
H21
H22

and in ffoplsaa-c.tdb file I have added

[ Eth ]
[ replace ]
C1  opls_13512.011  -0.18
[ add ]
34  H1  C1  C2
 opls_140   1.008   0.06
[ delete ]
H21
H22


After using pdb2gmx -f input.pdb I have obtained following message

Opening library file ffoplsaa.rtp
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/aminoacids.dat

Reading hexane.pdb...
Read 6 atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/xlateat.dat

26 out of 26 lines of xlateat.dat converted succesfully
Analyzing pdb file
There are 1 chains and 0 blocks of water and 3 residues with 6 atoms
  chain  #res #atoms
  1 ' ' 3  6
All occupancies are one
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/ffoplsaa.atp

Atomtype 817
Reading residue database... (ffoplsaa)
Opening library file ffoplsaa.rtp
Residue 59
Sorting it all out...
Opening library file ffoplsaa.hdb
Error in hdb file: nah = 41
line = ''
Opening library file ffoplsaa-n.tdb
Opening library file ffoplsaa-c.tdb
Back Off! I just backed up topol.top to ./#topol.top.8#
Processing chain 1 (6 atoms, 3 residues)
There are 0 donors and 0 acceptors
There are 0 hydrogen bonds
Checking for duplicate atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/specbond.dat

5 out of 5 lines of specbond.dat converted succesfully
N-terminus: Eth
C-terminus: Eth
Segmentation fault

Could you specify what means Error in hdb file: nah = 41
line = '' and what is the problem with segmentation.

My input structure in .pdb file is as follows

ATOM  1  C1  Eth 1   1.000   1.540   0.000
ATOM  2  C2  Eth 1   2.456   2.041   0.000
ATOM  3  C1  Eth 2   2.456   3.581   0.000
ATOM  4  C2  Eth 2   3.912   4.083   0.000
ATOM  5  C1  Eth 3   3.912   5.623   0.000
ATOM  6  C2  Eth 3   5.368   6.124   0.000
END

Thank you for help in advance
Best regards
Zuzana


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Polyethylene

2008-02-18 Thread Zuzana Benkova

Dear GMX users,
I want to prepare .gro and topology file starting from the .pdb file  
for a polyethylene chain , whose residue is not defined in .rtp, .hdb.  
and .tdb files. I have defined a new residue in ffoplsaa.rtp file as  
follows


[ Eth ]
 [ atoms ]
   C1opls_136-0.1201
   H11   opls_140 0.0601
   H12   opls_140 0.0601
   C2opls_136-0.1202
   H21   opls_140 0.0602
   H22   opls_140 0.0602
 [ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
   C2+C1
 In ffoplsaa.hdb I defined it as

Eth 2
2   6   H1   C1  C2 +C1
2   6   H2   C2  C1 -C2

At this point I have to note that the presence of the last control  
atoms in this file has not influenced the process.

In ffoplsaa-n.tdb file I have added

[ Eth ]
[ replace ]
C2  opls_13512.011   -0.18
[ add ]
34  H2  C2  C1
 opls_140   1.008   0.06
[ delete ]
H21
H22

and in ffoplsaa-c.tdb file I have added

[ Eth ]
[ replace ]
C1  opls_13512.011  -0.18
[ add ]
34  H1  C1  C2
 opls_140   1.008   0.06
[ delete ]
H21
H22


After using pdb2gmx -f input.pdb I have obtained following message

Opening library file ffoplsaa.rtp
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/aminoacids.dat

Reading hexane.pdb...
Read 6 atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/xlateat.dat

26 out of 26 lines of xlateat.dat converted succesfully
Analyzing pdb file
There are 1 chains and 0 blocks of water and 3 residues with 6 atoms
  chain  #res #atoms
  1 ' ' 3  6
All occupancies are one
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/ffoplsaa.atp

Atomtype 817
Reading residue database... (ffoplsaa)
Opening library file ffoplsaa.rtp
Residue 59
Sorting it all out...
Opening library file ffoplsaa.hdb
Error in hdb file: nah = 41
line = ''
Opening library file ffoplsaa-n.tdb
Opening library file ffoplsaa-c.tdb
Back Off! I just backed up topol.top to ./#topol.top.8#
Processing chain 1 (6 atoms, 3 residues)
There are 0 donors and 0 acceptors
There are 0 hydrogen bonds
Checking for duplicate atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/specbond.dat

5 out of 5 lines of specbond.dat converted succesfully
N-terminus: Eth
C-terminus: Eth
Segmentation fault

Could you specify what means Error in hdb file: nah = 41
line = '' and what is the problem with segmentation.

My input structure in .pdb file is as follows

ATOM  1  C1  Eth 1   1.000   1.540   0.000
ATOM  2  C2  Eth 1   2.456   2.041   0.000
ATOM  3  C1  Eth 2   2.456   3.581   0.000
ATOM  4  C2  Eth 2   3.912   4.083   0.000
ATOM  5  C1  Eth 3   3.912   5.623   0.000
ATOM  6  C2  Eth 3   5.368   6.124   0.000
END

Thank you for help in advance
Best regards
Zuzana


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Polyethylene

2008-02-17 Thread Zuzana Benkova

Dear GMX users,
I want to prepare .gro and topology files starting from the .pdb file  
for a polyethylene chain , whose residue is not defined in .rtp, .hdb.  
and .tdb files. I have defined a new residue in ffoplsaa.rtp file as  
follows


[ Eth ]
 [ atoms ]
   C1opls_136-0.1201
   H11   opls_140 0.0601
   H12   opls_140 0.0601
   C2opls_136-0.1202
   H21   opls_140 0.0602
   H22   opls_140 0.0602
 [ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
   C2+C1
 In ffoplsaa.hdb I defined it as

Eth 2
2   6   H1   C1  C2 +C1
2   6   H2   C2  C1 -C2

At this point I have to note that the presence of the last control  
atoms in this file has not influenced the process.

In ffoplsaa-n.tdb file I have added

[ Eth ]
[ replace ]
C2  opls_13512.011   -0.18
[ add ]
34  H2  C2  C1
 opls_140   1.008   0.06
[ delete ]
H21
H22

and in ffoplsaa-c.tdb file I have added

[ Eth ]
[ replace ]
C1  opls_13512.011  -0.18
[ add ]
34  H1  C1  C2
 opls_140   1.008   0.06
[ delete ]
H11
H12

After using pdb2gmx -f input.pdb I have obtained following message

Opening library file ffoplsaa.rtp
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/aminoacids.dat

Reading hexane.pdb...
Read 6 atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/xlateat.dat

26 out of 26 lines of xlateat.dat converted succesfully
Analyzing pdb file
There are 1 chains and 0 blocks of water and 3 residues with 6 atoms
  chain  #res #atoms
  1 ' ' 3  6
All occupancies are one
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/ffoplsaa.atp

Atomtype 817
Reading residue database... (ffoplsaa)
Opening library file ffoplsaa.rtp
Residue 59
Sorting it all out...
Opening library file ffoplsaa.hdb
Error in hdb file: nah = 41
line = ''
Opening library file ffoplsaa-n.tdb
Opening library file ffoplsaa-c.tdb
Back Off! I just backed up topol.top to ./#topol.top.8#
Processing chain 1 (6 atoms, 3 residues)
There are 0 donors and 0 acceptors
There are 0 hydrogen bonds
Checking for duplicate atoms
Opening library file  
/usr/local/share/gromacs/gromacs-3.3/share/top/specbond.dat

5 out of 5 lines of specbond.dat converted succesfully
N-terminus: Eth
C-terminus: Eth
Segmentation fault

Could you specify what means Error in hdb file: nah = 41
line = '' and what is the problem with segmentation.

My input structure in .pdb file is as follows

ATOM  1  C1  Eth 1   1.000   1.540   0.000
ATOM  2  C2  Eth 1   2.456   2.041   0.000
ATOM  3  C1  Eth 2   2.456   3.581   0.000
ATOM  4  C2  Eth 2   3.912   4.083   0.000
ATOM  5  C1  Eth 3   3.912   5.623   0.000
ATOM  6  C2  Eth 3   5.368   6.124   0.000
END

Thank you for help in advance
Best regards
Zuzana

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


RE: [gmx-users] Temperature increase

2007-05-04 Thread Zuzana Benkova
Dear Mark, thank you very much for prompt reply, I am sorry but I did
mistake in my first e-mail I have run an NVE simulation of decamer of glycin
first with PBC (coulombtype=PME)and later without PBC (coulombtype=Cut-off).
Zuzana

-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Mark Abraham
Sent: Thursday, May 03, 2007 3:46 PM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] Temperature increase

Zuzana Benkova wrote:
> I have run an NVE simulation of decamer of glycin first with PBC 
> (coulombtype=Cut-off) and later without PBC (coulombtype=PME).

You can't have a non-periodic system with PME.

> The 
> gen-temp was set to 300K with option gen-vel=yes. The temperature start 
> to raise during simulation and finaly the following messages appeared in 
> log file
> 
> Shake did not converge in 1000 steps
> 
> i mi  j mj  before   after   should be
> 
>66  14.01 67   1.01 0.10100 0.10318 0.10100
> 
> Constraint error in algorithm Shake at step 980454
> 
> Wrote pdb files with previous and current coordinates
> 
> How can I cope with this?

We can't tell, because the above description is not self-consistent.

In any case you haven't defined your topology, or minimized your 
structure or equilibrated properly. Then some bad structural element is 
causing a region to have too much energy which is breaking structures. 
Search the mailing list archives for suggested equilibration regimes, or 
read relevant literature examples, or do tutorial material.

Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Temperature increase

2007-05-03 Thread Zuzana Benkova
I have run an NVE simulation of decamer of glycin first with PBC
(coulombtype=Cut-off) and later without PBC (coulombtype=PME). The gen-temp
was set to 300K with option gen-vel=yes. The temperature start to raise
during simulation and finaly the following messages appeared in log file

Shake did not converge in 1000 steps

i mi  j mj  before   after   should be

   66  14.01 67   1.01 0.10100 0.10318 0.10100

Constraint error in algorithm Shake at step 980454

Wrote pdb files with previous and current coordinates

How can I cope with this?

Tahnk you in advance

Zuzana

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] (no subject)

2007-05-03 Thread Zuzana Benkova
test

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] polyethelenes

2006-07-04 Thread Zuzana Benkova








Hello gmx-users,

I have only little experiences with Gromacs and I am
trying to run MD simulations on polyethelenes. I tried to modify the
ffoplsaa.rtp and ffoplsaa.hdb where I have added my own residui. The addition
in ffoplsaa.rtp is as follows

 

; Ethyle - CH3-CH2 the first unit in polyethylenes

[ TTH ]

 [ atoms ]

    C1  
opls_135    -0.180    1

   HA1  
opls_140 0.060    1

   HA2  
opls_140 0.060    1

   HA3  
opls_140 0.060    1

    C2  
opls_136    -0.120    2

   HB1  
opls_140 0.060    2

   HB2  
opls_140 0.060    2

 [ bonds ]

   C1    HA1

   C1    HA2

   C1    HA3

   C1    C2

   C2    HB1

   C2    HB2

 

; Ethylene - CH2-CH2 central parts of polyetylens
(CnHm)

[ ETH ]

 [ atoms ]

    C1  
opls_135    -0.180    1

   HA1  
opls_140 0.060    1

   HA2  
opls_140 0.060    1

    C2  
opls_136    -0.120    2

   HB1  
opls_140 0.060    2

   HB2  
opls_140 0.060    2

 [ bonds ]

   C1    HA1

   C1    HA2

   C1    C2

   C2    HB1

   C2    HB2

 

; Ethyle - CH2-CH3 the last unit in polyethylenes

[ THT ]

 [ atoms ]

    C1  
opls_135    -0.180    1

   HA1  
opls_140 0.060    1

   HA2  
opls_140 0.060    1

    C2  
opls_136    -0.120    2

   HB1  
opls_140 0.060    2

   HB2  
opls_140 0.060    2

   HB3  
opls_140 0.060    2

 [ bonds ]

   C1    HA1

   C1    HA2

   C1    C2

   C2    HB1

   C2    HB2

   C2    HB3

and the addition in ffoplsaa.hdb file is as follows

 

TTH 2

3  
4   HA 
C1  C2
HB +C2

2  
6   HB 
C2  C1 +C1

ETH 2

2  
6   HA  
C1  -C2    C2

2  
6   HB  
C2  C1    +C1

THT 2

2  
6   HA  
C1  -C2    C2

3  
4   HB  
C2  C1
-C2    HA

 

I have also added these three residui (TTH, ETH, and
THT) into aminoacids.dat file. 

 

 

My attempt to convert the pdb file into the gro and
top files using pdb2gmx I got following error messages 

 

Error in hdb file: nah = 38

line = ''

 

and

 

Fatal error:



Atom
HB not found in residue TTH1 while adding hydrogens

The
pdb file is attached in this mail. I would be thankful for any help. Thank you
in advance.

Greetings

Zuzana










Polyethylene.pdb
Description: Binary data
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php