24小时热门版块排行榜

返回列表

gtssongchi

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 4373.1
散金: 136
红花: 2
帖子: 157
在线: 99.5小时
虫号: 1137504
注册: 2010-11-02
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 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是什么含义呢，它们之间是什么关系，负的压强是什么含义呢？

回复此楼