[gmx-users] Re: np command with GROMACS 4.0.4

2009-04-22 Thread Vitaly V. Chaban
If seems in the gmx versions above 4.0 '-np' is not used. If I am not mistaken you should just point out the number of nodes in your queueing system while submitting the job. Vitaly > I am trying to do my position restrained dynamic simulation on GROMACS > 4.0.4, and I want to use 4 nodes on th

[gmx-users] RE: transport properties

2009-04-28 Thread Vitaly V. Chaban
Please read about obtaining diffusion constant from VAC in the gromacs wiki. I wrote an individual article on that topic. Best, Vitaly Subject: RE: [gmx-users] transport properties > To: Discussion list for GROMACS users > Message-ID: <569363.99887...@web36304.mail.mud.yahoo.com> >

[gmx-users] Re: How to freeze the atoms?

2009-06-23 Thread Vitaly V. Chaban
nged. >   How can I really fix the atoms? Look forward > for your suggestions! > All the atoms which relate to your frozen groups should be non-movable during MD. >the coordinates of C and O in carbonyl groups changed. Only ones of CHO changed? And the CNT itself is nonmovable

Re: 回复:[gmx-users] Re: How to freeze the atoms ?

2009-06-23 Thread Vitaly V. Chaban
le atoms in > carbonyl groups take charges. > > Thank you for your time! > > > - 原始邮件 - > 发件人:Vitaly V. Chaban > 收件人:gmx-us...@gromacs.org > 主题:[gmx-users] Re: How to freeze the atoms? > 日期:2009-6-23 18:04:35 > > > > > Dear all > >

Re: 回复:Re: 回复:[gmx-users] Re: How to fre eze the atoms?

2009-06-24 Thread Vitaly V. Chaban
Car" and the group name of carbonyl in >> the ndx file is "Car", too. Is there anything wrong with the groups? >> >> Is it possible that the frozen atoms move a short distance due to the >> strong electrostatic forces. C in the CNT has no charge, while atoms

[gmx-users] Re: dmso topology

2009-06-24 Thread Vitaly V. Chaban
I used the DMSO model from dmso.itp for some test runs with the gromacs versions before 4.0 and it gave me the same diffusion constant (~1.1) as reported in the original paper. So, the topology seems to be OK in general. Please also see the paper of Skaf (it seems of 1997) for DMSO force field...

[gmx-users] Re: How to freeze the atoms?

2009-06-25 Thread Vitaly V. Chaban
I do not see topology file in your archive. 2009/6/25 > The input files for EM of the job are in the attachment. > > > > Thank you for your time! > > > - 原始邮件 - > 发件人:Vitaly V. Chaban > 收件人:toby10222...@sina.com, gmx-users@gromacs.org > 主题:Re: 回复:Re:

[gmx-users] C6 constant units

2009-06-26 Thread Vitaly V. Chaban
Hi, In what units C6 constants are present in the force field files? Thank you, Vitaly ___ 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 bef

[gmx-users] Re: gmx-users Digest, Vol 62, Issue 158

2009-06-30 Thread Vitaly V. Chaban
example of them? Thanks. > > The examples are in the topology folder of the gromacs installation. As for me they are quite self-explanatory. -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@un

[gmx-users] RE: How to freeze the atoms?

2009-07-05 Thread Vitaly V. Chaban
D for some time to relax the system completely. -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, mob.: +38-097-8259698 Hi > > > > I

[gmx-users] xdr library. Just for an example.

2009-07-07 Thread Vitaly V. Chaban
Hi, Can anybody share any working code which uses xdr-1.1 library? - Just to consult. Thank you very much in advance, Vitaly ___ gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archiv

[gmx-users] Re: number of coordinates in coordinate file does not match topology

2009-07-16 Thread Vitaly V. Chaban
his error before adding > ions. > > Can anybody help me regarding same and suggest what could be possible > reason for this ? > > Thanks, > > Regards, > I guess you just missed five atoms in your gro file. Does the number of solvent molecules in the coordinate file equal

[gmx-users] Re: number of coordinates in coordinate file does not match topology

2009-07-16 Thread Vitaly V. Chaban
ve this error.I am getting this error before adding > >> ions. > >> > >> Can anybody help me regarding same and suggest what could be possible > >> reason for this ? > >> > >> Thanks, > >> > >> Regards, &

[gmx-users] Re: number of coordinates in coordinate file does not match topology

