Dear Users,
First of all thank you Mark, yes you are right. it was problem of force
field and topoloy file mismatch.
I have made changes in the ffG43a1 residue toplogy file for DMSO., but
still getting the same error. Please let me know if i am doing some
mistake.
*******************changes made***********************
[ DMSO ]
[ atoms ]
SD SD 0.13900 0
OD OD -0.45900 0
CD1 CD 0.16000 0
CD2 CD 0.16000 0
[ bonds ]
SD OD gb_38
SD CD1 gb_39
SD CD2 gb_39
OD CD1 gb_44
OD CD2 gb_44
CD1 CD2 gb_45
[ angles ]
; ai aj ak gromos type
[ impropers ]
; ai aj ak al gromos type
[ dihedrals ]
; ai aj ak al gromos type
*********************DMSO before changes**********************
[ DMSO ]
[ atoms ]
SDmso SDMSO 0.13900 0
ODmso ODMSO -0.45900 0
CDms1 CDMSO 0.16000 0
CDms2 CDMSO 0.16000 0
[ bonds ]
SDmso ODmso gb_38
SDmso CDms1 gb_39
SDmso CDms2 gb_39
ODmso CDms1 gb_44
ODmso CDms2 gb_44
CDms1 CDms2 gb_45
[ angles ]
; ai aj ak gromos type
[ impropers ]
; ai aj ak al gromos type
[ dihedrals ]
; ai aj ak al gromos type
***********************************************************
dmso.itp file which i am using for dmso is following.
********************************************************
[ moleculetype ]
; name nrexcl
DMSO 2
[ atoms ]
; nr type resnr residu atom cgnr charge
1 SD 1 DMSO SD 1 0.139
2 CD 1 DMSO CD1 1 0.16
3 OD 1 DMSO OD 1 -0.459
4 CD 1 DMSO CD2 1 0.16
[ bonds ]
; ai aj funct b0 kb
1 2 1 0.195 376560 ; kb from Methionine
1 3 1 0.153 502080 ; kb from C=O bond
1 4 1 0.195 376560
[ angles ]
; ai aj ak funct theta kt
3 1 2 1 106.75 460.240
3 1 4 1 106.75 460.240
2 1 4 1 97.40 460.240 ; kt from Methionine
[ dihedrals ]
; ai aj ak al funct xi kxi
1 2 4 3 2 35.264 334.720 ; default params
Any suggestions, hints or helps are most appreciated.
Thanks in advance,
With Regards
anil
>
> Message: 1
> Date: Sat, 12 May 2007 16:59:40 +0530 (IST)
> From: "Anil Kumar" <[EMAIL PROTECTED]>
> Subject: [gmx-users] Error while running grompp for em (MD) in Urea To:
[email protected]
> Message-ID: <[EMAIL PROTECTED]>
Content-Type: text/plain;charset=iso-8859-1
>
> Dear All,
>
> Can anyone help me in running simulation in explicit urea.
> When i was using urea+h2o.gro and urea.itp and tried to do simulation in
explicit urea.
> After doing genbox, i am not able to do the grompp for energy
minimization
> and getting the error which is given below.
>
> I have added masses of all the atoms of urea in atommass.dat file, but
still getting the same error.
> Please help me, where does am i making the mistakes.
>
> processing coordinates...
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> Warning: atom names in base.top and after_sol.gro don't match (HW2 -
HW3)
> Warning: atom names in base.top and after_sol.gro don't match (OW - OW1)
Warning: atom names in base.top and after_sol.gro don't match (HW1 -
HW2)
> (more than 20 non-matching atom names)
> WARNING 3 [file "base.top", line 808]:
> 6150 non-matching atom names
> atom names from base.top will be used
> atom names from after_sol.gro will be ignored
>
> double-checking input for internal consistency...
> renumbering atomtypes...
> converting bonded parameters...
> # BONDS: 945
> # G96BONDS: 128
> # ANGLES: 1215
> # G96ANGLES: 190
> # PDIHS: 1134
> # IDIHS: 485
> # LJ14: 1242
> # SETTLE: 2050
> ERROR 4 [file "base.top", line 808]:
> atom C1 (Res UREA-11) has mass 0
>
> ERROR 5 [file "base.top", line 808]:
> atom O2 (Res UREA-11) has mass 0
>
> ERROR 6 [file "base.top", line 808]:
> atom N3 (Res UREA-11) has mass 0
>
> ERROR 7 [file "base.top", line 808]:
> atom H4 (Res UREA-11) has mass 0
>
> ERROR 8 [file "base.top", line 808]:
> atom H5 (Res UREA-11) has mass 0
>
> ERROR 9 [file "base.top", line 808]:
> atom N6 (Res UREA-11) has mass 0
>
> ERROR 10 [file "base.top", line 808]:
> atom H7 (Res UREA-11) has mass 0
>
> -------------------------------------------------------
> Program grompp, VERSION 3.3.1
> Source code file: fatal.c, line: 416
>
> Fatal error:
> Too many warnings, grompp terminated
> -------------------------------------------------------
>
> "I'm An Oakman" (Pulp Fiction)
>
> Thanks in advance.
>
> With Warm Regards
> anil
> --
> (¨`.´¨) Always
> `.¸(¨`.´¨) Keep
> (¨`.´¨)¸.´ Smiling!
> `.¸.´
> «´`.(*.¸(`.¸ ¸.´)¸.*).´`»
> «´¨*.¸¸. * ANIL *.¸¸.*¨`»
> «´`.(¸.´(¸.* *.¸)`.¸).´`»
>
> ANIL KUMAR(Research Scholar),
> Bio-Organic Lab No-336(2nd Floor),
> Dept. of Chemistry,I.I.T.Bombay,Powai,
> Mumbai-400076,
> Ph. No.-022-25764780(Lab)
> -----------------------------------------
> Residence:-
> Hostel#1,Room#297,IIT Bombay,Powai,
> Mumbai-400076,Ph.No.:-022-25721017(Hostel)
> -------------------------------------------------
> Web:http://chemanil.googlepages.com/
> --------------------------
> "Education is a progressive discovery of our ignorance"
>
> - Will Durant
> ----------------------------
>
>> Send gmx-users mailing list submissions to
>> [email protected]
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://www.gromacs.org/mailman/listinfo/gmx-users
>> or, via email, send a message with subject or body 'help' to
>> [EMAIL PROTECTED]
>> You can reach the person managing the list at
>> [EMAIL PROTECTED]
>> When replying, please edit your Subject line so it is more specific
than "Re: Contents of gmx-users digest..."
>> Today's Topics:
>> 1. Re: CNT Potassium doping with GROMACS (George Abadir)
>> 2. Re: Cannot find fftw3f (Mark Abraham)
>> 3. Re: suggest small revision to genbox (bharat v. adkar)
>> 4. Re: suggest small revision to genbox (David van der Spoel) 5.
Error while running MD in DMSO! (Anil Kumar)
>> ----------------------------------------------------------------------
Message: 1
>> Date: Fri, 11 May 2007 18:52:43 -0700
>> From: George Abadir <[EMAIL PROTECTED]>
>> Subject: Re: [gmx-users] CNT Potassium doping with GROMACS
>> To: [EMAIL PROTECTED], Discussion list for GROMACS users
>> <[email protected]>
>> Message-ID: <[EMAIL PROTECTED]>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>> Another point: is it possible to use pdb2gmx to do the job instead of
x2top? I tried and used the GROMOS96 43b1 vacuum force field (since the
doping by potassium is in vacuum), but i got the error"Residue 'CNT'
not
>> found in residue topolgy database'. I went to the 'ffg43b1.rtp' file in
the /share/top folder and defined it as '[CNT]
>> [atoms]
>> C C 0 0'
>> but still get the same error.
>> Your help is much appreciated,
>> Thank you in advance,
>> George
>> George Abadir wrote:
>>> Hi,
>>> I want to make a simulation of CNT in an environment of potassium
>>> atoms to see how they will be adsorbed on the CNT surface (this is a
way of doping carbon nanotubes). I built a .pdb file with a CNT and
some potassium atoms with random coordinates around the nanotube. When
I run x2top to have my topology file, I get the error:"No forcefield
type for atom K with 0 bond". The problem is that the atoms of
potassium should indeed be not bonded to anything at the beginning of
the simulation. Does anybody have a way to describe such a system and
avoid this error?
>>> Thank you very much in advance,
>>> Regards,
>>> George
>>> _______________________________________________
>>> gmx-users mailing list [email protected]
>>> 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
>> ------------------------------
>> Message: 2
>> Date: Sat, 12 May 2007 12:22:41 +1000
>> From: Mark Abraham <[EMAIL PROTECTED]>
>> Subject: Re: [gmx-users] Cannot find fftw3f
>> To: Discussion list for GROMACS users <[email protected]>
>> Message-ID: <[EMAIL PROTECTED]>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>> Omovie, Sheyore John wrote:
>>> Hi,
>>> Even after following the instructions in the HowTo Cygwin and gromacs, I
>>> still come up
>>> with the configure error:Cannot find fftw3f.
>> The most likely thing is that you haven't actually followed the
instructions precisely...
>> http://www.gromacs.org/pipermail/gmx-users/2006-April/021147.html In
particular you need "--enable-float" with configure for fftw.
>>> I ahve also tried specifying the header file directory, it still doesn't
>>> work.
>> No, that wouldn't help.
>>> I would appreciate any advice on overcoming this problem.
>> In the absence of more information from you, all I can suggest is
trying
>> again :-) In extremis, re-installing cygwin has helped another user. Mark
>> ------------------------------
>> Message: 3
>> Date: Sat, 12 May 2007 12:05:44 +0530 (IST)
>> From: "bharat v. adkar" <[EMAIL PROTECTED]>
>> Subject: Re: [gmx-users] suggest small revision to genbox
>> To: Discussion list for GROMACS users <[email protected]>
>> Message-ID: <[EMAIL PROTECTED]>
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
>> genbox works fine for me. i use 3.3.1 version. it works fine with and
without script.
>> but i don't understand how will u get rid of backup files in the permanent
>> directory from whatever u have mentioned, because gromacs creates the
output topology in the same location where it reads the input topology
from. so temperory directory thing will work only if u make an link of
top
>> file from permanent directory and use this linked version as input,
then
>> move the new top from temp to permanent directory.
>> bharat
>> On Fri, 11 May 2007, David Mobley wrote:
>>> All,
>>> I'm trying to write some wrappers to set up a project using standard
gromacs tools. Anyway, it looks to me like genbox is not handling its
path names correctly. In particular, here's what I did:
>>> - Make temporary directory in which to run genbox
>>> - Move to temporary directory, run genbox (reading input files from
permanent directory; writing input files to permanent directory) -
Change back to starting directory; delete temporary directory When I
do this, I get a blank output topology file in my permanent directory.
>>> It appears what's going on is that the file temp.top is written into
my temporary directory, but never moved into the permanent directory;
I suspect genbox is not handling the absolute/relative path issue
correctly here?
>>> I can work around, but this is just a bit annoying -- when I script
it, I end up with a bunch of "#temp.top.1#" etc files in my working
directory if I use editconf in the normal way, which is why I was
using the temporary directory thing to begin with.
>>> Thanks,
>>> David
>>> _______________________________________________
>>> gmx-users mailing list [email protected]
>>> 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
>> ------------------------------
>> Message: 4
>> Date: Sat, 12 May 2007 11:13:00 +0200
>> From: David van der Spoel <[EMAIL PROTECTED]>
>> Subject: Re: [gmx-users] suggest small revision to genbox
>> To: Discussion list for GROMACS users <[email protected]>
>> Message-ID: <[EMAIL PROTECTED]>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>> David Mobley wrote:
>>> All,
>>> I'm trying to write some wrappers to set up a project using standard
gromacs tools. Anyway, it looks to me like genbox is not handling its
path names correctly. In particular, here's what I did:
>>> - Make temporary directory in which to run genbox
>>> - Move to temporary directory, run genbox (reading input files from
permanent directory; writing input files to permanent directory) -
Change back to starting directory; delete temporary directory When I
do this, I get a blank output topology file in my permanent directory.
>>> It appears what's going on is that the file temp.top is written into
my temporary directory, but never moved into the permanent directory;
I suspect genbox is not handling the absolute/relative path issue
correctly here?
>>> I can work around, but this is just a bit annoying -- when I script
it, I end up with a bunch of "#temp.top.1#" etc files in my working
directory if I use editconf in the normal way, which is why I was
using the temporary directory thing to begin with.
>> Hi David,
>> plz submit a bugzilla with all input and scripts.
>> --
>> David.
>> ________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
>> Husargatan 3, Box 596, 75124 Uppsala, Sweden
>> phone: 46 18 471 4205 fax: 46 18 511 755
>> [EMAIL PROTECTED] [EMAIL PROTECTED] http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
------------------------------
>> Message: 5
>> Date: Sat, 12 May 2007 14:03:11 +0530 (IST)
>> From: "Anil Kumar" <[EMAIL PROTECTED]>
>> Subject: [gmx-users] Error while running MD in DMSO!
>> To: [EMAIL PROTECTED], [email protected]
>> Message-ID: <[EMAIL PROTECTED]>
Content-Type: text/plain;charset=iso-8859-1
>> Dear All,
>> I am running my simulations in explicit MeOH and DMSO. In case of MeOH, it
>> is running fine but while running in DMSO.
>> After running genbox i am having problem in running grompp for energy
minimization.
>> Error which i got is given below.
>> Cleaning up temporary file gromppkkHvWC
>> Fatal error: Atomtype 'SD' not found!
>> Please, anyone let me know why it is happening when i had used the
dmso.gro and dmso.itp files which is available by default to top
folder.
>> Thanks in advance.
>> With Warm Regards
>> Anil
>> --
>> (¨`.´¨) Always
>> `.¸(¨`.´¨) Keep
>> (¨`.´¨)¸.´ Smiling!
>> `.¸.´
>> «´`.(*.¸(`.¸ ¸.´)¸.*).´`»
>> «´¨*.¸¸. * ANIL *.¸¸.*¨`»
>> «´`.(¸.´(¸.* *.¸)`.¸).´`»
>> ANIL KUMAR(Research Scholar),
>> Bio-Organic Lab No-336(2nd Floor),
>> Dept. of Chemistry,I.I.T.Bombay,Powai,
>> Mumbai-400076,
>> Ph. No.-022-25764780(Lab)
>> -----------------------------------------
>> Residence:-
>> Hostel#1,Room#297,IIT Bombay,Powai,
>> Mumbai-400076,Ph.No.:-022-25721017(Hostel)
>> -------------------------------------------------
>> Web:http://chemanil.googlepages.com/
>> --------------------------
>> "Education is a progressive discovery of our ignorance"
>> - Will Durant
>> ----------------------------
>> ------------------------------
>> _______________________________________________
>> gmx-users mailing list
>> [email protected]
>> http://www.gromacs.org/mailman/listinfo/gmx-users
>> Please search the archive at http://www.gromacs.org/search before posting!
>> End of gmx-users Digest, Vol 37, Issue 44
>> *****************************************
>
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 12 May 2007 19:18:46 +0100 (BST)
> From: jagannath mondal <[EMAIL PROTECTED]>
> Subject: RE: [gmx-users] problem with freeze-group : Large VCM(Berk
> Hess)
> To: Discussion list for GROMACS users <[email protected]>
> Message-ID: <[EMAIL PROTECTED]>
> Content-Type: text/plain; charset=iso-8859-1
>
> Hi Berk,
> Thanks a lot for your reply. According to your
> suggestion, I did not remove the cmm for the WHOLE
> molecule ---- i.e I used 'cmm-mode=None' for the
> entire system. (Did you mean that I should have
> removed the restriction only for freeze-group or for
> whole system?)
> But now I have two problems:
> 1) the simulation is now going for longer time but
> finally again crashes.
> 2) I freezed the mainchain because I did not want
> the mainchain
>
> move at all and I wanted the only side-chain move. But
> , when I am
>
> visualising the trajectory in VMD, I am finding that
> whole molecule is
>
>
> tumbling including the main-chain i.e the whole
> molecule is changing
>
> its position. As ususal, I did not use pressure
> coupling . Do you think that the movement of whole
> molecule is expected during freeze group simulation?
> should I modify something else in my parameter file?
> I am again giving the .mdp file.
> --- Berk Hess <[EMAIL PROTECTED]> wrote:
>
>> >From: jagannath mondal <[EMAIL PROTECTED]>
>> >Reply-To: Discussion list for GROMACS users
>> <[email protected]>
>> >To: [email protected]
>> >Subject: [gmx-users] problem with freeze-group :
>> Large VCM
>> >Date: Fri, 11 May 2007 00:19:53 +0100 (BST)
>> >
>> >Hi Gromacs user,
>> >
>> > I am a gromacs beginner struggling with
>> >freeze-group simulation.
>> > I was trying to simulate a
>> beta-peptide(un-natural
>> >peptide ) by relaxing only the side-chains but I am
>> >getting error regarding large VCM and
>> >The system has only 1 peptide (14-residue ) and No
>> >solvent. So, for this purpose , I generated a
>> >freeze-group which contains all the main-chain
>> atoms.
>> >Initially I was using pressure-coupling and it was
>> >giving error in simulation. Later I found many
>> >discussion on this freeze-group simulation in
>> >user-archive and manual and accordingly
>> > I did not use pressure-coupling in my simulation
>> >and before the simulation, I minimised my peptide
>> >using steep
>> > integrator and then with the minimised structure
>> I
>> >tried a MD run. But after 40 ps, the mdrun crashes
>> >with complaint about nsgrid and large VCM:
>> You should not remove com motion when using freeze
>> groups,
>> since the com is no longer free to move.
>> We should let grompp print a warning for this.
>> Berk.
> _________________________________________________________________
>> Play online games with your friends with Messenger
>> http://www.join.msn.com/messenger/overview
>> _______________________________________________
>> gmx-users mailing list [email protected]
>> 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
>
>
>
> Office firewalls, cyber cafes, college labs, don't allow you to
> download CHAT? Click here:
> http://in.messenger.yahoo.com/webmessengerpromo.php
>
>
> ------------------------------
>
> Message: 3
> Date: Sun, 13 May 2007 09:12:01 +1000
> From: Mark Abraham <[EMAIL PROTECTED]>
> Subject: Re: [gmx-users] Error while running MD in DMSO!
> To: Discussion list for GROMACS users <[email protected]>
> Message-ID: <[EMAIL PROTECTED]>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Anil Kumar wrote:
>> Dear All,
>> I am running my simulations in explicit MeOH and DMSO. In case of MeOH, it
>> is running fine but while running in DMSO.
>> After running genbox i am having problem in running grompp for energy
minimization.
>> Error which i got is given below.
>> Cleaning up temporary file gromppkkHvWC
>> Fatal error: Atomtype 'SD' not found!
>> Please, anyone let me know why it is happening when i had used the
dmso.gro and dmso.itp files which is available by default to top
folder.
>
> They have to match the force field you're using, and the structure file
you're using.
>
> Mark
>
>
> ------------------------------
>
> Message: 4
> Date: Sun, 13 May 2007 09:16:15 +1000
> From: Mark Abraham <[EMAIL PROTECTED]>
> Subject: Re: [gmx-users] Error while running grompp for em (MD) in
> Urea
> To: Discussion list for GROMACS users <[email protected]>
> Message-ID: <[EMAIL PROTECTED]>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Anil Kumar wrote:
>> Dear All,
>> Can anyone help me in running simulation in explicit urea.
>> When i was using urea+h2o.gro and urea.itp and tried to do simulation
in
>> explicit urea.
>> After doing genbox, i am not able to do the grompp for energy
>> minimization
>> and getting the error which is given below.
>> I have added masses of all the atoms of urea in atommass.dat file, but
still getting the same error.
>
> This file isn't used for making topologies.
>
>> Please help me, where does am i making the mistakes.
>> processing coordinates...
>> Warning: atom names in base.top and after_sol.gro don't match (OW -
OW1)
>> Warning: atom names in base.top and after_sol.gro don't match (HW1 - HW2)
>> Warning: atom names in base.top and after_sol.gro don't match (HW2 - HW3)
>
> Well, if these don't match you're doing something wrongly, but I can't
tell what.
>
>> double-checking input for internal consistency...
>> renumbering atomtypes...
>> converting bonded parameters...
>> # BONDS: 945
>> # G96BONDS: 128
>> # ANGLES: 1215
>> # G96ANGLES: 190
>> # PDIHS: 1134
>> # IDIHS: 485
>> # LJ14: 1242
>> # SETTLE: 2050
>> ERROR 4 [file "base.top", line 808]:
>> atom C1 (Res UREA-11) has mass 0
>> ERROR 5 [file "base.top", line 808]:
>> atom O2 (Res UREA-11) has mass 0
>> ERROR 6 [file "base.top", line 808]:
>> atom N3 (Res UREA-11) has mass 0
>> ERROR 7 [file "base.top", line 808]:
>> atom H4 (Res UREA-11) has mass 0
>> ERROR 8 [file "base.top", line 808]:
>> atom H5 (Res UREA-11) has mass 0
>> ERROR 9 [file "base.top", line 808]:
>> atom N6 (Res UREA-11) has mass 0
>> ERROR 10 [file "base.top", line 808]:
>> atom H7 (Res UREA-11) has mass 0
>
> Your force field's .rtp file needs entries for these atoms, probably.
>
> Mark
>
>
> ------------------------------
>
> Message: 5
> Date: Sun, 13 May 2007 09:23:26 +1000
> From: Mark Abraham <[EMAIL PROTECTED]>
> Subject: Re: [gmx-users] Re: Parameterizing a force field (Mark
> Abraham)
> To: Discussion list for GROMACS users <[email protected]>
> Message-ID: <[EMAIL PROTECTED]>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Lee-Ping wrote:
>> Thank you both for the response! I looked at the -j option in mdrun,
but it looks like that has about as much documentation as ffscan (i.e.
none at all). I also fruitlessly checked the mailing list archives for
the past year and a half; there are some questions about ffscan, but
none of which contained instructions on how to run it.
>> The next step is to look inside the gromacs source code, but there are
a
>> lot of different files and I'm not sure as to which pieces of source
code describe how the input files are parsed. If you could tell me
which files to look at, I would appreciate it a lot. Thanks again!
>
> Mostly, files that parse files end in io.c
>
>> Specifically:
>> 1) What piece of code do I look at to see how mdrun parses the wham.gct
file (turned on by the -j switch)?
>
> Don't know. Start with src/kernel/mdrun.c and see what it calls to parse
this file.
>
>> 2) Which file describes how ffscan parses its inputs?
>
> Start with src/kernel/ffscan.c and see what it calls :-)
>
>> Also: If ffscan performs a grid scan, is it really able to fit 100+
parameters in a reasonable period of time (couple of hours)? Would a
numerical optimization procedure for the parameters be a possible
improvement over this?
>
> ffscan doesn't do a fit, per se, it just does a grid scan and reports
suitable values (see man ffscan). It'd have to be heaps faster than your
alternative scripts calling grompp and mdrun to do your series of single
point energy evaluations.
>
> Mark
>
>
> ------------------------------
>
> _______________________________________________
> gmx-users mailing list
> [email protected]
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before
posting!
>
> End of gmx-users Digest, Vol 37, Issue 45
> *****************************************
>
_______________________________________________
gmx-users mailing list [email protected]
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
