| 5 | 1/1 | 返回列表 |
| 查看: 2527 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
離子液體模擬:如何從別的軟件移植力場到Gromacs(itp文件的含義以及如何寫itp文件)
|
|||
|
我要做純的離子液體盒子,使用的力場是從文獻(xiàn)里下的,是DL-POLY軟件的格式(附件1),我準(zhǔn)備把它寫成gromacs的itp文件,然后我看了下gromacs manual 上關(guān)于itp文件的說明(附件2),里面有一個itp文件的例子,其內(nèi)容和我下的東西看起來差不多,但是覺得里面好像說的都是分子結(jié)構(gòu)啊,而且例子里在top文件里就沒有再給結(jié)構(gòu)了,只是給了分子數(shù)目和include了itp,反而是有這么一行引用力場的文件: #include "gmx.ff/forcefield.itp" 那么manual上給的例子,itp文件是力場還是分子結(jié)構(gòu)? 還有如果我把力場寫好了,還要用高斯優(yōu)化離子對結(jié)構(gòu)做一個pdb,然后做成top再genbox,然后把力場.itp給include進(jìn)去,是這樣嘛?還是像這個例子里寫上分子數(shù)目就行了? 另外,在哪里可以找到itp文件里各個參數(shù)的說明?我得弄清楚各項的含義是什么,網(wǎng)上搜不到...最好有教你怎么寫itp文件的例子 謝謝啦 |
xuexijisuan |
|
1、先讀gromacs的手冊,仔細(xì)閱讀itp文件的那一章,這部分內(nèi)容說的比其他任何人說的都權(quán)威詳細(xì) 2、你可以打開任意一個例子文件的itp文件,用寫字板打開就能看到文件的內(nèi)容和格式。 3、個人愚見:structure 是坐標(biāo)文件,每個原子的xyz坐標(biāo),一般的gro,pdb文件;force field 就是力場了嘛,就是告訴你原子之間的成鍵情況,以及鍵長鍵能大小,非鍵相互作用等等。力場一般你只需要從軟件里選擇就行了,如果你模擬的對象比較特殊,可能需要你自己去尋找新的力場文件。 |
|
比如我的 gro 文件中的一部分: 36DPPC NC3 1 18.208 7.007 59.608 36DPPC PO4 2 18.422 6.811 59.976 36DPPC GL1 3 18.660 7.089 60.226 36DPPC GL2 4 18.587 7.416 60.230 36DPPC C1A 5 18.642 7.077 60.691 36DPPC C2A 6 18.725 7.215 61.110 36DPPC C3A 7 18.652 7.173 61.537 36DPPC C4A 8 18.643 6.878 61.865 36DPPC C1B 9 18.764 7.568 60.577 36DPPC C2B 10 18.676 7.668 60.980 36DPPC C3B 11 18.981 7.644 61.338 36DPPC C4B 12 18.973 7.685 61.800 這是一個DPPC分子的初始坐標(biāo)文件(初始模型),有了這個之后,就需要選擇合適的力場文件來模擬,就需要拓?fù)湮募?.top) 拓?fù)湮募袝嬖V你一共有多少個分子,每個個DPPC分子中,哪兩個原子之間成鍵,哪三個成角等等力場信息。而每個原子的坐標(biāo),隨著模擬的時間都在變化的,但是分子中的各個原子間的拓?fù)潢P(guān)系不會變。這樣不重復(fù)啊。 不知我這樣理解是否正確。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +5 | DLkz1314. 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑 一志愿070303有機化學(xué) +5 | 蘿卜燉青菜 2026-03-28 | 6/300 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +14 | jhybd 2026-03-23 | 19/950 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚 2026-03-27 | 11/550 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|