2009-07-16 Thread Vitaly V. Chaban
> -- > >> >> > >> >> I checked gromacs wiki page > >> >> also( > >> >> > >> > http://oldwiki.gromacs.org/index.php/Errors#Number_of_coordinates_

[gmx-users] Re: installing grace

2009-07-20 Thread Vitaly V. Chaban
In what way do you install this software: from source or from packages? Do you have the appropriate variable telling the grace where the lesstiff is installed? If you install lesstiff from the package it should have a name like "lesstiff-dev" (not just lesstiff). ~ Vitaly -- Vitaly

[gmx-users] Re: installing grace

2009-07-20 Thread Vitaly V. Chaban
Hello Shamik, Why not to try from rpm packages? It should be easier especially if this is the fisrt time of doint that. ~ Vitaly -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua

[gmx-users] Re: Groups based on distance criteria

2009-07-20 Thread Vitaly V. Chaban
V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University, Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, mob.: +38-097-8259698 http://www-rmn.univer.kharkov.ua/chaban.html

Re: [gmx-users] Re: Groups based on distance criteria

2009-07-20 Thread Vitaly V. Chaban
Hi Vasya, Glad to hear from you here! I did not work with -rerun yet - it's a useful option I see now. On Mon, Jul 20, 2009 at 5:39 PM, Vasilii Artyukhov wrote: > > > 2009/7/20 Vitaly V. Chaban > >> Hi Matteus, >> >> Yeah, dynamic groups would be a po

[gmx-users] Re: comparison of gromacs in CENTOS and UBUNTU

2009-07-21 Thread Vitaly V. Chaban
senseless. The linux distribution will not affect gromacs performance. -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, mob.: +38-0

[gmx-users] Re: PDB and GRO files

2009-07-21 Thread Vitaly V. Chaban
What about .gro files? > > I appreciate your help for solve my basic problems. > > Thanks, Jamie > Jamie, Please see the help for genbox, genconf and editconf utlities of the gromacs package. It seems you should be interested in Protein Data Bank as well: http://www.rcsb.org/pdb/h

[gmx-users] Re: PDB and GRO files

2009-07-21 Thread Vitaly V. Chaban
etc). 3. The nanotubes files (coordinates and topologies) as well as a number of already equilibrated systems were uploaded by me to the topology section of the old site. I guess they are still accessible. You may also want to look for Christopher Stiles' wonderful page with many details.

[gmx-users] Re: reliability of NOTES

2009-07-21 Thread Vitaly V. Chaban
when I changed the thermostat to nose-hoover, it said that it can not > be used with Berendsen barostat. Why? > > C'est la vie... -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@un

Re: [gmx-users] PDB and GRO files

2009-07-21 Thread Vitaly V. Chaban
Yeah, we are friendly. :) Don't hesitate to ask if you experience future problems. ~ Vitaly On Wed, Jul 22, 2009 at 1:55 AM, Jamie Seyed wrote: > Thank you for your help and comments. > -Jamie > > > > ___ gmx-users mailing listgmx-users@gromacs.

[gmx-users] Re: gmx-users Digest, Vol 63, Issue 110

2009-07-21 Thread Vitaly V. Chaban
comes out that non-running windows interfere gromacs run. -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, mob.: +38-097-8259698 http://www-rmn.un

[gmx-users] Re: genconf

2009-07-22 Thread Vitaly V. Chaban
file > to make sure no particles are overlapping with the original box? Or no > particle is coming from the new one in the old one? > > > > Payman > > > > *From:* gmx-users-boun...@gromacs.org [mailto: > gmx-users-boun...@gromacs.org] *On Behalf Of *Vitaly V. Chaban >

[gmx-users] Re: genconf

2009-07-22 Thread Vitaly V. Chaban
om > becomes infinity. This will cause the simulation to crash! > > > > Payman > > > > *From:* Vitaly V. Chaban [mailto:vvcha...@gmail.com] > *Sent:* July 22, 2009 2:20 AM > *To:* Payman Pirzadeh; gmx-users@gromacs.org > *Subject:* Re: genconf > > > >

[gmx-users] Re: trr and xtc

2009-07-22 Thread Vitaly V. Chaban
Jamie, To have a trajectory in TRR you should just change > > nstxout=0 > nstvout=0 > > to non-zeros. -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@

[gmx-users] Users software contribution. Old site.

2009-07-29 Thread Vitaly V. Chaban
Hi, I can't download the software form the users contributions section located on the old gromacs website (http://oldwww.gromacs.org/component/option,com_docman/task,cat_view/gid,24/dir,DESC/order,name/limit,5/limitstart,35/). As pressing "download" it seems to open binary file in the browser wind

