 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2540 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
fengshaxia21鐵蟲(chóng) (初入文壇)
|
[求助]
SPC水分子模型，求密度有些許問(wèn)題 已有1人參與
|
|
|
RT，我用MS建立了一個(gè)6000分子的SPC水分子模型。但是得到的在293K的密度總是0.96這樣，距離實(shí)際值或者SPC水分子的實(shí)驗(yàn)值都有0.2～0.3的差距（雖然看上去不算大，不過(guò)這沒(méi)關(guān)系嗎？） 我想知道這是什么導(dǎo)致的，density=mass(all)/(參數(shù)*vol)，但是mass是不變的，密度的體積就只跟體積有關(guān)了。為什么在溫度到達(dá)293附近的時(shí)候波動(dòng)還比較大（270～300上下波動(dòng)），密度的波動(dòng)（0.95～0.97上下波動(dòng)）。請(qǐng)大牛們幫我分析一下到底哪里有問(wèn)題。 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 …… |
回復(fù)此樓» 猜你喜歡
070300化學(xué)求調(diào)劑
已經(jīng)有9人回復(fù)
-
085404 22408 315分
已經(jīng)有6人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有22人回復(fù)
-
347求調(diào)劑
已經(jīng)有10人回復(fù)
-
354求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿a區(qū)211，085601-307分求調(diào)劑
已經(jīng)有11人回復(fù)
-
284求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有18人回復(fù)
-
362求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- TIP3P 和 SPC/E 適用的體系
已經(jīng)有7人回復(fù)
cuso4852cuoh
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 223.4
- 散金: 20
- 帖子: 170
- 在線: 49.6小時(shí)
- 蟲(chóng)號(hào): 1481571
- 注冊(cè): 2011-11-08
- 專(zhuān)業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
9樓2018-12-27 21:18:42
liuasia
銅蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 102.6
- 散金: 200
- 帖子: 69
- 在線: 122.6小時(shí)
- 蟲(chóng)號(hào): 3001880
- 注冊(cè): 2014-02-27
- 性別: GG
- 專(zhuān)業(yè): 藝術(shù)美學(xué)
2樓2014-04-09 16:10:49
laughingjudy
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 185.8
- 帖子: 33
- 在線: 31.5小時(shí)
- 蟲(chóng)號(hào): 2738781
- 注冊(cè): 2013-10-20
- 專(zhuān)業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
3樓2014-06-13 14:59:07
4樓2014-06-18 21:18:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +9 | 柳溪邊 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 安徽大學(xué)專(zhuān)碩生物與醫(yī)藥專(zhuān)業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +5 | 美味可樂(lè)雞翅 2026-03-26 | 5/250 |
|
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +3 | 鐵憨憨123425 2026-03-31 | 3/150 |
|
|
|
[考研] 看看我！301求調(diào)劑 081700 +3 | 一個(gè)小明同學(xué) 2026-03-24 | 3/150 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 材料專(zhuān)碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專(zhuān)業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|