 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1691 | 回復(fù): 4 | |||||
wzjing103新蟲 (初入文壇)
|
[求助]
Lammps模擬苯酚和水在PDMS膜中的擴(kuò)散，msd沒有增加
|
|
我做水和苯酚在聚二甲基硅氧烷膜中的擴(kuò)散對比，以下是in文件中設(shè)置的部分，為什么水和苯酚的msd根本沒有增加，只是稍微有波動？（本人是初學(xué)者，剛接觸幾個月） 體系中有3個聚二甲基硅氧烷分子，苯酚和水各五個 pair_style hybrid lj/charmm/coul/charmm 8.0 10.0 lj/cut/tip4p/long 6 7 5 8 0.15 12.0 bond_style harmonic angle_style charmm kspace_style pppm/tip4p 0.001 neighbor 2.0 bin neigh_modify delay 0 every 1 check yes timestep 1.0 min_style fire minimize 1.0e-16 1.0e-16 50000 100000 fix limit all nve fix speed all langevin 293 293 10 482799 run 500000 unfix limit fix relax all nvt temp 293 293 100 run 1000000 |
回復(fù)此樓» 猜你喜歡
283分材料與化工求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料調(diào)劑
已經(jīng)有6人回復(fù)
-
材料340分調(diào)劑
已經(jīng)有6人回復(fù)
-
321求調(diào)劑 一志愿 浙江工業(yè)大學(xué)生物醫(yī)藥
已經(jīng)有6人回復(fù)
-
085801 總分275 本科新能源 求調(diào)劑
已經(jīng)有19人回復(fù)
-
285求調(diào)劑
已經(jīng)有13人回復(fù)
-
280求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有27人回復(fù)
-
071000生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
26考研調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助LAMMPS關(guān)于MSD計(jì)算的問題
已經(jīng)有17人回復(fù)
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wzjing103: 金幣+2, ★有幫助 2015-04-27 16:13:19
wzjing103(abinitio代發(fā)): 金幣+2, 鼓勵交流 2015-04-28 07:42:17
感謝參與，應(yīng)助指數(shù) +1
wzjing103: 金幣+2, ★有幫助 2015-04-27 16:13:19
wzjing103(abinitio代發(fā)): 金幣+2, 鼓勵交流 2015-04-28 07:42:17
|
很難說，你的設(shè)置看上去有點(diǎn)奇怪。正常體系的MD模擬，一般是先minimize，然后NVT下平衡，然后NVE下采樣。 而你的步驟，在minimize之后，卻用了fix langevin + fix nve 你這樣做有什么特殊的理由嗎？是別人文獻(xiàn)里也這么做的？ fix langevin一般是對隱式溶劑(implicit solvent）下的模擬使用的，也就是說在模擬中沒有水分子，所有水分子被當(dāng)成均一的介質(zhì)，用隨機(jī)力來模擬溶質(zhì)分子與水分子的碰撞。但我看你的pair_style里有 lj/cut/tip4p，看起來是有水分子存在的，所以我覺得用fix langevin不是非常合理。 另外, fix langevin不能和fix nvt混用，所以你在unfix nve后也需要unfix langevin。 你的MSD計(jì)算是怎么設(shè)置的？ |
2樓2015-04-27 13:26:14
wzjing103
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 19.2
- 帖子: 18
- 在線: 8.4小時
- 蟲號: 3621957
- 注冊: 2014-12-31
- 性別: GG
- 專業(yè): 污染控制化學(xué)
|
我是想模擬水和苯酚的布朗運(yùn)動所以用了fix langevin，我對模擬的過程了解的不太多，msd計(jì)算用的是 compute m1 phenol msd/molecule compute m2 h2o msd/molecule variable msd1 equal c_m1[1][4] variable msd2 equal c_m1[2][4] variable msd3 equal c_m1[3][4] variable msd4 equal c_m1[4][4] variable msd5 equal c_m1[5][4] variable msd6 equal c_m2[1][4] variable msd7 equal c_m2[2][4] variable msd8 equal c_m2[3][4] variable msd9 equal c_m2[4][4] variable msd10 equal c_m2[5][4] variable step equal step fix msd all print 10000 "${step} ${msd1} ${msd2} ${msd3} ${msd4} ${msd5} ${msd6} ${msd7} ${msd8} ${msd9} ${msd10}" file all.msd screen no |
3樓2015-04-27 14:37:43
silence-jin
禁蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
4樓2015-12-16 21:36:20
5樓2015-12-18 19:02:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wzjing103 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283分材料與化工求調(diào)劑 +10 | 羅KAKA 2026-04-02 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 312 化工或制藥調(diào)劑 +3 | 小小墨123 2026-04-02 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +9 | 汲亦昊 2026-03-29 | 9/450 |
|
|
|
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
|
[考研]
 一志愿廈門大學(xué)化學(xué)工程（專碩）-數(shù)二英二406分-求調(diào)劑
+5
|
廈大化工 2026-04-01 | 5/250 |
|
|
|
[考研] 0710生物學(xué)，325求調(diào)劑 +3 | mkkkkkl 2026-04-01 | 3/150 |
|
|
|
[考研] 379求調(diào)劑 +3 | ?苦瓜不苦 2026-04-01 | 3/150 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +10 | 腳滑的守法公民 2026-03-27 | 10/500 |
|
|
|
[考研] 0703一志愿南師大334求調(diào)劑 +4 | seven7yu 2026-03-30 | 4/200 |
|
|
|
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +9 | 袁sy 2026-04-01 | 9/450 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[考研] 349求調(diào)劑 +6 | zwjjjjjj 2026-03-31 | 6/300 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
|
[考研] 086000調(diào)劑 +5 | 7901117076 2026-03-26 | 5/250 |
|
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +3 | Acesczlo 2026-03-29 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|