[gmx-users] Re: Users software contribution. Old site.

2009-07-29 Thread Vitaly V. Chaban
Thanks, Justin. I succeeded eventually with your link. Vitaly On Wed, Jul 29, 2009 at 11:54 PM, Justin A. Lemkul wrote: > > > Vitaly V. Chaban wrote: >> >> Hi, >> >> I can't download the software form the users contributions section >> lo

[gmx-users] Re: CNT

2009-07-30 Thread Vitaly V. Chaban
a.n2t > > WARNING: all CONECT records are ignored > > Opening library file > /homea/scocchi/GROMACS/share/gromacs/top/ffoplsaa.n2t > > There are 23 name to type translations > > Generating bonds from distances... > > atom 0 > > > > ... and at this point it

[gmx-users] Re: CNT

2009-07-30 Thread Vitaly V. Chaban
ogies > www.dti.supsi.ch > SUPSI - University of Applied Science of Southern Switzerland > www.supsi.ch > > Via Cantonale - Galleria 2 > 6928 Manno > Switzerland > tel: +41 (0)58 666 66 46 > fax: +41 (0)58 666 66 20 > > > -Messaggio originale- > Da: Vi

[gmx-users] Re: R: CNT

2009-07-30 Thread Vitaly V. Chaban
icimsi.ch > DTI - Department of Innovative Technologies > www.dti.supsi.ch > SUPSI - University of Applied Science of Southern Switzerland > www.supsi.ch > > Via Cantonale - Galleria 2 > 6928 Manno > Switzerland > tel: +41 (0)58 666 66 46 > fax: +41 (0)58 666 66 20

Re: [gmx-users] Re: CNT

2009-07-30 Thread Vitaly V. Chaban
boundary conditions for infinite molecules > seem to be handled correctly... > > 2009/7/30 Vitaly V. Chaban >> >> Giulio, >> >> I suggest you to try x2top which is spread with gromacs-3.3.1. The >> newer versions of x2top are not sain. >> >> Vi

[gmx-users] Re: R: R: CNT

2009-07-30 Thread Vitaly V. Chaban
rsity of Applied Science of Southern Switzerland > www.supsi.ch > > Via Cantonale - Galleria 2 > 6928 Manno > Switzerland > tel: +41 (0)58 666 66 46 > fax: +41 (0)58 666 66 20 > > -Messaggio originale- > Da: Vitaly V. Chaban [mailto:vvcha...@gmail.com] > Inviato: gio

[gmx-users] Re: R: R: R: CNT

2009-07-30 Thread Vitaly V. Chaban
If you unzipped/untarred the archive properly it should be some system already prepared to run. So you can just start grompp and mdrun to begin MD. > > By the way, I had already seen that web page... I downloaded and unzipped the > swcnt file, but I can't understand what kind of file it is. > __

[gmx-users] Re: topolbuild1_2_1.tgz decompressing problem

2009-08-02 Thread Vitaly V. Chaban
> > Hello, > > I downloaded "topolbuild1_2_1.tgz" from the URL: > http://www.gromacs.org/@api/deki/files/40/=topolbuild1_2_1.tgz > It seems you are right. I tried the same and also found the archive broken. Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V

[gmx-users] Re: itp file

2009-08-08 Thread Vitaly V. Chaban
Hi, Did you try pdb2gmx, x2top to make .top and then export to .itp? Sometimes also it appears easier to make the topology by hand if your structure is not big. Please see the topologies for some molecules in the gromacs archive to use them as an example. Vitaly > > Dear gmx users, > I need some

Re: [gmx-users] Re: itp file

2009-08-09 Thread Vitaly V. Chaban
tried but I got error messages for both cases. > Actually the structure is big and contains a c60 part. Do you have any > advice? I could not find in the topology part of archive... > Thanks/Jamie > > > On Sun, Aug 9, 2009 at 2:54 AM, Vitaly V. Chaban wrote: >> >> Hi,

Re: [gmx-users] Re: itp file

2009-08-10 Thread Vitaly V. Chaban
that something > else?? I appreciate your help... > Thanks a lot/Jamie > > On Mon, Aug 10, 2009 at 2:18 AM, Vitaly V. Chaban > wrote: >> >> Jamie, >> >> What error messages do they provide? What force fields do you use? >> There were many rumors about the

Re: [gmx-users] Re: itp file

