Dear gromacs users, I am trying to simulate a 225aa protein at 300K in water, with OPLSAA force filed, tip4p water model, using the parameters below (for pressure equilibration) during position restrained md.
integrator = md dt = 0.002 nsteps = 50000 nstcomm = 10 nstlist = 5 ns_type = grid pbc = xyz rlist = 1.40 coulombtype = PME rcoulomb = 1.40 vdwtype = Switch rvdw_switch = 0.9 rvdw = 1.00 fourierspacing = 0.12 pme_order = 4 ewald_rtol = 1e-5 optimize_fft = yes Tcoupl = nose-hoover tc_grps = Protein Non-Protein tau_t = 0.4 0.4 ref_t = 300 300 nh-chain-length = 1 gen_vel = no continuation = yes constraints = all-bonds constraint-algorithm = LINCS lincs-order = 4 lincs-iter = 1 lincs-warnangle = 30 when I calculate the densities, the average came out to be 1025.91 kg/m^3, Could anyone please clarify me the reason for this density. Thanking you With regards Kavya
-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists