| 查看: 139 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
[求助]用g_energy得到的pressure.xvg文件问题
|
|||
|
在gromacs中，我用g_energy命令得到的压力随时间变化的数据，对于这结果我有两个问题： 1.在pressure.xvg文件中有两列数据，第一列是时间，第二列是压力值，那么这个压力值的单位是什么？bar？ 2. 在md.mdp文件中我设置的ref_p= 1.0，可是由下面这幅图看，压力值的偏差好像比较大，并且在执行g_energy命令时得到下面结果 Energy Average RMSD Fluct. Drift Tot-Drift ------------------------------------------------------------------------------- Pressure (bar) 7.84002 736.433 736.433 -0.000402568 -1.00642 压力值的平均值不是1，而同样的模拟得到的温度的变化的平均值就是我设定的380K，为什么压力的平均值不在1.0 atm？这是什么原因呢？还是我的设置有问题？下面是我的输入文件，请大家指教！谢谢！ ; md in vacuo ; ; Input file ; cpp = /usr/bin/cpp constraints = none integrator = md tinit = 0 dt = 0.001 nsteps = 2000000 init_step = 0 nstcomm = 1 nstxout = 1000 nstvout = 1000 nstfout = 1000 nstcheckpoint = 1000 ns_type = grid pbc = xyz rlist = 1.2 nstlist = 5 morse = no coulombtype = PME rcoulomb = 1.2 ; Spacing for the PME/PPPM FFT grid fourierspacing = 0.12 ; FFT grid size, when a value is 0 fourierspacing will be used fourier_nx = 0 fourier_ny = 0 fourier_nz = 0 pme_order = 4 vdw_type = Shift rvdw = 1.0 rvdw_switch = 0.9 DispCorr = EnerPres Tcoupl = berendsen tc-grps = System tau_t = 0.1 ref_t = 380 Pcoupl = berendsen Pcoupltype = isotropic tau_p = 1.0 compressibility = 4.5e-5 ref_p = 1.0 gen_vel = yes gen_temp = 380 gen_seed = 173529 [ Last edited by zzgyb on 2008-12-31 at 15:45 ] |
回复此楼» 猜你喜欢
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
-
物理化学论文润色/翻译怎么收费?
已经有199人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习 催化剂设计方向
已经有0人回复
-
法国里昂大学招生CSC 资助博士生
已经有0人回复
-
法国里昂大学接收CSC 资助访问学者
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）
已经有1人回复
2楼2008-10-03 17:44:42