2009-08-10 Thread Vitaly V. Chaban
t; 2) with Christopher Stiles method, why I got that fatal error and how I can > solve it to run a correct md simulation?? > Many Thanks in advance/Jamie > > On Mon, Aug 10, 2009 at 10:45 AM, Vitaly V. Chaban > wrote: >> >> Jamie, >> >> The best solution is

Re: [gmx-users] Re: itp file

2009-08-10 Thread Vitaly V. Chaban
mixed the ff and I put ff.itp on the top of my topology file... > would you please let me know what is causing the problem?? I appreciate your > help.../Jamie > > On Mon, Aug 10, 2009 at 4:29 PM, Vitaly V. Chaban > wrote: >> >> Dear Jamie, >> >> > 1) when I h

Re: [gmx-users] Re: itp file

2009-08-10 Thread Vitaly V. Chaban
t;> the case I am going to undo it to see what happens... I will let you know. >> Thanks/Jamie >> On Mon, Aug 10, 2009 at 5:05 PM, Vitaly V. Chaban >> wrote: >>> >>> Why did you comment out [defaults] in ffoplsaa.itp? The general idea >>> is that the

Re: [gmx-users] Re: itp file

2009-08-10 Thread Vitaly V. Chaban
nclude "c80ch3.itp" > [system] > fullerene in water > [molecules] > C80CH3   1 > SOL   1000 > > my command is: > grompp -f md.mdp -c f.gro -p topol.top -o pr.tpr  -maxwarn 10 > > Thank you for helping me out to solve the problem/ Regard

[gmx-users] Re: A bug for MPI of GROMACS?

2009-08-11 Thread Vitaly V. Chaban
system is a polymer > (dendrimer) without pbc. I wonder if this is a bug for MPI of GROMACS. Could > you give me some hints? Thanks a lot for your attention! > -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Uk

[gmx-users] Re: xdrfile library, read co-ordinates from .xtc

2009-08-11 Thread Vitaly V. Chaban
&step_xtc, &time_xtc, box_xtc, > x_xtc, &prec_xtc); >  printf("%f %f %f\n",x_xtc[0],x_xtc[1],x_xtc[2]); > > The coordinates in the variable x_xtc seems to be random. How do I access > the proper coordinates from the trajectory? > > I would appreciate any help.

[gmx-users] Fwd: xdrfile library, read co-ordinates from .xtc

2009-08-11 Thread Vitaly V. Chaban
than me can provide such one, I will be very grateful. Thanks, Vitaly -- Forwarded message -- From: Vitaly V. Chaban Date: Tue, Aug 11, 2009 at 11:11 PM Subject: Re: xdrfile library, read co-ordinates from .xtc To: gmx-users@gromacs.org, Lipi Thukral Hi Lipi, What do you

[gmx-users] Re: [gmx-developers] Fwd: xdrfile library, read co-ordinates from .xtc

2009-08-12 Thread Vitaly V. Chaban
Janne, Jussi, Thank you very much! Vitaly On Wed, Aug 12, 2009 at 11:35 AM, Janne Blomqvist wrote: > Vitaly V. Chaban wrote: >> >> On the same topic, I have a question to the community: does anybody >> try to make *so*-library (or *dll* on windows) of xdrfile in order

[gmx-users] Re: solvation

2009-08-13 Thread Vitaly V. Chaban
Jamie, I didn't understand why you meant saying "waters far from the outside of the pore"... Far from outside, then where exactly? Maybe I just imagine your box incorrectly. To fill the pore with a solvent try to "deceive" genbox a little: 1) make a box (see editconf) in such a way that it contai

Re: [gmx-users] Re: solvation

2009-08-13 Thread Vitaly V. Chaban
> first without introducing my molecule... > I appreciate to get more information and some help... Thanks a lot/Jamie > > > On Thu, Aug 13, 2009 at 12:21 PM, Vitaly V. Chaban > wrote: >> >> Jamie, >> >> I didn't understand why you meant saying &quo

Re: [gmx-users] Re: solvation

2009-08-13 Thread Vitaly V. Chaban
     Fatal error: >       One of the box vectors has become shorter than twice the > cut-off length or box_yy-|box_zy| or box_zz has become smaller than the > cut-off. >  --- > Thanks in advance/Jamie > > > On Thu, Aug 13, 2009 at 2:10 PM, Vitaly

[gmx-users] Re: solvation

2009-08-13 Thread Vitaly V. Chaban
Tsjerk, I understood the situation in such a way that MD was not run at all, so Jamie visualizes *gro* file obtained with *genbox*... Jamie, didn't you look at this: http://www.ime.unicamp.br/~martinez/packmol/ to prepare your system? Vitaly > >> Thanks for the answer. When I said far away me

