| 查看: 178 | 回复: 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 ] |
回复此楼» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有139人回复
-
祈福青基
已经有34人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
2楼2008-10-03 17:44:42