Don't forget to zero out the torsional term in the force field. The 1-4
interactions will account for some of the profile, so the torsion term is
just a correction for the 1-4 inadequacies. Then you can do the same scan
with amber as you did with Gaussian and fit the difference to the fourier
terms that define the torsion.
Good luck!
Jason
On Thu, Nov 26, 2009 at 10:09 AM, case <case.biomaps.rutgers.edu> wrote:
> On Thu, Nov 26, 2009, Nahoum Anthony wrote:
> >
> > I want to parameterize a ligand for which the default torsion term given
> > in the .frcmod file by antechamber is inadequate. I've used Gaussian to
> > do a torsion scan and get an energy plot for the full rotation and I
> > want to fit that energy plot using Amber's force field equation. I've
> > read several papers where people have done this sort of things, but I
> > can never quite get if they fit the Gaussian energy with the full force
> > field equation, allowing only the torsional parameters to vary or if
> > they try to fit only the torsional part of the force field equation to
> > the Gaussian results. What is the correct procedure, if any ?
>
> The former: you need to compare the total energy from Gaussian to the total
> energy from the force field.
>
> ...good luck...dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
---------------------------------------
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Graduate Student
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Nov 26 2009 - 08:00:02 PST
