[gmx-users] free energy caculation error
Szilárd Páll
pall.szilard at gmail.com
Wed Apr 9 13:52:55 CEST 2014
lmgtfy: http://bit.ly/1lJYVGQ ;)
--
Szilárd
On Wed, Apr 9, 2014 at 12:41 PM, jane <348363213 at qq.com> wrote:
> Hi all,
> I want to caculate the free energy of a molecule, but when I run the MD, it
> appears the error: "There is no domain decomposition for 8 nodes that is
> compatible with the given box and a minimum cell size of 4.44444 nm. Change
> the number of nodes or mdrun option -rdd or -dds".However, even though i
> added the -rdd and -dds, the error is same as the previous. my mdp file as
> follows.
>
> ; we'll use the sd integrator
> ;run parameters
> integrator = sd
> dt = 0.002 ; ps
> nsteps = 5000 ; total 1ns
> nstcomm = 100
>
> ; turn off trajectory writing
> nstxout = 0
> nstvout = 0
> ; cut-offs at 1nm
> rlist = 1.0
> nstlist = 10
> coulombtype = pme
> vdw-type = cut-off
> rvdw = 1.0
> constraints = all-bonds
> ; set temperature to 300K
> tcoupl = v-rescale
> tc-grps = system
> tau-t = 0.2
> ref-t = 300
> ; and pressure to 1 bar
> pcoupl = parrinello-rahman
> ref-p = 1
> compressibility = 4.5e-5
> tau-p = 5
> ; and set the free energy parameters
> free-energy = yes
> couple-moltype = RNA_chain_A
> init-lambda = 0
> ; these 'soft-core' parameters make sure we never get overlapping
> ; charges as lambda goes to 0
> sc-power = 1
> sc-sigma = 0.3
> sc-alpha = 1.0
> ; we still want the molecule to interact with itself at lambda=0
> couple-intramol = no
> couple-lambda1 = vdwq
> couple-lambda0 = none
> foreign-lambda = 0 0.2 0.4 0.6 0.8 0.9 1
> nstdhdl = 10
> separate-dhdl-file = yes
> dhdl-derivatives = yes
> dh_hist_size = 5
> dh_hist_spacing = 50
>
> thanks for advances.
>
>
> --
> View this message in context: http://gromacs.5086.x6.nabble.com/free-energy-caculation-error-tp5015758.html
> Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
> --
> 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