Re: Force field design

From: Axel Kohlmeyer (akohlmey_at_gmail.com)
Date: Tue Jul 22 2014 - 11:14:16 CDT

On Tue, Jul 22, 2014 at 11:40 AM, Kenno Vanommeslaeghe
<kvanomme_at_rx.umaryland.edu> wrote:
> Dear Abishek,
>
> In that case, you should try to make your prof read this:
> http://www.charmm.org/ubbthreads/ubbthreads.php?ubb=showflat&Number=32743#Post32743
> He or she basically gave you an impossible task. I don't want to discourage

i would not only call this an impossible task, but also outright
irresponsible behavior of the PI. not only to his student, but the
community in general. i can only encourage students to own up to this
and *demand* to be properly trained. it is not the job of the
community to compensate for this, especially when this is piled on top
of software and parameters that are provided at no cost and contain
the collected knowledge and experience of many, many people. PIs
exhibiting such irresponsible behavior should be named in public so
that a) students are made aware of the risk they may be taking and b)
their publications including simulation data, can be put under
particularly close scrutiny, since they are likely to be faulty and
misinterpreted.

axel.

> you, but it's the truth. Even if after a lot of sweat and tears, you somehow
> succeed to get a result out of NAMD, it will be hard to make sure you
> avoided all the nontrivial pitfalls and your results are relevant.
>
> Cheers,
>
> Kenno.
>
>
>
> On 07/22/2014 06:33 AM, Abhishek TYAGI wrote:
>>
>> Dear Kenno,
>>
>> Thanyou for your detail reply and instructions, actually I am alone doing
>> work on NAMD in my lab, even in my university, my Prof do not have
>> subsequent understanding of NAMD and even I don't have so any interaction to
>> Post doc/RA or other students, therefore, I directly ask you for the
>> problem. I have to do it myself and try to understand the things and do them
>> carefully.
>>
>> I hope in future you understand my situation.
>>
>> regards
>>
>> Abhishek
>>
>> ________________________________________
>> From: owner-namd-l_at_ks.uiuc.edu <owner-namd-l_at_ks.uiuc.edu> on behalf of
>> Kenno Vanommeslaeghe <kvanomme_at_rx.umaryland.edu>
>> Sent: Monday, July 21, 2014 11:47 PM
>>
>> To: namd-l_at_ks.uiuc.edu
>> Subject: Re: namd-l: Force field design
>>
>> Dear Abhishek,
>>
>> This list is full of very knowledgeable people who are very good at
>> resolving people's errors, but even the greatest guru cannot help you if
>> you don't tell them what exactly your problem is.
>>
>> Since I get the impression you're relatively new to this, allow me to link
>> you to CHARMM's "parameter set discussion" forum guide lines:
>> http://www.charmm.org/ubbthreads/ubbthreads.php?ubb=showflat&Number=29297
>> I'm NOT suggesting you ask your question there. The CHARMM forums are for
>> CHARMM problems and the regulars don't like it if you post VMD/NAMD
>> problems. For NAMD, this mailinglist it the correct support channel.
>> Rather, the reason I'm giving you this link is that most of the good
>> practices it outlines are valid in *any* technical community, including
>> this list. In a nutshell, you are advised to search the web, consult the
>> valid resources (tutorials, FAQs...), experiment a bit yourself and
>> consult your advisor/mentor/senior postdoc before asking a question. Then,
>> if you must ask a question, take plenty of time formulating it, thinking
>> of it from the point of view of the expert trying to answer it. Make it
>> concise but include all relevant information.
>>
>> Good luck,
>>
>> Kenno.
>>
>> On 07/19/2014 08:44 AM, Abhishek TYAGI wrote:
>>>
>>> Dear Kenno,
>>>
>>> Thankyou for the suggestion, I had read all the information provided so
>>> far and gone through the links in detail. But still I have a silly question,
>>> right now I have a.mol2 file and I got CGenFF a.str file, and for psfgen If
>>> I use a.str file it shows error, therefore I used CHARM-GUI to get psf and
>>> pdb file, I tried doe quick md on C-GUI it gives some topology and prm file,
>>> I used them for my simulation and got error, and for other files top and prm
>>> I still confused how to proceed.
>>>
>>> regards
>>> Abhi
>>> ________________________________________
>>> From: owner-namd-l_at_ks.uiuc.edu <owner-namd-l_at_ks.uiuc.edu> on behalf of
>>> Kenno Vanommeslaeghe <kvanomme_at_rx.umaryland.edu>
>>> Sent: Thursday, July 17, 2014 1:24 AM
>>> To: namd-l_at_ks.uiuc.edu
>>> Subject: Re: namd-l: Force field design
>>>
>>> On 07/15/2014 12:55 PM, Abhishek TYAGI wrote:
>>>>
>>>> 1. I tried CGenFF and got a *.str file but the data in the file is
>>>> different from my input file i.e. mol2 file, the pdb format is different
>>>> from the *.str file.
>>>
>>>
>>> The str file contains a force field for your molecule, while the mol2 and
>>> pdb formats contain a geometry. How could the data not be different???
>>>
>>> The str file is a concatenated topology and parameter file. NAMD knows to
>>> ignore the topology section and can read the file as it is.
>>> http://mackerell.umaryland.edu/~kenno/cgenff/program.html#namd
>>> I believe psfgen does the same (read the topology and ignore the
>>> parameters) though I'm less sure about that.
>>>
>>>> 2. I tried to get force field from swissparam website and received the
>>>> different file as *.par, *.top file, I tried MD on this field it runs
>>>> successful, but the question is that is this correct approach?
>>>
>>>
>>> It is strongly encouraged to search the archive before posting a
>>> question.
>>> If you had, you would have discovered the use of SwissParam for
>>> condensed-phase calculations is discouraged, and found advise on better
>>> ways to do it:
>>>
>>> http://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l.2013-2014/0846.html
>>>
>>> http://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l.2013-2014/0893.html
>>>
>>> http://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l.2013-2014/1250.html
>>>
>>> Christopher Mayne's advice in this thread is sound.
>>>
>>
>

-- 
Dr. Axel Kohlmeyer  akohlmey_at_gmail.com  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.

This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:22:38 CST