 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1069 | 回復(fù): 2 | ||
[求助]
關(guān)于GEMC模擬相平衡的一些問題。
|
|
我使用的是MCCCS Towhee來模擬純物質(zhì)的共存曲線，每次模擬還是碰到很多類似問題，現(xiàn)在想就我碰到的問題請教一下大家： 1、體積的設(shè)定：兩相體積的設(shè)定，文獻Sample Text給的建議是總體積的設(shè)定最好保證20%的粒子數(shù)目在氣相。不知道大家是怎么看這個建議的。在低溫的時候20%的粒子在氣相會造成氣相的體積非常之大，有的時候汽液相體積之比在10以上了這樣做會不會影響模擬？ 2、各種概率的設(shè)置：很多文獻（Ilja Siepm ann的系列文章中都有提到）上寫的是體積和粒子交換（swap）這兩個移動概率要保證在5~10個cycles之內(nèi)分別有一次被接受； 我是這樣做的，體積概率 P 的設(shè)置：N（分子數(shù)）*nsteps（循環(huán)數(shù)）*0.5（接受概率）*P（體積變化的概率）=nsteps（循環(huán)數(shù)）/（5~10）； swap的概率這個起初我是隨便給的0.3左右，先跑個10000cycles再看看接受的swap數(shù)目{兩部分（box1到box2+box2到box1）}的大小，再根據(jù)情況調(diào)試，比如：10000cycles有3000次接受。說明swap概率大了，需要適當(dāng)調(diào)��；500次被接受說明需要適當(dāng)增大；剩下的概率均分給其他的移動。但是，問題出現(xiàn)在低溫的情況，低溫的swap就算給到了65%也滿足不了5~10個cycles之內(nèi)分別有一次被接受，并且有些文獻有說，除去體積和swap的概率，剩下的概率要不要小于50%。不知道大家是怎么分配的。 3、模擬的平衡與否的判斷條件：看文獻寫的是有明顯的不變的汽液密度，說明達到了平衡。但是我模擬的時候汽液密度一直是波動的，而且波動還是挺大的，我想是不是我的步數(shù)不夠，就一直接著跑，但是看起過來他就是達不到平衡（收斂）。文獻上都說的一般都是10000~25000達到平衡，接著另外20000~40000產(chǎn)出。我都快10w步還沒看到平衡。。。。。。一開始接觸MC時候，不是說看體系的總能量減小嗎，Metropolis抽樣法不是朝著能量小的構(gòu)型去構(gòu)建馬爾科夫鏈的嗎?為什么結(jié)果總是看到總能量（box1+box2）在不斷的減小又增大（有時差一個數(shù)量級）。 4、總覺得小木蟲上很少有關(guān)于分子模擬在汽液平衡這方面的交流，不知道隊伍啥時候能壯大.... |
回復(fù)此樓» 猜你喜歡
求材料調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有6人回復(fù)
-
286求調(diào)劑
已經(jīng)有9人回復(fù)
-
焦慮
已經(jīng)有12人回復(fù)
-
344求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有9人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有18人回復(fù)
-
314求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 大學(xué)物化相平衡問題
已經(jīng)有4人回復(fù)
- Lammps計算油水界面，油相密度、平衡等問題求指導(dǎo)
已經(jīng)有6人回復(fù)
- 相平衡 關(guān)于半透膜的問題
已經(jīng)有10人回復(fù)
- 一些關(guān)于模擬汽液相平衡的問題（towhee新手）
已經(jīng)有4人回復(fù)
- 問一個關(guān)于聚合物體系平衡的問題
已經(jīng)有25人回復(fù)
- NPT系綜下判斷是否平衡的標志
已經(jīng)有12人回復(fù)
- 擴散系數(shù)的計算及體系平衡的判斷
已經(jīng)有4人回復(fù)
- Aspen Plus 軟件中氣液平衡問題
已經(jīng)有9人回復(fù)
- 徑向分布函數(shù)的問題
已經(jīng)有13人回復(fù)
- 使用Towhee做NVT模擬時，氣液兩盒子的hmatrix 大小該怎么設(shè)置？ 有什么原則或公式么？
已經(jīng)有7人回復(fù)
- 關(guān)于氣液平衡的簡單問題
已經(jīng)有24人回復(fù)
- 相平衡部分，小問題求助
已經(jīng)有4人回復(fù)
- 關(guān)于towhee 回答zhang_jaj的問題
已經(jīng)有44人回復(fù)
- topol文件不識別小分子
已經(jīng)有8人回復(fù)
- MS金屬氧化物表面吸附有機物平衡問題
已經(jīng)有4人回復(fù)
- 用熱力學(xué)計算的平衡轉(zhuǎn)化率只有實驗測得的反應(yīng)轉(zhuǎn)化率的一半。。計算是否會有問題
已經(jīng)有5人回復(fù)
- 使用towhee進行GEMC模擬相平衡的MC運動定義問題..
已經(jīng)有15人回復(fù)
- 【求助】關(guān)于相平衡數(shù)據(jù)的資料
已經(jīng)有4人回復(fù)
- 【求助】氣液平衡的數(shù)據(jù)
已經(jīng)有5人回復(fù)
2樓2014-07-23 09:28:37
草田露
銀蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 227.1
- 散金: 24
- 帖子: 26
- 在線: 60.3小時
- 蟲號: 3080499
- 注冊: 2014-03-22
- 性別: MM
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
3樓2015-09-07 09:04:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 992648352 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +6 | 隔壁陳先生 2026-03-12 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +3 | 想上岸的鯉魚 2026-03-18 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +3 | 1孫悟空 2026-03-17 | 3/150 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|