[gmx-users] tabulated dihedral potential

Mark Abraham mark.j.abraham at gmail.com
Tue Jan 19 18:50:16 CET 2016


Hi,

Can't tell, but you should test such things by using mdrun -rerun to write
energies and forces for a trivial input that you can verify by hand (e.g.
four atoms with a single interaction)

Mark

On Tue, Jan 19, 2016 at 6:01 PM Parvez Mh <parvezmh89 at gmail.com> wrote:

> (Repost)
> Dear all:
>
> I am trying to use tabulated dihedral potential in my simulation. As a test
> case, i tried with something known case like octane. I know the
> Ryckaert-Bellemans function coefficient of dihedral. So i tried to create a
> table from this known case using following fortran code
>
> program table
>     implicit none
>     double precision::a0,a1,a2,a3,a4,a5
>     double precision::x
>     double precision::fx,dfx
>     character(len=15)::filename
>     integer::angle_deg
>     a0=8.3945
>     a1=16.7787
>     a2=1.1332
>     a3=-26.3070
>     a4=0.000676
>     a5=-0.00001816
>     filename='table_d0.xvg'
>     open(unit=100,file=filename,action='write',status='replace')
>
>     do angle_deg=-180,180,1
>
>         x=angle_deg*0.0174533    !converting degree to radian
>
>         fx=a0-a1*cos(x)+a2*(cos(x))**2-a3*(cos(x))**3+a4*(cos(x))**4 &
>         -a5*(cos(x))**5
>
>         dfx=sin(x)*(a1-2*a2*cos(x)+3*a3*(cos(x))**2-4*a4*(cos(x))**3+ &
>         5*a5*(cos(x))**4)     !analytical differentiation of
> Ryckaert-Bellemans function
>
>         write(100,*)angle_deg,fx,-dfx
>     end do
> end program table
>
> I used force constant equal to 1 in topology file.
> In npt equlibration, my simulation crashed with lot of LINCS error.
> Would somebody help me figure out where i did the mistake.
>
> Regards
> Masrul
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
>


More information about the gromacs.org_gmx-users mailing list