 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1690 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
wzjing103新蟲 (初入文壇)
|
[求助]
Lammps模擬苯酚和水在PDMS膜中的擴散，msd沒有增加
|
|||
|
我做水和苯酚在聚二甲基硅氧烷膜中的擴散對比，以下是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ù)此樓» 猜你喜歡
805600專碩材料與化工348分求調(diào)劑
已經(jīng)有6人回復(fù)
-
重慶大學(xué)材料與化工085600，初試370+，求求調(diào)劑建議
已經(jīng)有5人回復(fù)
-
309求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北交大材料工程，總分358
已經(jīng)有4人回復(fù)
-
302求調(diào)劑一志愿北航070300，本科鄭大化學(xué)
已經(jīng)有11人回復(fù)
-
化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有4人回復(fù)
-
336求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助LAMMPS關(guān)于MSD計算的問題
已經(jīng)有17人回復(fù)
wzjing103
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 19.2
- 帖子: 18
- 在線: 8.4小時
- 蟲號: 3621957
- 注冊: 2014-12-31
- 性別: GG
- 專業(yè): 污染控制化學(xué)
|
我是想模擬水和苯酚的布朗運動所以用了fix langevin，我對模擬的過程了解的不太多，msd計算用的是 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
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗: +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時
- 蟲號: 556471
- 注冊: 2008-05-11
- 性別: GG
- 專業(yè): 理論和計算化學(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è)置看上去有點奇怪。正常體系的MD模擬，一般是先minimize，然后NVT下平衡，然后NVE下采樣。 而你的步驟，在minimize之后，卻用了fix langevin + fix nve 你這樣做有什么特殊的理由嗎？是別人文獻里也這么做的？ fix langevin一般是對隱式溶劑(implicit solvent）下的模擬使用的，也就是說在模擬中沒有水分子，所有水分子被當(dāng)成均一的介質(zhì)，用隨機力來模擬溶質(zhì)分子與水分子的碰撞。但我看你的pair_style里有 lj/cut/tip4p，看起來是有水分子存在的，所以我覺得用fix langevin不是非常合理。 另外, fix langevin不能和fix nvt混用，所以你在unfix nve后也需要unfix langevin。 你的MSD計算是怎么設(shè)置的？ |
2樓2015-04-27 13:26:14
silence-jin
禁蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
4樓2015-12-16 21:36:20
5樓2015-12-18 19:02:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +5 | dianeeee 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +7 | 雙馬尾痞老板2 2026-04-01 | 7/350 |
|
|
[考研] 270求調(diào)劑 +7 | 小杰pp 2026-03-31 | 8/400 |
|
|
|
[考研] 261求調(diào)劑 +3 | 明仔· 2026-04-01 | 3/150 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +7 | 哇呼哼呼哼 2026-03-26 | 7/350 |
|
|
|
[考研] 310分求調(diào)劑 +4 | 成功上岸wang 2026-04-01 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
|
[考研] 0856 335分 +9 | cccchenso 2026-03-29 | 9/450 |
|
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +4 | 陳陽坤 2026-03-31 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 求佛 +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|