Hi, Dear users,I have created a box of polymer with about 50 molecules. I'm performing NPT runs to reach to the pressure of 1 bar. After a long time of simulation (15 ns) the average pressure is still about -0.5 and the density has saturated to around 700 kg/m^3 which should become 1200 kg/m^3 at the pressure of one bar. These are the related parameters in .mdp file: ; Pressure coupling is onpcoupl = berendsen ; Pressure coupling on in NPTpcoupltype = isotropic ; uniform scaling of box vectorstau_p = 0.2 ; time constant, in psref_p = 1.0 ; reference pressure, in barcompressibility = 4.5e-5 ; isothermal compressibility of water, bar^-1 Does anyone know how can I fix this? Which parameters should I change to reach the desired pressure and density?Thanks in advance.-Nima Azar -- 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-requ...@gromacs.org.
