| 5 | 1/1 | 返回列表 |
| 查看: 2522 | 回復(fù): 4 | |||
[求助]
離子液體模擬:如何從別的軟件移植力場到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ù)的說明?我得弄清楚各項(xiàng)的含義是什么,網(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ù)啊。 不知我這樣理解是否正確。 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩調(diào)劑 +5 | 椰椰。 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|