| 7 | 1/1 | 返回列表 |
| 查看: 2113 | 回復(fù): 6 | ||||
趙紅霞鐵桿木蟲 (著名寫手)
|
[求助]
聚合物的pdb文件生成psf文件用于NAMD動力學(xué)計算 已有1人參與
|
|
想用全原子做聚合物(聚電解質(zhì))在溶液里的構(gòu)象,因為體系大,materials studio算起來慢,想在NAMD里算。在生成psf文件時遇到問題,請各位支招。 1. 先在materials studio構(gòu)建好了聚合物分子,導(dǎo)出為pdb格式。 2. 下一步生成psf文件,不知道哪種方法好用,我暫且先選擇在VMD里借助psfgen轉(zhuǎn)化。pgn文件如下: package require psfgen topology par_all27_prot_lipid_na.inp segment U {pdb PCE1_period30nm11.pdb} coordpdb PCE1_period30nm11.pdb U guesscoord writepdb PCE1_period30nm11.pdb writepsf PCE1_period30nm11.psf 問題來了,因為charmm或者amber力場主要針對蛋白和核酸體系,里面規(guī)定了標準的殘基名,根據(jù)殘基來規(guī)定殘基里每個原子的力場參數(shù)。 我的polymer的殘基名無法識別。請問如何解決? PS:雖然charmm或者amber力場主要針對蛋白和核酸體系,但是我看到文獻里仍然有人拿它們做別的polymer。那么psf如何生成呢? |
鐵桿木蟲 (著名寫手)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
鐵桿木蟲 (著名寫手)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
我不太明白你的意思。lammps里當(dāng)然可以用charmm力場?墒悄愕捏w系里分子的類型如果和charmm力場不匹配,你怎么調(diào)用charmm力場參數(shù)來做模擬? 我不清楚你是怎么生成lammps輸入文件的。如果是用charmm力場,一般通過charmm2lammps.pl這個腳本來生成輸入文件。即便使用這個腳本,你也需要charmm格式的topology和parameter文件,以及psf文件。如果你的pdb文件里原子類型和charmm力場topology文件里的原子類型根本不匹配,是根本不可能生成輸入文件的。 你現(xiàn)在的問題就是,你pdb文件里原子的類型和charmm力場的topology里原子還有分子類型不匹配,所以無法生成psf文件。 首先,你要搞清楚你要用什么力場,如果你要用charmm力場,那就把你的pdb文件里原子類型改成charmm的原子類型。 其次,你要確定charmm的topology文件里是否有你想模擬的分子。如果沒有,你應(yīng)該參照charmm topology文件的格式,自己建一個新的topology文件,文件中應(yīng)該包含你要模擬的高分子的原子質(zhì)量,成鍵,鍵角,二面角等信息。并注意topology文件里的原子名和residue name與你的pdb文件相匹配。至于有沒有腳本可以自動生成topology文件,我不太清楚,請自行查找。 |
鐵桿木蟲 (著名寫手)
新蟲 (初入文壇)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[基金申請]
有必要更換申報口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|