[gmx-users] Re: help needed for using xdr library

2009-08-17 Thread Vitaly V. Chaban
Hi, 1. It seems you also need to include xdrfile.h. 2. Are the library files located in the same directory with your source files when you link? VItaly > > Recently I am trying to use the xdr library downloaded from gromacs org to > read xtc files. However when I link my program with the library

[gmx-users] gmx-list digest

2009-08-20 Thread Vitaly V. Chaban
Hi all, In what way is gromacs list digest actually formed? How the system decides that the time to form and send the list digest comes? During the last days I got many digests (more than 5 ones a day) many of which contained only 1 or 2 messages, that seems quite inconvenient. Maybe it's reasona

[gmx-users] Re: gmx-users Digest, Vol 64, Issue 138

2009-08-21 Thread Vitaly V. Chaban
Jamie, How do you like g_vellacc or g_dipoles? :) Best, Vitaly > > Dear all, > Is there any tool in gromacs that can calculate time correlation > function ??? R(x,y,z). Many Thanks in Advance/Jamie ___ gmx-users mailing listgmx-users@gromacs.org

[gmx-users] Re: water simulation

2009-08-25 Thread Vitaly V. Chaban
--- > Program grompp, VERSION 4.0.4 > Source code file: topio.c, line: 415 > > Fatal error: > Syntax error - File topol.top, line 4 > Last line read: > '[ atomtypes ]' > Invalid order for directive atomtypes >

[gmx-users] Re: wierd behavior of mdrun

2009-08-26 Thread Vitaly V. Chaban
r shows that > the job is running, also the processors are 100% engaged while I have no > outputs! A standard guess: what about trying to run the single-processor job on the same cluster? Does it run OK? -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National Univer

[gmx-users] Re: wierd behavior of mdrun

2009-08-27 Thread Vitaly V. Chaban
aymon Pirzadeh wrote: > Yes, > it works when it is run on one processor interactively! > > On Thu, 2009-08-27 at 09:23 +0300, Vitaly V. Chaban wrote: >> > I made a .tpr file for my md run without any problems (using the bottom >> > mdp file). My job submission scrip

[gmx-users] Re: SIMULATION OF GOLD NANOPARTICLE

2009-08-28 Thread Vitaly V. Chaban
Hi George, I know some attempts of gold nanosize particles simulation published during the last years. Did you explore the literature? ~ Vitaly > > Hi GROMACS users > > Does anyone have any information on the simulation of gold nanoparticles with > GROMACS? > > Thank you in advance > > George _

[gmx-users] Re: gmx-users Digest, Vol 64, Issue 200

2009-08-30 Thread Vitaly V. Chaban
> Do you have any suggestions on what options to use in the configure command > when installing GROMACS for parallel run? Is it possible to change these > options at a later time? > > Thank you in advance > > George I suppose "--enable-mpi" will be enough to make

[gmx-users] Re: Binary mixture

2009-09-01 Thread Vitaly V. Chaban
Hi, Your objects look very familiar for me. :) One can't make .gro from .itp. You should first make .pdb with coordinates. By the way gromacs works with .pdb directly. Good luck, Vitaly > Thanks a lot for your reply. Actually I will tell you the > problem. I want a binary mixt

[gmx-users] Re: Binary mixture

2009-09-01 Thread Vitaly V. Chaban
So what's your goal? If you have .pdb of one molecule, it's easy to disseminate it using genconf or genbox. On Tue, Sep 1, 2009 at 3:49 PM, Abhishek Banerjee wrote: > hi, >    Actually I checked what you told me to do. I used genbox -ci DMF.pdb > -box 10 10 10 -nmol 100. Then I am getting a ra

[gmx-users] Re: Binary mixture

2009-09-01 Thread Vitaly V. Chaban
Banerjee, I understand your task but I still did not understand your problem on this way. First you should have the .pdb with methanol molecule (coordinates) and the .pdb with DMF (coordinates). Then use genconf or genbox to make a box of pure methanol (or DMF) and add another species with genbox.

[gmx-users] Re: MPI on Windows

2009-09-02 Thread Vitaly V. Chaban
It seems gromacs can be compiled on windows directly if you a bit lucky. And what's the sense of this task? Do you have a cluster running windows? Vitaly > Hi all > > I would like to ask if it is possible to use MPI when GROMACS is installed on > Windows using Cygwin. If yes, is there any spec

