<DIV>Dear gmxers,</DIV>
<DIV>I try to perform NPT MD to obtain the density of system. But the value of density does not stablize even after over 40 ns. The content of mdp file can be seen below. How to deal with it? Please give me some hints. Thanks a lot for any replies.</DIV>
<DIV> </DIV>
<DIV>.mdp file: </DIV>
<DIV>title = </DIV>
<DIV>cpp = /lib/cpp<BR>include = -I../top<BR>define = ;{-DFLEXIBLE,-DPOSRES} <BR>integrator = md<BR>dt = 0.001<BR>comm_mode = Linear ;{Linear,Angular,No}<BR>nstcomm = 10<BR>nstcalcenergy = 1<BR>nsttcouple = 1<BR>nstpcouple = 1<BR>comm_grps = <BR>nsteps = 25000000<BR>nstxout = 1000<BR>nstvout = 1000<BR>nstlog = 1000<BR>nstenergy = 1000<BR>nstxtcout = 1000<BR>nstlist = 10<BR>ns_type = grid ;{simple,grid}<BR>rlist = 1.0<BR>coulombtype = pme<BR>fourierspacing = 0.12 ; [0.12 nm]<BR>optimize_fft = yes ;{no,yes}<BR>rcoulomb = 1.0<BR>rvdw = 1.0<BR>DispCorr = EnerPres ;{no,EnerPres,Ener}<BR>pbc = xyz ; Periodic Boundary Conditions {xyz,no,xy}<BR>tcoupl = berendsen ;no,berendsen,nose-hoover,v-rescale<BR>tc-grps = Other <BR>tau_t = 0.1 <BR>ref_t = 600 <BR>Pcoupl = berendsen ;{no,berendsen,Parrinello-Rahman}<BR>Pcoupltype = isotropic ;{isotropic,semiisotropic;anisotropic;surface-tension}<BR>tau_p = 0.1<BR>compressibility = 4.5e-5<BR>ref_p = 1.0 ;the unit is bar(1 bar =10^5 Pa)<BR>gen_vel = no <BR>gen_temp = 600<BR>gen_seed = 173529<BR>constraints = none ;{none,hbonds,all-bonds,h-angles,all-angles}<BR>constraint_algorithm = LINCS ;{LINCS,SHAKE}<BR>periodic_molecules = no ;{no,yes}</DIV>
<DIV>annealing = no ;{no,single,periodic}<BR>annealing_npoints = 11<BR>annealing_time = 0 100 200 300 400 500 600 700 800 900 1000<BR>annealing_temp = 600 580 560 540 520 500 520 540 560 580 600<BR></DIV>
<DIV>Best regards,</DIV>
<DIV>Dr. Chaofu Wu</DIV>
<DIV><sign signid="0">
<DIV style="FONT-FAMILY: Arial Narrow; COLOR: #909090; FONT-SIZE: 12px">------------------</DIV>
<DIV style="FONT-FAMILY: Verdana; COLOR: #000; FONT-SIZE: 14px">
<DIV><A name=OLE_LINK1><I style="mso-bidi-font-style: normal"><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US>Department of Chemistry and </SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US>Materials Science</SPAN></I></A><I style="mso-bidi-font-style: normal"><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US></SPAN></I></DIV>
<DIV><I style="mso-bidi-font-style: normal"><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US>Hunan </SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US>University of </SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US>Humanities, Science and Technology</SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US>, </SPAN></I></DIV>
<DIV><I style="mso-bidi-font-style: normal"><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US></SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US>Loudi 417000</SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US>, </SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US>the </SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA" lang=EN-US>People's Republic of China</SPAN><SPAN style="LINE-HEIGHT: 200%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt; mso-fareast-font-family: 宋体; mso-ansi-language: EN-US; mso-fareast-language: ZH-CN; mso-bidi-language: AR-SA" lang=EN-US> (P.R. China)</SPAN></I></DIV></DIV></sign></DIV>
<DIV> </DIV>