| 查看: 1412 | 回复: 3 | ||
gtssongchi木虫 (小有名气)
|
[求助]
fix npt 压强设置为0，但结果却不为零，是正常的吗？
|
|
大家好，我在做一个lammps练习的时候，采用了NPT系综，设置fix 1 all npt temp 300 300 0.1 iso 0 0 1 drag 1，输出的温度在300左右，但是压强却差得比较远，这是为什么呢，是Pdamp设置的太大了吗，手册上说这个值设为1000*timestep比较好，所以我这里设置为1。 Step Lx Ly Lz Press Pxx Pyy Pzz PotEng Temp 0 40.5 40.5 40.5 2496.1233 2446.9902 2534.6541 2506.7256 -13440 300 1000 40.634006 40.634006 40.634006 1119.3254 1175.1523 1093.0898 1089.7341 -13289.697 323.77631 2000 40.639934 40.639934 40.639934 427.59205 315.46653 610.36608 356.94355 -13288.932 303.86628 3000 40.639834 40.639834 40.639834 620.40903 665.98372 721.1629 474.08047 -13286.379 303.00734 4000 40.642068 40.642068 40.642068 521.48289 516.80643 572.54414 475.09809 -13286.259 307.38948 5000 40.633951 40.633951 40.633951 831.76741 819.89919 800.00968 875.39336 -13287.407 299.53837 6000 40.630463 40.630463 40.630463 1334.2687 1495.4409 1448.6485 1058.7169 -13283.273 303.65216 7000 40.623885 40.623885 40.623885 1432.9403 1533.4104 1290.1301 1475.2802 -13286.511 295.65281 8000 40.639712 40.639712 40.639712 634.52083 547.05804 743.16154 613.34291 -13285.981 301.84658 9000 40.644657 40.644657 40.644657 293.13866 349.41597 359.53622 170.46378 -13286.033 295.47994 10000 40.647726 40.647726 40.647726 173.46971 80.18083 386.83148 53.396831 -13284.09 300.44559 11000 40.638741 40.638741 40.638741 584.10248 712.26887 479.18644 560.85214 -13285.847 296.85101 12000 40.630857 40.630857 40.630857 1068.3334 1148.9693 1043.5886 1012.4424 -13285.143 297.78126 13000 40.635333 40.635333 40.635333 1047.9483 1115.2864 1155.3644 873.19403 -13281.729 298.43361 14000 40.641403 40.641403 40.641403 456.52097 413.91916 596.25831 359.38545 -13287.485 300.51631 15000 40.640346 40.640346 40.640346 571.15573 579.90767 485.6046 647.95493 -13285.752 296.11894 16000 40.643871 40.643871 40.643871 424.02343 224.95595 411.42919 635.68515 -13284.767 303.06084 17000 40.641585 40.641585 40.641585 455.81834 589.43501 334.88122 443.1388 -13286.495 299.85166 18000 40.632633 40.632633 40.632633 997.39983 916.15558 987.70674 1088.3372 -13285.251 305.99268 19000 40.621975 40.621975 40.621975 1738.1364 1525.3629 1874.4026 1814.6438 -13284.087 296.45112 20000 40.630799 40.630799 40.630799 962.72247 1026.9758 864.28067 996.91093 -13289.452 301.50417 还想请问一下，这里的Press 、Pxx、 Pyy、 Pzz是什么含义呢，它们之间是什么关系，负的压强是什么含义呢？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
re_searcher
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 329.8
- 帖子: 42
- 在线: 21.9小时
- 虫号: 5530986
- 注册: 2017-02-02
- 专业: 凝聚态物性 II ：电子结构
2楼2017-09-26 09:25:31
3楼2017-09-26 19:09:42
gtssongchi
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4373.1
- 散金: 136
- 红花: 2
- 帖子: 157
- 在线: 99.5小时
- 虫号: 1137504
- 注册: 2010-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2017-09-27 07:55:20
10