 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1692 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
wzjing103新蟲 (初入文壇)
|
[求助]
Lammps模擬苯酚和水在PDMS膜中的擴(kuò)散，msd沒有增加
|
|||
|
我做水和苯酚在聚二甲基硅氧烷膜中的擴(kuò)散對(duì)比，以下是in文件中設(shè)置的部分，為什么水和苯酚的msd根本沒有增加，只是稍微有波動(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 |
回復(fù)此樓» 猜你喜歡
362求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑 302分初試 0854
已經(jīng)有4人回復(fù)
-
299求調(diào)劑
已經(jīng)有4人回復(fù)
-
266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有4人回復(fù)
-
315求調(diào)劑
已經(jīng)有7人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有3人回復(fù)
-
316求調(diào)劑
已經(jīng)有16人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有13人回復(fù)
-
材料與化工306分找調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助LAMMPS關(guān)于MSD計(jì)算的問題
已經(jīng)有17人回復(fù)
silence-jin
禁蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
4樓2015-12-16 21:36:20
lsloneil
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +49 - 模擬EPI: 2
- 應(yīng)助: 148 (高中生)
- 金幣: 8640.9
- 散金: 507
- 紅花: 30
- 帖子: 505
- 在線: 598.4小時(shí)
- 蟲號(hào): 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, 鼓勵(lì)交流 2015-04-28 07:42:17
感謝參與，應(yīng)助指數(shù) +1
wzjing103: 金幣+2, ★有幫助 2015-04-27 16:13:19
wzjing103(abinitio代發(fā)): 金幣+2, 鼓勵(lì)交流 2015-04-28 07:42:17
|
很難說，你的設(shè)置看上去有點(diǎn)奇怪。正常體系的MD模擬，一般是先minimize，然后NVT下平衡，然后NVE下采樣。 而你的步驟，在minimize之后，卻用了fix langevin + fix nve 你這樣做有什么特殊的理由嗎？是別人文獻(xiàn)里也這么做的？ fix langevin一般是對(duì)隱式溶劑(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小時(shí)
- 蟲號(hào): 3621957
- 注冊: 2014-12-31
- 性別: GG
- 專業(yè): 污染控制化學(xué)
|
我是想模擬水和苯酚的布朗運(yùn)動(dòng)所以用了fix langevin，我對(duì)模擬的過程了解的不太多，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
5樓2015-12-18 19:02:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請]
請問共同通訊和共同一作的認(rèn)可度問題
 10+4
|
psa1234 2026-04-01 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) +4 | goldtt 2026-04-02 | 6/300 |
|
|
|
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +5 | 小江小江江江 2026-04-02 | 5/250 |
|
|
|
[考研] 085600，320分求調(diào)劑 +5 | 大饞小子 2026-04-02 | 5/250 |
|
|
|
[考研] 一志愿同濟(jì)大學(xué)323分（080500）求調(diào)劑 +6 | yikeniu 2026-04-01 | 6/300 |
|
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +6 | 1段時(shí)間 2026-03-31 | 7/350 |
|
|
|
[考研] 292分，材料與化工，申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-27 | 7/350 |
|
|
|
[考研] 270調(diào)劑 +7 | maxjxbsk 2026-04-02 | 7/350 |
|
|
|
[考研] 303求調(diào)劑 +8 | DLkz1314. 2026-03-30 | 8/400 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +15 | minht 2026-03-29 | 16/800 |
|
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|