| 查看: 5008 | 回复: 11 | ||
[求助]
NPT下的压力 已有4人参与
|
|
大家好！最近碰到个问题想请教一下!打扰了！ 我想在NPT系综下对水弛豫一段时间，压力给定1个大气压。但是不知道为什么压力跳跃很厉害，体系远远偏离给定的P。并且还出现负值。大家能不能帮我看看到底哪里出错了？ 另外怎么会出现负的压力呢？真的想不明白。。。 #in文件 units metal atom_style full bond_style harmonic angle_style harmonic processors 3 2 2 read_data MSD_293.lammps_data group Water type 1 2 group Oxygen type 1 group hydrogen type 2 pair_style lj/cut/coul/long 15.0 kspace_style pppm 10e-4 pair_coeff 1 1 0.006738853 3.166 pair_coeff 1 2 0.0 0.0 pair_coeff 2 2 0.0 0.0 velocity all create 293.0 87287 dist gaussian rot yes units box neigh_modify every 5 delay 10 check yes page 1000000 one 10000 fix 1 all shake 0.0001 10 0 b 1 a 1 fix 2 all npt temp 293. 293. 0.1 iso 1.01325 1.01325 1.0 fix 3 all efield 0.0 0.0 0.0 thermo_style custom step pe ke etotal pxx pyy pzz lx ly lz temp press vol thermo 1000 timestep 0.001 run 250000 #log文件 Step Pxx Pyy Pzz Press 0 181188.5 183019.4 203592.6 189266.8 1000 -1751.91 -1655.94 -1439.9 -1615.92 2000 -595.486 -744.252 -360.922 -566.887 3000 -484.424 -187.144 -554.685 -408.751 4000 -78.0598 -66.1142 87.11058 -19.0211 5000 121.5207 -145.694 -50.427 -24.8668 6000 217.5663 47.1163 161.9961 142.2262 7000 -214.931 -35.5806 -99.7462 -116.753 8000 -14.925 -95.959 -314.102 -141.662 9000 146.2439 -91.2819 151.249 68.73699 10000 -163.258 -484.921 -472.64 -373.606 11000 71.24688 -233.925 27.50841 -45.0566 12000 270.4836 211.9671 81.42193 187.9576 13000 -54.2666 -10.4728 18.88591 -15.2845 14000 -162.172 -69.7858 -71.7508 -101.236 15000 210.7572 34.77109 202.5613 149.3632 16000 -15.7762 -42.257 -190.674 -82.9024 17000 -78.0043 423.1222 -111.782 77.77874 18000 -4.91048 -282.987 -139.973 -142.624 19000 -176.231 206.6059 9.343061 13.23946 20000 -76.6895 272.7583 252.8252 149.6313 21000 76.41657 166.0637 -156.414 28.68863 22000 -567.576 -99.6713 45.87968 -207.123 23000 -82.8798 152.1289 -195.899 -42.2166 24000 -491.898 180.091 -365.879 -225.895 25000 183.6895 128.4664 261.8275 191.3278 26000 362.8234 262.1277 81.00841 235.3199 27000 -97.8331 -222.905 -14.8982 -111.879 28000 14.95781 -8.2984 2.938879 3.199427 29000 102.8558 120.8838 285.0376 169.5924 30000 103.8935 -96.6333 94.75559 34.00526 31000 10.09829 -305.575 -158.684 -151.387 32000 0.789599 -113.531 -23.0684 -45.2698 33000 107.5283 49.42627 286.5654 147.84 34000 296.2923 282.9671 200.8988 260.0527 |
回复此楼» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有22人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用Lammps模拟计算水的密度
已经有10人回复
- lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′）
已经有14人回复
- gromacs npt出现问题
已经有3人回复
- Gromacs模拟温度耦合和压力耦合选择算法的问题，感谢！！！
已经有4人回复
- lammps中的npt系综下逐步降温问题
已经有3人回复
- lammps中界面的建立
已经有9人回复
- 分子模拟中MSD的问题讨论
已经有25人回复
- 动力学NPT压力
已经有25人回复
- nvt下模拟之后获得压力的问题
已经有5人回复
- MS模拟气体扩散时的参数设置问题
已经有4人回复
- NPT系综下判断是否平衡的标志
已经有12人回复
- 关于NVT系统下的压力设定
已经有10人回复
- 第一原理动力学的系综选择，和各种系综的输入文件Example
已经有13人回复
- Gibbs系综和NPT系综的区别
已经有4人回复
- 【讨论】npt压力控制问题
已经有10人回复
- 【求助】NPT分子动力学程序
已经有5人回复
- 【求助】用NAMD 做NPT模拟遇到下面的问题
已经有10人回复
- 【讨论】GCMC系综中化学势如何确定
已经有10人回复
- 【求助】Discover模块下的压力单位GPa和平常的压力单位一样吗？
已经有3人回复
abinitio
版主 (著名写手)
- 模拟EPI: 2
- 应助: 221 (大学生)
- 贵宾: 0.093
- 金币: 7365.2
- 散金: 3917
- 红花: 57
- 帖子: 1317
- 在线: 485.2小时
- 虫号: 368545
- 注册: 2007-05-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2015-04-29 07:41:19
3楼2015-04-29 10:31:56
克里斯蒂亚诺
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 204.1
- 散金: 17
- 帖子: 16
- 在线: 41.7小时
- 虫号: 2967274
- 注册: 2014-02-13
- 性别: GG
- 专业: 传热传质学
4楼2015-06-23 10:06:57
5楼2015-06-24 09:48:02
|
本帖内容被屏蔽 |
6楼2015-06-30 16:04:44
7楼2015-07-01 09:51:10
8楼2015-07-01 14:52:05
9楼2015-07-01 14:55:00
10楼2015-07-01 17:30:15
20