[gmx-users] Re: hi questions about installing GROMACS

2009-09-04 Thread Vitaly V. Chaban
fied and no makefile found. Stop." > > What should i do for this? Should i save the GROMACS Directory somewhere > else? > > Hi, What does your configure script finish with? -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National Univer

[gmx-users] Re: DMF topology file

2009-09-05 Thread Vitaly V. Chaban
automatically. The second way is to make topology by hand. If one has enough experience it can be of some pleasure. :) Good luck! Vitaly -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha

[gmx-users] Re: System has non-zero total charge: 2.570000e+00

2009-09-12 Thread Vitaly V. Chaban
rst question is: what charge (integer) of the system do you expect? May the value you report here occur due to the wrong force field you apply? -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@unive

[gmx-users] Re: System has non-zero total charge: 2.570000e+00

2009-09-13 Thread Vitaly V. Chaban
In a word, do you use a respectable FF or a instantly made one? :) ~ Vitaly On Sun, Sep 13, 2009 at 11:06 AM, Tsjerk Wassenaar wrote: > Hi, > > On Sun, Sep 13, 2009 at 8:55 AM, Vitaly V. Chaban wrote: >>> before energy minimization step , I performed the preprosessing step

[gmx-users] imprper dihedrals

2009-09-21 Thread Vitaly V. Chaban
yntax. Thanks a lot! Vitaly -- Vitaly V. Chaban, Ph.D. (ABD) School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, cell.: +38-097-8259698 http://www-rmn.univer.kharkov.ua/chaban

[gmx-users] n2t files in gromacs

2009-09-22 Thread Vitaly V. Chaban
Hi, Could anybody please suggest me what is the syntax of N2T files? The only sentence I found in the manual tells that it is "atom to atomtype translation table". What do the numbers mean? Thanks, Vitaly ___ gmx-users mailing listgmx-users@gromacs

[gmx-users] viscovity with a help of g_energy

2009-09-22 Thread Vitaly V. Chaban
Hi, I which file should one look for the value of shear viscosity as using "g_energy -vis" to calculate it? Thanks, Vitaly ___ gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive

[gmx-users] Improper dihedrals, which way is correct?

2009-09-23 Thread Vitaly V. Chaban
Hi, I have got confused with setting up improper dihedrals in the gromacs topologies. In the manual (ver. 4.0) I see the section named [ impropers ] (page 104) " [ impropers ] C O1 O2 CA " but in the example system of the peptide (share/gromacs/tutor/speptide) one finds defining impropers using t

Re: [gmx-users] viscovity with a help of g_energy

2009-09-23 Thread Vitaly V. Chaban
Justin, It looks especially strange for me that all these files contain only the columns of data but not the resulting shear viscosity value. Vitaly On Tue, Sep 22, 2009 at 8:26 PM, Justin A. Lemkul wrote: > > > Vitaly V. Chaban wrote: >> >> Hi, >> >> I which f

[gmx-users] Re: Improper dihedrals, which way is correct?

2009-09-23 Thread Vitaly V. Chaban
l.top: [ dihedrals ] atoms func angle constant 1 2 6 7 2180 4.602 Will this line do what I want here? Vitaly On Wed, Sep 23, 2009 at 1:53 PM, Vitaly V. Chaban wrote: > Hi, > > I have got confused with setting up improper dihedrals in the gromacs > topologies. In the ma

[gmx-users] Re: Improper dihedrals, which way is correct?

2009-09-23 Thread Vitaly V. Chaban
Justin, I think force field specific syntax is important for the utilities like pdb2gmx or x2top but not for grompp which works with topol.top. Here, I do not use any "standard" force field but define everything myself from the very beginning. I believe this format (in topol.top) should be FF-inde

[gmx-users] viscosity and acflen in g_energy

2009-09-25 Thread Vitaly V. Chaban
Hi, When using g_energy to calculate viscosity ("g_energy -vis"), in enecorr.xvg the ACF of stress tensor is output. How can one specify the "acflen" of this ACF? I try the key "g_energy -vis -acflen XXX" but the length of the ACF is still equal to the half of the trajectory length. Thanks, Vital

[gmx-users] 1-4 interactions

2009-09-25 Thread Vitaly V. Chaban
Hi, Is it possible in gromacs to make 1-4 interaction energy reduced as compared to the intermolecular interations between the same centers? In other words, I want the INTRAmolecular interaction A-B to be by some factor smaller than INTERmolecular A-B interactions in the box. Thanks, Vitaly _

[gmx-users] Re: viscosity and acflen in g_energy

