24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者NgWYong将赠送您 20 个金币

NgWYong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 101.6
散金: 27
帖子: 42
在线: 10.6小时
虫号: 29819190
注册: 2022-06-26
专业: 零件成形制造

[求助] lammps原子扩散过程的NPT控压问题

各位老师好：
我模拟异种原子之间的扩散，首先在1K下用NPT弛豫（fix  2 all npt temp 1 1 0.2 iso 1 1 2 drag 0.2），温度压力控制正常
log：Step       Temp       TotEng       Press          Lz
   9200 1.0285191    -18942.729    0.8686156    90.000694
   9300 1.0278229    -18942.752    1.6054926    90.010722
   9400 1.0444687    -18942.773    1.8562374    90.023797
   9500 1.0280457    -18942.799    1.2425338    90.041611
   9600 1.012377    -18942.823    0.4046528    90.059331
   9700 0.98485296 -18942.842    0.0093609708 90.072429
   9800 0.97200394 -18942.849    0.37225422    90.082723
   9900 0.96710447 -18942.834    1.1952155    90.09377
   10000 0.94556821 -18942.815    1.8227487    90.105221
但后面用NPT将温度升高之后（fix 2 all npt temp 1 1123 0.2 iso 1 1 5 drag 0.2），压力变化毫无规律，这正常吗？
log：Step       Temp       TotEng       Press          Lz
   15000 1082.9175    -17458.712    -1610.4027    92.623064
   15100 1104.0539    -17449.548    1225.136    92.620972
   15200 1116.4221    -17437.688    4216.4516    92.622847
   15300 1113.3924    -17437.27    4253.7854    92.639257
   15400 1101.657    -17456.56    2379.0338    92.670386
   15500 1107.5108    -17463.981    564.36355    92.709329
   15600 1121.5267    -17451.407    -345.28204    92.750105
   15700 1119.7161    -17437.794    -2084.121    92.788302
   15800 1132.8069    -17430.172    -4655.8022    92.81755
   15900 1099.6295    -17426.075    -4663.1892    92.829421
   16000 1124.084    -17414.53    -6356.8777    92.819158
温度升高后用NPT系综模拟扩散时（fix 2 all npt temp 1123 1123 0.2 iso 1 1 10 drag 0.2），压力还是变化很大
Step       Temp       TotEng       Press          Lz
299500 1140.9535    -17413.696    1493.379    92.250118
299600 1112.8234    -17426.355    1146.7604    92.243936
299700 1121.5769    -17442.101    -1090.2118    92.240042
299800 1111.9851    -17442.116    -113.55944    92.233639
299900 1122.0422    -17436.086    141.25512    92.226788
300000 1138.9037    -17422.323    -659.76276    92.22097
1.请问下这是正常的压力变化吗？还是写的不对呢？
2.我先用NVT升温再用NPT控压也没啥效果
本人刚接触不久，请各位老师指点一下，谢谢

回复此楼