 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1537 | 回復(fù): 1 | |||
[交流]
chgdiff.pl 處理電荷差分密度 已有1人參與
|
|
本來想用chgdiff.pl處理電荷差分密度, 但處理三個文件時有個小bug的, 浪費半天時間 OMG 寫下來供借鑒 1. 因為原子信息保留了file1的信息, 應(yīng)該倒過來, file1-file2 2. printf OUT輸出前加個空格, 保持與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, 注意刪除原子符號那行 齊活! |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
朱全喜» 猜你喜歡
材料學(xué)碩333求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有7人回復(fù)
-
收08調(diào)劑生
已經(jīng)有6人回復(fù)
-
材料專碩找調(diào)劑
已經(jīng)有5人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有5人回復(fù)
-
化工專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
274求調(diào)劑
已經(jīng)有5人回復(fù)
-
調(diào)劑
已經(jīng)有3人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有15人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有16人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 急需會對chgdiff.pl 腳本對電荷密度做差的高手幫忙！
已經(jīng)有4人回復(fù)
- 差分電荷密度圖的疑惑
已經(jīng)有41人回復(fù)
- CHGCAR的文件含義
已經(jīng)有7人回復(fù)
- P4vasp畫PDOS，請教一下，謝謝。
已經(jīng)有14人回復(fù)
- 如何用xcrysden畫電荷密度圖
已經(jīng)有8人回復(fù)
- vasp如何計算差分電荷密度
已經(jīng)有7人回復(fù)
- 差分電荷密度的計算
已經(jīng)有11人回復(fù)
朱全喜
新蟲 (小有名氣)
- 應(yīng)助: 19 (小學(xué)生)
- 金幣: 2977.9
- 紅花: 12
- 帖子: 205
- 在線: 132.6小時
- 蟲號: 1972785
- 注冊: 2012-09-04
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
2樓2015-05-27 16:28:40
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 shixha 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-24 | 3/150 |
|
|
|
[考研] 340求調(diào)劑 +3 | 話梅糖111 2026-03-24 | 3/150 |
|
|
|
[考研]
 0854 考研調(diào)劑 招生了！AI 方向
+5
|
pk3725069 2026-03-19 | 17/850 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|