2009-09-26 Thread Vitaly V. Chaban
Mark, In "g_energy -h" I see nothing abous viscosity except the option "-vis". The command line is as I wrote earlier: "g_energy -vis -acflen 1001" The system is liquid water. The simulation time is 2 000 ps. The file I look into is enecorr.xvg. What the error of some sort do you mean? Vitaly

[gmx-users] Re: viscosity and acflen in g_energy

2009-09-28 Thread Vitaly V. Chaban
> Viscosity in g_energy is computed from the pressure tensor which is > stored in the energy file. Unfortunately is seems like the acflen is > ignored in this case. Please submit a bugzilla. Done. > On the other hand you need all the data you can get for this quantity, If the system is very big

[gmx-users] Re: viscosity and acflen in g_energy

2009-09-28 Thread Vitaly V. Chaban
> Viscosity in g_energy is computed from the pressure tensor which is > stored in the energy file. Unfortunately is seems like the acflen is > ignored in this case. Please submit a bugzilla. Done. > On the other hand you need all the data you can get for this quantity, > because it converges extr

[gmx-users] Re: Generating a .gro file from .pdb file

2009-09-28 Thread Vitaly V. Chaban
Hi, editconf -f waters.pdb -o conf.pdb As for CHARMM, you may select the atoms of interest and incorporate to your topology (see chapter 5 of the manual) Also maybe some script for the conversion exist, try to search. Vitaly > > Hello everyone, >  Thanks to Justin Lemkul and Thomas Schlesier fo

[gmx-users] Re: Simulation of a protein in a water/DMSO box

2009-09-28 Thread Vitaly V. Chaban
Hi Simone, I would advise you to insert the needed number of waters (genbox -cp conf.gro -cs water.gro -o conf1.gro) and then do the same with DMSO ((genbox -cp conf1.gro -cs dmso.gro -o conf2.gro). Then just equilibrate. Vitaly > > I need to run a simulation of a protein in a water/DMSO mixtur

Re: [gmx-users] Re: Simulation of a protein in a water/DMSO box

2009-09-28 Thread Vitaly V. Chaban
Justin, Simeone will insert only "the needed number of waters" but not the full box. :) To do so he will use an option "-maxsol XXX" of genbox. Vitaly On Mon, Sep 28, 2009 at 5:08 PM, Justin A. Lemkul wrote: > > > Vitaly V. Chaban wrote: >> >> Hi Simon

[gmx-users] Re: g_energy and viscosity

2009-10-01 Thread Vitaly V. Chaban
> Viscosity in g_energy is computed from the pressure tensor which is > stored in the energy file. Unfortunately is seems like the acflen is > ignored in this case. Please submit a bugzilla. > > On the other hand you need all the data you can get for this quantity, > because it converges extremely

[gmx-users] Re: Advantage of NPT over NVT?

2009-10-14 Thread Vitaly V. Chaban
If your model gives a correct fluid density with NPT, use NPT. If it does not, then use NVT. This is a quick practical guide, not not read much. :) Vitaly ___ gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-u

[gmx-users] Re: Advantage of NPT over NVT?

2009-10-14 Thread Vitaly V. Chaban
> Vitaly V. Chaban wrote: >> If your model gives a correct fluid density with NPT, use NPT. If it >> does not, then use NVT. >> This is a quick practical guide, not not read much. :) > > Sounds dangerous :-) If the model produces a wrong density, you'd want &

[gmx-users] Re: Creating topology for a small organic molecule

2009-10-25 Thread Vitaly V. Chaban
would advise to start with X2TOP and CHAPTER 5 of the manual. Vitaly -- Vitaly V. Chaban, Ph.D. School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, cell.: +38-097-825

[gmx-users] Polarizable models

2009-10-29 Thread Vitaly V. Chaban
-- Vitaly V. Chaban, Ph.D. School of Chemistry V.N. Karazin Kharkiv National University Svoboda sq.,4, Kharkiv 61077, Ukraine email: cha...@univer.kharkov.ua,vvcha...@gmail.com skype: vvchaban, cell.: +38-097-8259698 http://www-rmn.univer.kharkov.ua/chaban.html === !

[gmx-users] Re: Number of solvent molecules

2009-10-30 Thread Vitaly V. Chaban
> Dear users, > > How is it possible to fix the number of solvent molecules and the size of the > box used so that I can maintain exactly the same environment for different > runs. > For eg, if I have 84000 water solvent molecules and the dimensions of 6.7 > 6.73 4.71 nm; how can I use the same

