| 7 | 1/1 | 返回列表 |
| 查看: 2055 | 回復(fù): 6 | ||
[求助]
dump 輸出原子名稱
|
|
各位好: lammps定義原子時用數(shù)字1,2,3,導(dǎo)致在vmd中顯示的name也為1,2,3,請問怎么才能使lammps結(jié)果在vmd顯示的name為真實名字名稱N,H,O? 8660473.png |
金蟲 (著名寫手)
木蟲 (正式寫手)
|
現(xiàn)成的沒有,不過思路簡單: 1、把lammps軌跡載入VMD 2、做循環(huán),直到你的軌跡全部循環(huán)完畢 3、保存成PDB,或者XYZ 格式的文件,其中就有坐標(biāo)信息 去看看VMD所使用的幾個 tcl 命令 for set sel [atomselect top "your atoms that will be saved"] $sel writepdb yourfilename.pdb 或者直接輸出到另一個文件中 set outfile [open $yourfilename.dat w] $sel get {x y z} puts $outfile $sel 思路如此,去查手冊搞懂各個命令,完善此腳本即可。 |
|
昨天正好研究了下這個。你可以看下 dump_modify ?梢允謩有薷牡。dump custom 里面可以選擇輸出元素 element 。 不過不修改的話都是默認(rèn)的 C .(如dump 1 all custom 1 MCM.mol id x y z element )。而 dump_modify 可以修改。如有3種原子 Si O H, 可以寫 dump_modify 1 element Si O H . 不過這個會不會寫到自動生成的輸出格式我就沒試過了。另外補充一點, 自定義輸出的格式定義,實際除了 %d %g 也可以用 %s 就是用來處理個,如對應(yīng)上面那個dump 可以寫 dump_modify 1 format "%5d %9.5g %9.5g %10s"。 如果對文件格式比較熟的話,應(yīng)該可以用這個方法實現(xiàn)。 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | yzysaa 2026-03-21 | 4/200 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|