求助含1 M NaCl 的水溶液里Na+和Cl-的RDF计算的峰高问题

noir

我模拟了一个含有 1 M NaCl 的水溶液在 298 K，1 atm 下的体系，经过 min，nvt，npt 后，跑了 1 ns。计算此平衡状态下密度为 1026 kg/m^3，做 Na-Ow 的 RDF，均与文献一致，但是做 Na-Cl 的 RDF，峰位置与文献接近，但是峰高却远高于文献值。
想问一下这是怎么回事，或者是计算 RDF 的方法有问题？跑平衡的命令为

1. gmx grompp -f min.mdp -o min -pp min -po min -c solution.gro
   
2. gmx mdrun -v -deffnm min
   
3. gmx grompp -f nvt.mdp -o nvt -pp nvt -po nvt -c min -t min
   
4. gmx mdrun -deffnm nvt
   
5. gmx grompp -f npt.mdp -o npt -pp npt -po npt -c nvt -t nvt
   
6. gmx mdrun -deffnm npt
   
7. gmx grompp -f prd.mdp -o prd -pp prd -po prd -c npt -t npt
   
8. gmx mdrun -deffnm prd

复制代码

计算 RDF 的命令为

1. gmx make_ndx -f prd.gro -o solution.ndx
   
2. gmx rdf -f prd.xtc -s prd.tpr -o gr-Na-Cl.xvg -cn Nr-Na-Cl.xvg -n solution.ndx
   
3. gmx rdf -f prd.xtc -s prd.tpr -o gr-Na-wat.xvg -cn Nr-Na-wat.xvg -n solution.ndx

复制代码

prd.tpr (53.84 KB, 下载次数 Times of downloads: 0)

2019-9-26 10:48 上传 Uploaded

点击下载
Click to download

prd.xtc (1.28 MB, 下载次数 Times of downloads: 1)

2019-9-26 10:48 上传 Uploaded

点击下载
Click to download

prd.gro (103.7 KB, 下载次数 Times of downloads: 1)

2019-9-26 10:48 上传 Uploaded

点击下载
Click to download

prd.mdp (836 Bytes, 下载次数 Times of downloads: 3)

2019-9-26 10:46 上传 Uploaded

点击下载
Click to download

nvt.mdp (640 Bytes, 下载次数 Times of downloads: 3)

2019-9-26 10:45 上传 Uploaded

点击下载
Click to download

npt.mdp (619 Bytes, 下载次数 Times of downloads: 3)

2019-9-26 10:45 上传 Uploaded

点击下载
Click to download

min.mdp (440 Bytes, 下载次数 Times of downloads: 2)

2019-9-26 10:45 上传 Uploaded

点击下载
Click to download

topol.top (1.59 KB, 下载次数 Times of downloads: 8)

2019-9-26 10:45 上传 Uploaded

点击下载
Click to download

solution.gro (67.65 KB, 下载次数 Times of downloads: 2)

2019-9-26 10:45 上传 Uploaded

点击下载
Click to download

sobereva

用的力场，水模型、计算非键作用的方法和设定等等都会明显影响这种问题。光看你的文件看不出什么