[gmx-users] Shear viscosity

2009-10-30 Thread Vitaly V. Chaban
Hi all, I am trying to master the art of shear viscosity calculation using "g_energy -vis". Plotting the first column in visco.xvg over time, the nearly straight line is obtained. Is it OK? How can one get the viscosity value from these data? My simulation details for this task are: NVT, 298K, shi

Re: [gmx-users] Polarizable models

2009-10-31 Thread Vitaly V. Chaban
On Thu, Oct 29, 2009 at 1:14 PM, David van der Spoel wrote: > Vitaly V. Chaban wrote: >> >> Hi, >> >> Is there any place to read how gromacs treats polarizable models? >> The only thing I found in the official manual is that shells are >> denoted as

Fwd: [gmx-users] Polarizable models

2009-10-31 Thread Vitaly V. Chaban
And one more question: is this the only "polarizability scheme" currently incorporated in gromacs? Is there any place to read about it except the source code? On Thu, Oct 29, 2009 at 1:14 PM, David van der Spoel wrote: > Vitaly V. Chaban wrote: >> >> Hi, >> &g

[gmx-users] Re: Charge group moving too much

2009-11-10 Thread Vitaly V. Chaban
Hi, The easiest way to fix the problem is to decrease the current timestep. Also please read about -dd and -rdd keys of mdrun utility. Vitaly > One of my simulations just crashed after 660,000 steps with the > following message: > > Step 660560: > The charge group starting at atom 4100 moved th

[gmx-users] Re: Increase temperature to 550K

2009-11-10 Thread Vitaly V. Chaban
I think the problem is hidden in your water force field model. > The simulation system is merely water + one lysozyme. > I increase temperate to 550K. > > Then, the simulation broke. > The following message is shown. > > MX:hpc0011:Remote endpoint is closed, peer=00:60:dd:48:14:5b (hpc0010:0) > m

[gmx-users] Re: Increase temperature to 550K

2009-11-10 Thread Vitaly V. Chaban
> Vitaly V. Chaban wrote: >> I think the problem is hidden in your water force field model. > > On what basis? The user hasn't provided any information to suggest it's > other than a poorly-conditioned integration. On the basis of intuition. Most of the known water

[gmx-users] Re: Force Fields

2009-11-19 Thread Vitaly V. Chaban
The selection of FF depends on what you study. Perhaps, none of them will be useful for you. The existing ones differ at least by the procedures if generation although many parameters are very close. > > What are the main differences between the force fields used in gromacs?  In > which scenarios

[gmx-users] Heat of vaporization in gromacs

2009-12-03 Thread Vitaly V. Chaban
Hi, Please let me know which of the gromacs utilities is responsible for heat of vaporization calculation. Thanks in advance, Vitaly -- 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 b

[gmx-users] 1-4 interaction energies in g_energy

2009-12-03 Thread Vitaly V. Chaban
Hi, Please suggest why is the reason that 1-4 term is not displayed (to be selected for calculation) in the g_energy utility? I have two kinds of objects (with 25 and 5 atoms) in my system - with "nexcl" equal to 3 and 2. I guess the 1-4 term may not be displayed if there is no 1-4 interactions bu

[gmx-users] Re: 1-4 interaction energies in g_energy

2009-12-03 Thread Vitaly V. Chaban
>> Hi, >> >> Please suggest why is the reason that 1-4 term is not displayed (to be >> selected for calculation) in the g_energy utility? I have two kinds of >> objects (with 25 and 5 atoms) in my system - with "nexcl" equal to 3 >> and 2. I guess the 1-4 term may not be displayed if there is no 1-

[gmx-users] intramolecular LJ and Coulomb interactions

2009-12-05 Thread Vitaly V. Chaban
Hi, Is there a way to output only intramolecular interaction (LJ and Coulomb but not bond, angle, etc) energies (with a respect of exclusions, of course). It seems LJ-(SR) and Coulomb-(SR) terms in g_energy include all interaction, both inter- and intramolecular ones. Thanks, Vitaly -- gmx-users

[gmx-users] Re: intramolecular LJ and Coulomb interactions

2009-12-05 Thread Vitaly V. Chaban
On Sat, Dec 5, 2009 at 10:40 PM, Justin A. Lemkul wrote: > > > Vitaly V. Chaban wrote: >> >> Hi, >> >> Is there a way to output only intramolecular interaction (LJ and >> Coulomb but not bond, angle, etc) energies (with a respect of >> exclusions,

  1   2   3   >