 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1536 | 回復(fù): 1 | |||
[交流]
chgdiff.pl 處理電荷差分密度 已有1人參與
|
|
本來(lái)想用chgdiff.pl處理電荷差分密度, 但處理三個(gè)文件時(shí)有個(gè)小bug的, 浪費(fèi)半天時(shí)間 OMG 寫下來(lái)供借鑒 1. 因?yàn)樵有畔⒈Ａ袅薴ile1的信息, 應(yīng)該倒過(guò)來(lái), file1-file2 2. printf OUT輸出前加個(gè)空格, 保持與VASP格式一致 修改: 81c81 < $line1[$j] = $line1[$j]-$line2[$j]; --- > $line1[$j] = $line2[$j]-$line1[$j]; 83c83 < printf OUT " %18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; --- > printf OUT "%18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; chgdiff.pl CHGCAR1 CHGCAR2 mv CHGCAR_diff tmp chgdiff.pl tmp CHGCAR3 得到差分密度:CHGCAR1-CHGCAR2-CHGCAR3 用網(wǎng)上找的fortran處理得到plane-averaged difference electronic density, 注意刪除原子符號(hào)那行 齊活! |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
朱全喜» 猜你喜歡
申博26年
已經(jīng)有4人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿C9材料與化工專業(yè)總分300求調(diào)劑
已經(jīng)有5人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
299求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有4人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有3人回復(fù)
-
材料學(xué)碩，求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩331求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 急需會(huì)對(duì)chgdiff.pl 腳本對(duì)電荷密度做差的高手幫忙！
已經(jīng)有4人回復(fù)
- 差分電荷密度圖的疑惑
已經(jīng)有41人回復(fù)
- CHGCAR的文件含義
已經(jīng)有7人回復(fù)
- P4vasp畫PDOS，請(qǐng)教一下，謝謝。
已經(jīng)有14人回復(fù)
- 如何用xcrysden畫電荷密度圖
已經(jīng)有8人回復(fù)
- vasp如何計(jì)算差分電荷密度
已經(jīng)有7人回復(fù)
- 差分電荷密度的計(jì)算
已經(jīng)有11人回復(fù)
朱全喜
新蟲 (小有名氣)
- 應(yīng)助: 19 (小學(xué)生)
- 金幣: 2977.9
- 紅花: 12
- 帖子: 205
- 在線: 132.6小時(shí)
- 蟲號(hào): 1972785
- 注冊(cè): 2012-09-04
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2015-05-27 16:28:40
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 shixha 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 085602 289分求調(diào)劑
+3
|
WWW西西弗斯 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 材料292調(diào)劑 +6 | 橘頌思美人 2026-03-23 | 6/300 |
|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|