 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 6059 | 回復(fù): 44 | ||||
| 【獎(jiǎng)勵(lì)】 本帖被評價(jià)4次，作者043114076增加金幣 3.2 個(gè) | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[資源]
關(guān)于towhee 回答zhang_jaj的問題
|
||||
|
以下是zhang_jaj 的問題: 你好，我想用towhee做NVT下argon氣在surface面上nucleation過程，希望你給一些指導(dǎo)，首先我想把potential選為LJ，surface與particles也是LJ作用， 1）問題是這個(gè)surface是需要用particles來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)，我沒發(fā)現(xiàn)有那個(gè)參數(shù)相關(guān)的，或者是需要自己來加?xùn)|西嗎？ 5）還有這個(gè)nucleation barrier該怎么來實(shí)現(xiàn)？ 6）towhee不能并行，那么其實(shí)一個(gè)任務(wù)就是在一個(gè)process上處理的；如果我的個(gè)人電腦比服務(wù)器的主頻高，那么是不是用臺式機(jī)比服務(wù)器更快呢（因?yàn)榕_式機(jī)是4核的，它是不是能自己把towhee的任務(wù)來用4核一起處理？）？ 7）希望您能給點(diǎn)關(guān)于towhee代碼學(xué)習(xí)的經(jīng)驗(yàn)，就是對于我要做的這個(gè)方向的。 這是我的回答: 實(shí)在不敢當(dāng), 我現(xiàn)在也是學(xué)習(xí)階段, 因?yàn)闆]做過類似的體系, 對于你的問題只能有一些建議. 1）問題是這個(gè)surface是需要用particles來make嗎，如果是的話那么構(gòu)成它的particles該怎么固定 可以使用particle模型來模擬, 固定比較簡單, 在move設(shè)置中不讓這些粒子移動就行了. 2 ）或者，直接用towhee里面的LJ9-3 wall嗎，想知道這種wall是與particles怎么作用的，它不是用小粒子構(gòu)成的嗎？ 這個(gè)我以前沒有用過, 不過文獻(xiàn)當(dāng)中確實(shí)有人這樣做. 這種勢能的使用情況可以參考例子Fris_Walls, 另外上網(wǎng)搜了一下, 這幾個(gè)帖子和文章可能對你有用. 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 這種勢是粒子勢能轉(zhuǎn)化而來的, 第四個(gè)參考文獻(xiàn)有推導(dǎo)過程, 主要的優(yōu)勢是可以大大提高速度. 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)，我沒發(fā)現(xiàn)有那個(gè)參數(shù)相關(guān)的，或者是需要自己來加?xùn)|西嗎？ 這個(gè)是需要后處理才能得到的吧. towhee本身不是計(jì)算成核過程的. 應(yīng)該需要編寫代碼實(shí)現(xiàn). 不過MC 如何計(jì)算時(shí)間呢? 如果用MC步數(shù)是不是會有問題? 是不是用MD更合適一點(diǎn)? 5）還有這個(gè)nucleation barrier該怎么來實(shí)現(xiàn)？ 這個(gè)應(yīng)該是進(jìn)行不同溫度點(diǎn)的模擬, 作速率-溫度關(guān)系圖, 然后通過擬合熱力學(xué)公式得到吧. 6）towhee不能并行，那么其實(shí)一個(gè)任務(wù)就是在一個(gè)process上處理的；如果我的個(gè)人電腦比服務(wù)器的主頻高，那么是不是用臺式機(jī)比服務(wù)器更快呢（因?yàn)榕_式機(jī)是4核的，它是不是能自己把towhee的任務(wù)來用4核一起處理？）？ towhee的不能用多個(gè)核并行計(jì)算一個(gè)任務(wù), 而事實(shí)上也不需要. 一般來說, 模擬多個(gè)任務(wù)可以等價(jià)于模擬一個(gè)大任務(wù), MC的本質(zhì)決定它是這樣的性質(zhì), MD可不行. 所以盡量還是使用towhee提供的并行模式. 注意選擇不同的隨機(jī)種子. 7）希望您能給點(diǎn)關(guān)于towhee代碼學(xué)習(xí)的經(jīng)驗(yàn)，就是對于我要做的這個(gè)方向的。 關(guān)于代碼, 其實(shí)我也沒有太多經(jīng)驗(yàn). 我主要是對結(jié)果進(jìn)行后處理. towhee 本身是一個(gè)general Monte Carlo code, 不是專門為某些體系設(shè)計(jì)的, 但是對于你的體系應(yīng)該是可以勝任的. 主要針對體系的代碼, 差別也主要是建模和分析結(jié)果兩個(gè)過程, 核心代碼都是類似的. 我個(gè)人認(rèn)為如果時(shí)間有限, 不要陷到源代碼中去.當(dāng)然能夠多了解一些代碼對于軟件的使用也是非常有好處的. towhee 的代碼由于需要我簡單的讀過一些, 還是比較清楚的. 可以嘗試從總體上去學(xué)習(xí)下. 對于結(jié)果的后處理的代碼, 建議從修改towhee 提供的代碼開始. 使用語言以自己最熟悉的為準(zhǔn), 如果都湊合, 推薦 fortran90. 多學(xué)習(xí) linux bash 腳本也是很有好處的, 我現(xiàn)在簡單的數(shù)據(jù)處理基本都用它. 最后強(qiáng)調(diào), 不要陷到代碼中去, 能夠完成自己的用途就可以了, 越簡單越直接越好. 多多思考跟你體系相關(guān)的問題, 多查文獻(xiàn). ///希望能多多向你學(xué)習(xí)towhee知識，導(dǎo)師不管，自己一個(gè)人摸索。擔(dān)心畢業(yè)十大問題啊，謝謝了 大家基本都差不多, 堅(jiān)定信心, 可以做的很好. 其實(shí)模擬計(jì)算這個(gè)方向比起做實(shí)驗(yàn)還是比較透明的, 重復(fù)文獻(xiàn)一般沒有問題, 技術(shù)問題一般論壇上都會回答. 我也是在學(xué)習(xí)的過程, 大家多多討論, 相互學(xué)習(xí). [ Last edited by 043114076 on 2012-4-26 at 09:55 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 towhee精華 | towhee的運(yùn)用討論��！ |
» 本帖已獲得的紅花（最新10朵）
zhang_jaj» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有13人回復(fù)
-
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
309求調(diào)劑
已經(jīng)有17人回復(fù)
-
省雙一流重點(diǎn)一本大學(xué)招收調(diào)劑
已經(jīng)有5人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有7人回復(fù)
-
267求調(diào)劑
已經(jīng)有4人回復(fù)
-
359求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑
已經(jīng)有4人回復(fù)
-
08工科275分求調(diào)劑
已經(jīng)有4人回復(fù)
-
309分085801求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 關(guān)于分子模塊求助的一些建議
已經(jīng)有29人回復(fù)
- towhee軟件的安裝問題
已經(jīng)有4人回復(fù)
- Towhee能模擬氣體吸附等溫線嗎？
已經(jīng)有10人回復(fù)
- towhee編譯失敗
已經(jīng)有9人回復(fù)
- 運(yùn)行 TOWHEE 時(shí)出現(xiàn)的問題 ，希望有人幫我解決一下，非常感謝
已經(jīng)有8人回復(fù)
- 使用amber力場，一些參數(shù)無法獲得
已經(jīng)有5人回復(fù)
- 想做粒子和固體（金屬）相互作用，不知道lammps能否實(shí)現(xiàn)？
已經(jīng)有6人回復(fù)
- 求問分子模擬斷檢問題。。。
已經(jīng)有5人回復(fù)
- LAMMPS并行安裝問題
已經(jīng)有4人回復(fù)
- MS 進(jìn)行溶液分子動力學(xué)計(jì)算時(shí)在哪里設(shè)置力場參數(shù)
已經(jīng)有5人回復(fù)
- 關(guān)于氣體液化的問題~~~~~
已經(jīng)有13人回復(fù)
- 求用到AMBER force field的開源軟件
已經(jīng)有10人回復(fù)
- 水分子在外場下鍵角、鍵長的變化？
已經(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)于氣液平衡的簡單問題
已經(jīng)有24人回復(fù)
- 石墨烯 lammps
已經(jīng)有11人回復(fù)
- 石墨烯和下層的SiO2會成鍵嗎？
已經(jīng)有7人回復(fù)
- 使用towhee進(jìn)行GEMC模擬相平衡的MC運(yùn)動定義問題..
已經(jīng)有15人回復(fù)
- 【討論】關(guān)于爆炸鹽/浴鹽/沐浴球的配方
已經(jīng)有4人回復(fù)
- 【求助】弱問有關(guān)AMBER遠(yuǎn)程調(diào)用的問題？
已經(jīng)有3人回復(fù)
- 【求助】分子的LJ勢能參數(shù)的擬合
已經(jīng)有8人回復(fù)
41樓2012-09-17 15:22:48
2樓2012-04-26 19:15:47
|
你好，我是按照你之前一個(gè)帖子直接編譯的并行towhee-7版的，但是這個(gè)里面沒有那個(gè)towheebase的目錄（應(yīng)該是6版以前才有的把）。在做練習(xí)的時(shí)候，我把towhee_input文件的forcefield改成了自己的絕對路徑 /public/software/towhee/Forcefield/() 然后towhee_parallel如下 1 towhee.log /public/software/towhee/Examples/Canonical/Amber_Villin 之后執(zhí)行了mpirun -np 4 towhee 顯示出來的為 Number ofjobs:0 Simulations complete Total time:0.000059 seconds 而且沒有towhee.log文件 還有官網(wǎng)說安裝需要的autoconf與automake文件是不是可以不用？ example中的makefile文件用來編譯什么的？ 非常感謝！（因?yàn)橐恍┰驅(qū)嶒?yàn)室不能上網(wǎng)，只能晚上回來求助） [ 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ā)表 | |
|---|---|---|---|---|
