| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1688 | 回復(fù): 4 | ||||
wzjing103新蟲(chóng) (初入文壇)
|
[求助]
Lammps模擬苯酚和水在PDMS膜中的擴(kuò)散,msd沒(méi)有增加
|
|
我做水和苯酚在聚二甲基硅氧烷膜中的擴(kuò)散對(duì)比,以下是in文件中設(shè)置的部分,為什么水和苯酚的msd根本沒(méi)有增加,只是稍微有波動(dòng)?(本人是初學(xué)者,剛接觸幾個(gè)月) 體系中有3個(gè)聚二甲基硅氧烷分子,苯酚和水各五個(gè) 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 |

專家顧問(wèn) (正式寫(xiě)手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
很難說(shuō),你的設(shè)置看上去有點(diǎn)奇怪。正常體系的MD模擬,一般是先minimize,然后NVT下平衡,然后NVE下采樣。 而你的步驟,在minimize之后,卻用了fix langevin + fix nve 你這樣做有什么特殊的理由嗎?是別人文獻(xiàn)里也這么做的? fix langevin一般是對(duì)隱式溶劑(implicit solvent)下的模擬使用的,也就是說(shuō)在模擬中沒(méi)有水分子,所有水分子被當(dāng)成均一的介質(zhì),用隨機(jī)力來(lái)模擬溶質(zhì)分子與水分子的碰撞。但我看你的pair_style里有 lj/cut/tip4p,看起來(lái)是有水分子存在的,所以我覺(jué)得用fix langevin不是非常合理。 另外, fix langevin不能和fix nvt混用,所以你在unfix nve后也需要unfix langevin。 你的MSD計(jì)算是怎么設(shè)置的? |
新蟲(chóng) (初入文壇)
|
我是想模擬水和苯酚的布朗運(yùn)動(dòng)所以用了fix langevin,我對(duì)模擬的過(guò)程了解的不太多,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 |

禁蟲(chóng) (小有名氣)
|
本帖內(nèi)容被屏蔽 |
銀蟲(chóng) (小有名氣)
|
我也遇到同樣的問(wèn)題,MSD的運(yùn)算步長(zhǎng)沒(méi)達(dá)到當(dāng)初設(shè)計(jì)的值。 能否看看你的in文件,不勝感激。 zzz5170@163.com |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +12 | 小熊~~ 2026-03-31 | 12/600 |
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
[考研] 314求調(diào)劑 +4 | 溪云珂 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 圖鑒212 2026-03-30 | 5/250 |
|
|
[考研] 085410人工智能 初試316分 求調(diào)劑 +3 | 殘星拂曙 2026-03-31 | 3/150 |
|
|
[考研] 土木304求調(diào)劑 +3 | 兔突突突, 2026-03-31 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +10 | @taotao 2026-03-29 | 10/500 |
|
|
[考研] 【調(diào)劑】一志愿廈大生物與醫(yī)藥調(diào)劑 +3 | Echo蝦米 2026-03-31 | 3/150 |
|
|
[考研] 生物學(xué)308分求調(diào)劑(一志愿華東師大) +3 | 相信必會(huì)光芒萬(wàn)?/a> 2026-03-31 | 3/150 |
|
|
[考研] 求調(diào)劑:一志愿:南京大學(xué) 專業(yè):0705 總分320 ,本科985,四六級(jí)已過(guò) +3 | lfy760306 2026-03-31 | 3/150 |
|
|
[考研] 311(085601)求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
[考博] 26申博 +4 | 加油沖! 2026-03-26 | 4/200 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|