| 查看: 4771 | 回復(fù): 11 | ||
zdyhbz新蟲 (初入文壇)
|
[求助]
NPT下的壓力 已有4人參與
|
|
大家好!最近碰到個問題想請教一下!打擾了! 我想在NPT系綜下對水弛豫一段時間,壓力給定1個大氣壓。但是不知道為什么壓力跳躍很厲害,體系遠(yuǎn)遠(yuǎn)偏離給定的P。并且還出現(xiàn)負(fù)值。大家能不能幫我看看到底哪里出錯了? 另外怎么會出現(xiàn)負(fù)的壓力呢?真的想不明白。。。 #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 |
版主 (著名寫手)
新蟲 (初入文壇)
新蟲 (初入文壇)

|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|