 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 6050 | 回復(fù): 44 | ||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)4次，作者043114076增加金幣 3.2 個(gè) | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[資源]
關(guān)于towhee 回答zhang_jaj的問(wèn)題
|
||||
|
以下是zhang_jaj 的問(wèn)題: 你好，我想用towhee做NVT下argon氣在surface面上nucleation過(guò)程，希望你給一些指導(dǎo)，首先我想把potential選為L(zhǎng)J，surface與particles也是LJ作用， 1）問(wèn)題是這個(gè)surface是需要用particles來(lái)make嗎，如果是的話那么構(gòu)成它的particles該怎么固定 2 ）或者，直接用towhee里面的LJ9-3 wall嗎，想知道這種wall是與particles怎么作用的，它不是用小粒子構(gòu)成的嗎？ 3）模擬完成后，液體與固體的結(jié)構(gòu)判斷該怎么做，有什么程序代碼嗎？ 4）這個(gè)nucleation的rate該怎么算（是單位時(shí)間單位體積內(nèi)產(chǎn)生臨街晶核的個(gè)數(shù)），我想知道程序中怎么實(shí)現(xiàn)，我沒(méi)發(fā)現(xiàn)有那個(gè)參數(shù)相關(guān)的，或者是需要自己來(lái)加?xùn)|西嗎？ 5）還有這個(gè)nucleation barrier該怎么來(lái)實(shí)現(xiàn)？ 6）towhee不能并行，那么其實(shí)一個(gè)任務(wù)就是在一個(gè)process上處理的；如果我的個(gè)人電腦比服務(wù)器的主頻高，那么是不是用臺(tái)式機(jī)比服務(wù)器更快呢（因?yàn)榕_(tái)式機(jī)是4核的，它是不是能自己把towhee的任務(wù)來(lái)用4核一起處理？）？ 7）希望您能給點(diǎn)關(guān)于towhee代碼學(xué)習(xí)的經(jīng)驗(yàn)，就是對(duì)于我要做的這個(gè)方向的。 這是我的回答: 實(shí)在不敢當(dāng), 我現(xiàn)在也是學(xué)習(xí)階段, 因?yàn)闆](méi)做過(guò)類似的體系, 對(duì)于你的問(wèn)題只能有一些建議. 1）問(wèn)題是這個(gè)surface是需要用particles來(lái)make嗎，如果是的話那么構(gòu)成它的particles該怎么固定 可以使用particle模型來(lái)模擬, 固定比較簡(jiǎn)單, 在move設(shè)置中不讓這些粒子移動(dòng)就行了. 2 ）或者，直接用towhee里面的LJ9-3 wall嗎，想知道這種wall是與particles怎么作用的，它不是用小粒子構(gòu)成的嗎？ 這個(gè)我以前沒(méi)有用過(guò), 不過(guò)文獻(xiàn)當(dāng)中確實(shí)有人這樣做. 這種勢(shì)能的使用情況可以參考例子Fris_Walls, 另外上網(wǎng)搜了一下, 這幾個(gè)帖子和文章可能對(duì)你有用. http://jcp.aip.org/resource/1/jcpsa6/v67/i5/p2384_s1 http://jcp.aip.org/resource/1/jcpsa6/v73/i8/p4050_s1 http://lammps.sandia.gov/threads/msg24356.html http://www.sklogwiki.org/SklogWi ... ard-Jones_potential 這種勢(shì)是粒子勢(shì)能轉(zhuǎn)化而來(lái)的, 第四個(gè)參考文獻(xiàn)有推導(dǎo)過(guò)程, 主要的優(yōu)勢(shì)是可以大大提高速度. lammps 的帖子跟你的情況有點(diǎn)像, 使用MD做的, 可以借鑒. 3）模擬完成后，液體與固體的結(jié)構(gòu)判斷該怎么做，有什么程序代碼嗎？ 這個(gè)主要是分析RDF吧, 代碼在 utils 中有, 可以根據(jù)你的情況進(jìn)行修改. 4）這個(gè)nucleation的rate該怎么算（是單位時(shí)間單位體積內(nèi)產(chǎn)生臨街晶核的個(gè)數(shù)），我想知道程序中怎么實(shí)現(xiàn)，我沒(méi)發(fā)現(xiàn)有那個(gè)參數(shù)相關(guān)的，或者是需要自己來(lái)加?xùn)|西嗎？ 這個(gè)是需要后處理才能得到的吧. towhee本身不是計(jì)算成核過(guò)程的. 應(yīng)該需要編寫代碼實(shí)現(xiàn). 不過(guò)MC 如何計(jì)算時(shí)間呢? 如果用MC步數(shù)是不是會(huì)有問(wèn)題? 是不是用MD更合適一點(diǎn)? 5）還有這個(gè)nucleation barrier該怎么來(lái)實(shí)現(xiàn)？ 這個(gè)應(yīng)該是進(jìn)行不同溫度點(diǎn)的模擬, 作速率-溫度關(guān)系圖, 然后通過(guò)擬合熱力學(xué)公式得到吧. 6）towhee不能并行，那么其實(shí)一個(gè)任務(wù)就是在一個(gè)process上處理的；如果我的個(gè)人電腦比服務(wù)器的主頻高，那么是不是用臺(tái)式機(jī)比服務(wù)器更快呢（因?yàn)榕_(tái)式機(jī)是4核的，它是不是能自己把towhee的任務(wù)來(lái)用4核一起處理？）？ towhee的不能用多個(gè)核并行計(jì)算一個(gè)任務(wù), 而事實(shí)上也不需要. 一般來(lái)說(shuō), 模擬多個(gè)任務(wù)可以等價(jià)于模擬一個(gè)大任務(wù), MC的本質(zhì)決定它是這樣的性質(zhì), MD可不行. 所以盡量還是使用towhee提供的并行模式. 注意選擇不同的隨機(jī)種子. 7）希望您能給點(diǎn)關(guān)于towhee代碼學(xué)習(xí)的經(jīng)驗(yàn)，就是對(duì)于我要做的這個(gè)方向的。 關(guān)于代碼, 其實(shí)我也沒(méi)有太多經(jīng)驗(yàn). 我主要是對(duì)結(jié)果進(jìn)行后處理. towhee 本身是一個(gè)general Monte Carlo code, 不是專門為某些體系設(shè)計(jì)的, 但是對(duì)于你的體系應(yīng)該是可以勝任的. 主要針對(duì)體系的代碼, 差別也主要是建模和分析結(jié)果兩個(gè)過(guò)程, 核心代碼都是類似的. 我個(gè)人認(rèn)為如果時(shí)間有限, 不要陷到源代碼中去.當(dāng)然能夠多了解一些代碼對(duì)于軟件的使用也是非常有好處的. towhee 的代碼由于需要我簡(jiǎn)單的讀過(guò)一些, 還是比較清楚的. 可以嘗試從總體上去學(xué)習(xí)下. 對(duì)于結(jié)果的后處理的代碼, 建議從修改towhee 提供的代碼開(kāi)始. 使用語(yǔ)言以自己最熟悉的為準(zhǔn), 如果都湊合, 推薦 fortran90. 多學(xué)習(xí) linux bash 腳本也是很有好處的, 我現(xiàn)在簡(jiǎn)單的數(shù)據(jù)處理基本都用它. 最后強(qiáng)調(diào), 不要陷到代碼中去, 能夠完成自己的用途就可以了, 越簡(jiǎn)單越直接越好. 多多思考跟你體系相關(guān)的問(wèn)題, 多查文獻(xiàn). ///希望能多多向你學(xué)習(xí)towhee知識(shí)，導(dǎo)師不管，自己一個(gè)人摸索。擔(dān)心畢業(yè)十大問(wèn)題啊，謝謝了 大家基本都差不多, 堅(jiān)定信心, 可以做的很好. 其實(shí)模擬計(jì)算這個(gè)方向比起做實(shí)驗(yàn)還是比較透明的, 重復(fù)文獻(xiàn)一般沒(méi)有問(wèn)題, 技術(shù)問(wèn)題一般論壇上都會(huì)回答. 我也是在學(xué)習(xí)的過(guò)程, 大家多多討論, 相互學(xué)習(xí). [ Last edited by 043114076 on 2012-4-26 at 09:55 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 towhee精華 | towhee的運(yùn)用討論��！ |
» 本帖已獲得的紅花（最新10朵）
zhang_jaj» 猜你喜歡
0856求調(diào)劑
已經(jīng)有13人回復(fù)
-
0703化學(xué)
已經(jīng)有11人回復(fù)
-
329求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有8人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有12人回復(fù)
-
309求調(diào)劑
已經(jīng)有3人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有4人回復(fù)
-
330分求調(diào)劑
已經(jīng)有3人回復(fù)
-
299求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于分子模塊求助的一些建議
已經(jīng)有29人回復(fù)
- towhee軟件的安裝問(wèn)題
已經(jīng)有4人回復(fù)
- Towhee能模擬氣體吸附等溫線嗎？
已經(jīng)有10人回復(fù)
- towhee編譯失敗
已經(jīng)有9人回復(fù)
- 運(yùn)行 TOWHEE 時(shí)出現(xiàn)的問(wèn)題 ，希望有人幫我解決一下，非常感謝
已經(jīng)有8人回復(fù)
- 使用amber力場(chǎng)，一些參數(shù)無(wú)法獲得
已經(jīng)有5人回復(fù)
- 想做粒子和固體（金屬）相互作用，不知道lammps能否實(shí)現(xiàn)？
已經(jīng)有6人回復(fù)
- 求問(wèn)分子模擬斷檢問(wèn)題。。。
已經(jīng)有5人回復(fù)
- LAMMPS并行安裝問(wèn)題
已經(jīng)有4人回復(fù)
- MS 進(jìn)行溶液分子動(dòng)力學(xué)計(jì)算時(shí)在哪里設(shè)置力場(chǎng)參數(shù)
已經(jīng)有5人回復(fù)
- 關(guān)于氣體液化的問(wèn)題~~~~~
已經(jīng)有13人回復(fù)
- 求用到AMBER force field的開(kāi)源軟件
已經(jīng)有10人回復(fù)
- 水分子在外場(chǎng)下鍵角、鍵長(zhǎng)的變化？
已經(jīng)有4人回復(fù)
- ms模型導(dǎo)入lammps
已經(jīng)有8人回復(fù)
- lammps里的隨機(jī)種子是如何確定的？
已經(jīng)有12人回復(fù)
- lammps里log文件里一些數(shù)據(jù)的單位？
已經(jīng)有10人回復(fù)
- 使用Towhee做NVT模擬時(shí)，氣液兩盒子的hmatrix 大小該怎么設(shè)置？ 有什么原則或公式么？
已經(jīng)有7人回復(fù)
- 關(guān)于氣液平衡的簡(jiǎn)單問(wèn)題
已經(jīng)有24人回復(fù)
- 石墨烯 lammps
已經(jīng)有11人回復(fù)
- 石墨烯和下層的SiO2會(huì)成鍵嗎？
已經(jīng)有7人回復(fù)
- 使用towhee進(jìn)行GEMC模擬相平衡的MC運(yùn)動(dòng)定義問(wèn)題..
已經(jīng)有15人回復(fù)
- 【討論】關(guān)于爆炸鹽/浴鹽/沐浴球的配方
已經(jīng)有4人回復(fù)
- 【求助】弱問(wèn)有關(guān)AMBER遠(yuǎn)程調(diào)用的問(wèn)題？
已經(jīng)有3人回復(fù)
- 【求助】分子的LJ勢(shì)能參數(shù)的擬合
已經(jīng)有8人回復(fù)
|
你好，您知道towhee用intel的編譯器好呢還是PGI的好，因?yàn)槲铱垂俜浇o的是PGI的。還有我在運(yùn)行例子的時(shí)候，提交大于3個(gè)任務(wù)并行時(shí)（如16），從倒數(shù)第三個(gè)就會(huì)顯示出 Job # 14 Insufficient entres in towhee_parallel file 后面就沒(méi)了，但是前13個(gè)job都順利。不知道是怎么回事。 當(dāng)小于三時(shí)也會(huì)有最后一個(gè)job顯示Insufficient entres in towhee_parallel file 還有想問(wèn)下，mpich2的并行環(huán)境好像不能編譯是嗎？ 還有想問(wèn)下，如果在./configrue時(shí)，不給CC，F(xiàn)77等賦值，編譯好的不會(huì)是用GNU的編譯器吧？ 麻煩了，不甚感激！ |
5樓2012-05-04 22:00:07
2樓2012-04-26 19:15:47
|
你好，我是按照你之前一個(gè)帖子直接編譯的并行towhee-7版的，但是這個(gè)里面沒(méi)有那個(gè)towheebase的目錄（應(yīng)該是6版以前才有的把）。在做練習(xí)的時(shí)候，我把towhee_input文件的forcefield改成了自己的絕對(duì)路徑 /public/software/towhee/Forcefield/() 然后towhee_parallel如下 1 towhee.log /public/software/towhee/Examples/Canonical/Amber_Villin 之后執(zhí)行了mpirun -np 4 towhee 顯示出來(lái)的為 Number ofjobs:0 Simulations complete Total time:0.000059 seconds 而且沒(méi)有towhee.log文件 還有官網(wǎng)說(shuō)安裝需要的autoconf與automake文件是不是可以不用？ example中的makefile文件用來(lái)編譯什么的？ 非常感謝�。ㄒ�?yàn)橐恍┰驅(qū)嶒?yàn)室不能上網(wǎng)，只能晚上回來(lái)求助） [ Last edited by zhang_jaj on 2012-4-26 at 20:25 ] |
3樓2012-04-26 20:24:20
4樓2012-04-26 21:53:51
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +7 | 溫柔的晚安 2026-03-25 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|