Re: [AMBER] problem with tleap

From: FyD <fyd.q4md-forcefieldtools.org>
Date: Fri, 21 Jun 2013 09:11:17 +0200

Dear Rejwan,

> Thanks for your input. The G09 input file I used to extract RESP
> via Antechamber reads as following. Should I make changes and re-run
> Gaussian in the following line?
>
> #HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt iop(6/50=1)

I think you should read what these IOP means; in particular iop(6/42=6).
iop(6/33=2) is not documented but means 'print the MEP'.
See http://www.gaussian.com/g_tech/g_iops/ov6.htm

iop(6/42=6) is not the default algorithm defined by Singh & Kollman.
See: http://q4md-forcefieldtools.org/RED/
& in particular: http://q4md-forcefieldtools.org/RED/popup/net.php

regards, Francois


> --Link1--
> %chk=molecule
> #HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt iop(6/50=1)
>
> remark line goes here
>
> 0 1
> C 6.6215000000 5.1461000000 0.1163000000
> C 5.7259000000 4.2389000000 -0.3176000000
> C 6.2373000000 6.3368000000 0.9704000000
> --------------------------------
> H 11.0017000000 6.8296000000 2.6926000000
> H 10.5984000000 6.6920000000 -1.6004000000
> H 11.7998000000 9.1604000000 2.5323000000
> H 11.3920000000 9.0243000000 -1.7406000000
> H 11.9987000000 10.2632000000 0.3195000000
>
> g09.gesp
>
> g09.gesp
>
>
>
>
>
>
> --- On Tue, 6/18/13, Rejwan wrote:
>
> &gt; From: Rejwan
> &gt; Subject: Re: [AMBER] problem with tleap
> &gt; To: &quot;AMBER Mailing List&quot;
> &gt; Cc: jason.swails.gmail.com
> &gt; Date: Tuesday, June 18, 2013, 9:46 PM
> &gt; --- On Tue, 6/18/13, Jason
> &gt; Swails  wrote:
> &gt;
> &gt; &gt; From: Jason Swails
> &gt; &gt; Subject: Re: [AMBER] problem with tleap
> &gt; &gt; To: &quot;AMBER Mailing List&quot;
> &gt; &gt; Date: Tuesday, June 18, 2013, 5:09 PM
> &gt; &gt; On Tue, Jun 18, 2013 at 4:06 PM,
> &gt; &gt; Rejwan
> &gt; &gt; wrote:
> &gt; &gt;
> &gt; &gt; &gt; I&#39;m trying to make the topology file
> &gt; of a 
> &gt; &gt; ligand-protein system.
> &gt; &gt; &gt;
> &gt; &gt; &gt; In tLeap suit of Amber I fed both .prepin
> &gt; and .frcmod
> &gt; &gt; of the ligand,  I
> &gt; &gt; &gt; used RESP charge obtained via Gaussian 09
> &gt; Geometry
> &gt; &gt; Optimization for PREPIN
> &gt; &gt; &gt; file. The .frcmod file was created by
> &gt; parmchk routine.
> &gt; &gt; &gt;
> &gt; &gt;
> &gt; &gt; Did you look at the resulting frcmod file?  Did
> &gt; you see
> &gt; &gt; any statements like
> &gt; &gt; [ATTN: Needs Revision]?
> &gt;
> &gt; I LOOKED AT FRCMOD FILE. IT DOES NOT HAVE ANY WARNING AT
> &gt; ALL.
> &gt;
> &gt; &gt;
> &gt; &gt;
> &gt; &gt; &gt; The tLeap error showed only 2 out of 10
> &gt; torsion terms
> &gt; &gt; could be
> &gt; &gt; &gt; implemented. Could you please help how can
> &gt; I extract
> &gt; &gt; the missing Angle and
> &gt; &gt; &gt; Torsion terms ? Also where should I plug
> &gt; the terms (Is
> &gt; &gt; it in  the .frmod
> &gt; &gt; &gt; file ?) to rerun Leap ?  I&#39;m copying
> &gt; the error
> &gt; &gt; screen.
> &gt; &gt; &gt;
> &gt; &gt;
> &gt; &gt; You should put any missing parameters into the
> &gt; frcmod
> &gt; &gt; file.  However,
> &gt; &gt; parmchk should have done this for you (unless there
> &gt; were
> &gt; &gt; bonds between your
> &gt; &gt; ligand and protein, for instance, but this seems
> &gt; unlikely
> &gt; &gt; given that all
> &gt; &gt; missing parameters involve gaff atom types).  Even
> &gt; if
> &gt; &gt; parmchk couldn&#39;t find
> &gt; &gt; any parameters, it will put zeroes in for
> &gt; everything and
> &gt; &gt; indicate that it
> &gt; &gt; needs revision. (You should always check the frcmod
> &gt; file
> &gt; &gt; after parmchk
> &gt; &gt; creates it).  Therefore, the fact that parameters
> &gt; are
> &gt; &gt; still missing
> &gt; &gt; indicates that you&#39;ve probably made some kind
> &gt; of error (did
> &gt; &gt; you remember
> &gt; &gt; to load the frcmod into your tleap session via the
> &gt; [
> &gt; &gt; loadAmberParams ]
> &gt; &gt; command?)
> &gt;
> &gt; PARMCHK DID NOT HAVE ANY COMPLAIN AND OF COURSE I DID IMPORT
> &gt; FRCMOD INTO TLEAP VIA loadAmberParams
> &gt; &gt;
> &gt; __________________
> &gt;
> &gt;
> &gt; IF YOU WANT TO HAVE A LOOK AT THOSE FILES I CAN SEND THEM.
> &gt; HOWEVER, WHAT ELSE COULD BE THE PITFALL ? INTERESTINGLY I
> &gt; RUN LIGEAND-PROTEIN  MD  BEFORE BY SAME
> &gt; ANTECHAMBER PROTOCOL WITHOUT PROBLEM. THANKS.
> &gt;
> &gt; -rejwan
> &gt;
> &gt;
> &gt; &gt;
> &gt; &gt; --
> &gt; &gt; Jason M. Swails
> &gt; &gt; Quantum Theory Project,
> &gt; &gt; University of Florida
> &gt; &gt; Ph.D. Candidate
> &gt; &gt; 352-392-4032
> &gt; &gt; _______________________________________________
> &gt; &gt; AMBER mailing list
> &gt; &gt; AMBER.ambermd.org
> &gt; &gt; http://lists.ambermd.org/mailman/listinfo/amber
> &gt; &gt;
> &gt;



_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jun 21 2013 - 00:30:03 PDT
Custom Search