| 5 | 1/1 | 返回列表 |
| 查看: 3131 | 回復(fù): 15 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
hyff鐵蟲 (初入文壇)
|
[求助]
使用towhee進(jìn)行GEMC模擬相平衡的MC運(yùn)動定義問題..
|
|
|
towhee新手,想使用GEMC方法來模擬二元相平衡。 目前已經(jīng)建立了兩個(gè)盒子來模擬氣液兩相,BOX1初始為氣相分子200個(gè),BOX2初始為液相分子200個(gè)。 MC運(yùn)動函數(shù)部分是從example里的Amber_IsoPropanol這個(gè)基礎(chǔ)上嘗試改為二元體系的,根據(jù)說明將相關(guān)的pm2cbswmt,pmcbmt等都賦了兩個(gè)值,但試算之后發(fā)現(xiàn)只有BOX1中的氣相分子可以插入BOX2,BOX2中的液相分子無法插入BOX1,結(jié)果也就無法達(dá)到化學(xué)勢平衡。 求教使用過towhee進(jìn)行過GEMC模擬和相平衡模擬的高手指點(diǎn)交流.. |
木蟲 (正式寫手)

鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
|
貼一個(gè)我最新修改的,因?yàn)檫@個(gè)東西我們組之前完全沒人做過,我是第一個(gè)做的……然后都是按我自己的理解寫的,哪里不對的煩請指正一下。 pmvol 0.002d0 pmvlpr 1.0d0 rmvol 0.1d0 tavol 0.5d0 pm2boxcbswap 0.4d0 pm2cbswmt 0.5d0 1.0d0 pm2cbswpr 1.0d0 pmcb 0.6d0 pmcbmt 0.5d0 1.0d0 pmall 0.0d0 0.0d0 pmtracm 0.8d0 pmtcmt 0.5d0 1.0d0 rmtrac 0.5d0 tatrac 0.5d0 pmrotate 1.0d0 pmromt 0.5d0 1.0d0 rmrot 0.05d0 tarot 0.5d0 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|