| 10 | 1/1 | 返回列表 |
| 查看: 2518 | 回復(fù): 9 | ||
fengshaxia21鐵蟲 (初入文壇)
|
[求助]
SPC水分子模型,求密度有些許問題 已有1人參與
|
|
RT,我用MS建立了一個6000分子的SPC水分子模型。但是得到的在293K的密度總是0.96這樣,距離實際值或者SPC水分子的實驗值都有0.2~0.3的差距(雖然看上去不算大,不過這沒關(guān)系嗎?) 我想知道這是什么導(dǎo)致的,density=mass(all)/(參數(shù)*vol),但是mass是不變的,密度的體積就只跟體積有關(guān)了。為什么在溫度到達293附近的時候波動還比較大(270~300上下波動),密度的波動(0.95~0.97上下波動)。請大牛們幫我分析一下到底哪里有問題。 in文件和data文件如下: in文件: units real dimension 3 boundary p p p atom_style full pair_style lj/cut/coul/long 9.0 9.0 bond_style harmonic angle_style harmonic kspace_style pppm 1.0e-4 read_data water.lammps05 neighbor 2.0 bin neigh_modify delay 5 every 1 velocity all create 298 4928459 dist uniform variable density equal mass(all)*10/(6.0221417930*vol) timestep 0.01 thermo 2000 thermo_style custom step temp press vol etotal v_density fix fixShake all shake 0.0001 100 0 b 1 a 1 fix 1 all npt temp 293.0 293.0 100 iso 1.0 1.0 1000 drag 0.2 run 100000 unfix 1 fix 2 all nve dump 1 all custom 1000 dumpnve.lammpstrj id type xu yu zu run 10000 data文件: LAMMPS 2005 data file for water 6000 atoms 4000 bonds 2000 angles 0 dihedrals 0 impropers 2 atom types 1 bond types 1 angle types 0.000000000 39.107300000 xlo xhi 0.000000000 39.107300000 ylo yhi 0.000000000 39.107300000 zlo zhi Masses 1 15.999400 2 1.007970 Pair Coeffs 1 0.1553000000 3.1660000000 2 0.0000000000 0.0000000000 Bond Coeffs 1 450 1 Angle Coeffs 1 47.0000 109.4700 Atoms 1 1 1 -0.82 7.597767833 9.408946512 5.391758138 2 1 2 0.41 6.641116667 9.263121029 5.557195421 3 1 2 0.41 7.756749235 8.739511125 4.691137645 …… |
銅蟲 (小有名氣)

新蟲 (初入文壇)
銀蟲 (小有名氣)
新蟲 (小有名氣)
金蟲 (正式寫手)
新蟲 (小有名氣)
新蟲 (小有名氣)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +8 | 梨花珞晚風(fēng) 2026-03-17 | 8/400 |
|
|
[考研] 招08考數(shù)學(xué) +3 | laoshidan 2026-03-20 | 7/350 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|