[gmx-users] solvation free energy
Daniele Veclani
danieleveclani at gmail.com
Fri Jul 17 09:57:49 CEST 2015
Thank you for your answer.
I read the manual. As I understand it, I should be deleted, from .top file, the
section of non-bonded parameters. It's right?
Of course, after the VdW I'll change the electrostatic transformation.
Best regards
D.V.
2015-07-16 16:55 GMT+02:00 Hannes Loeffler <Hannes.Loeffler at stfc.ac.uk>:
> The couple-* parameters take already care of including the non-bonded
> terms internal to your molecule to correctly describe the transfer of M
> to vacuum. That's the point of those parameters so that you would not
> have to run an additional "correction" in vacuo. See the discussion in
> the manual (it's section 5.3 for Gromacs 4.6.x). You would also need to
> compute the electrostatic transformation.
>
> On Thu, 16 Jul 2015 16:26:44 +0200
> Daniele Veclani <danieleveclani at gmail.com> wrote:
>
> > I'm trying to calculate the solvation free energy of a molecule (M).
> > I have done:
> > M+water ---> dum+water
> >
> > Now I have to do:
> > M(vacuo) --> dum(vacuo)
> >
> > In this case I have a problem, in fact I find a DG = 0.00000 and
> > within the .xvg file there are only zeros.
> >
> > Where is the problem?
> >
> > I use gromacs 5.0.4 and my .mdp file (for lambda 00) is:
> > integrator = sd
> > tinit = 0
> > dt = 0.002
> > nsteps = 2500000
> > comm_mode = angular
> > nstcomm = 100
> > nstxout = 500
> > nstvout = 500
> > nstfout = 0
> > nstlog = 500
> > nstenergy = 500
> > nstxout-compressed = 0
> > cutoff-scheme = group
> > nstlist = 0
> > ns_type = simple
> > pbc = no
> > rlist = 0
> > coulombtype = cutoff
> > rcoulomb = 0
> > epsilon_r = 1
> > vdwtype = cutoff
> > rvdw = 0
> > DispCorr = no
> > fourierspacing = 0.12
> > pme_order = 6
> > ewald_rtol = 1e-06
> > epsilon_surface = 0
> > tcoupl = berendsen
> > tc_grps = system
> > tau_t = 0.1
> > ref_t = 300
> > Pcoupl = no
> > tau_p = 1.0
> > compressibility = 4.5e-05
> > ref_p = 1.0
> > free_energy = yes
> > init_lambda_state = 0
> > delta_lambda = 0
> > calc_lambda_neighbors = 1
> > vdw_lambdas = 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70
> > 0.80 0.90 1.00
> > coul_lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
> > 0.00 0.00 0.00
> > bonded_lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
> > 0.00 0.00 0.00
> > restraint_lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
> > 0.00 0.00 0.00
> > mass_lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
> > 0.00 0.00 0.00
> > temperature_lambdas = 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
> > 0.00 0.00 0.00
> > sc-alpha = 0.5
> > sc-coul = no
> > sc-power = 1.0
> > sc-sigma = 0.3
> > couple-moltype = QUI
> > couple-lambda0 = vdw
> >
> > couple-lambda1 = none
> > couple-intramol = no
> > nstdhdl = 10
> > gen_vel = no
> > constraints = all-bonds ; we only have C-H bonds here
> > constraint-algorithm = lincs
> > continuation = no
> > lincs-order = 12
> >
> > Best regards
> > Daniele Veclani.
>
> --
> 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