| 20 | 1/1 | 返回列表 |
| 查看: 5500 | 回復(fù): 19 | |||||
| 【獎勵】 本帖被評價8次,作者nkleof增加金幣 7.5 個 | |||||
[資源]
【分享】處理DOSCAR文件的腳本和程序
|
|||||
|
根據(jù)JMS NRL的split_dos腳本稍做修改,使得能夠處理從s到d所有18個軌道(包括自旋上下)的DOSCAR,生成每個原子的軌道態(tài)密度,然后自己寫了一個C程序,可以把同類原子的p軌道和d的t2g,eg軌道分別求和生成新的軌道態(tài)密度文件。適用于對含過渡金屬元素體系的電子態(tài)密度進(jìn)行處理,可以很方便的得到p軌道和t2g,eg軌道的總態(tài)密度。 使用方法: 僅限linux平臺 1、把三個腳本程序解壓后放到一個/bin目錄并把路徑添加到環(huán)境變量里(或者把三個文件都放到要處理的數(shù)據(jù)所在目錄)。 2、鍵入split_dos.ksh運行分解DOSCAR的腳本,得到體系總的和每個原子的分波態(tài)密度。 3、鍵入sumdos *1 *2 *3執(zhí)行求p以及t2g,eg軌道態(tài)密度和,其中*1, *2, *3為三個參數(shù),依次為要處理的元素符號(如Fe),起始原子序號,最終原子序號(比如Fe原子在POSCAR里面是第5到13號原子,則*2和*3分別為5 13)。 http://d.namipan.com/d/bd4ee2bb9 ... 5c2334a0697fb150000 注意!。。。簊umdos的使用僅限于t2g和eg分別簡并得很好的情況,即它們內(nèi)部不再分裂,因為程序只是簡單的把dxy, dyz, dzx和dz2, dx2分別加起來!如果有較強退簡并情況,如晶體形變造成t2g和eg再分裂,請勿用此程序! [ Last edited by nkleof on 2010-5-24 at 11:14 ] |
納米技術(shù)與能源及模擬 | VASP | 各種各樣的軟件 |
|
我看過他的回帖,實際上就是,假如把dxy, dyz, dzx, dx2, dz2看作5個線性無關(guān)項,那么當(dāng)vasp所采用的坐標(biāo)系不和所定義的x, y, z重合時,輸出的5個軌道可以看做是這5個軌道的線性組合。我不確定vasp是不是會采用一套很特別的坐標(biāo)系,但是我通常建模型都是用內(nèi)坐標(biāo)也就是direct模式,所以得到的結(jié)果應(yīng)該是xyz就是晶體場坐標(biāo)。 我寫的那個sumdos程序只是對特定的幾列進(jìn)行簡單加和,所以只適用于特定情況,否則不保證加出來的t2g和eg正確。我在尖晶石結(jié)構(gòu)中使用過,是適用于內(nèi)坐標(biāo)建模的尖晶石結(jié)構(gòu)的。 因此提醒注意,如果使用那個程序,應(yīng)在使用前看看10和11,12和13,16和17這三組列數(shù)據(jù)(分別對應(yīng)自旋向上和向下)是不是對應(yīng)簡并的t2g三個軌道態(tài)密度,14和15,18和19這兩組列數(shù)據(jù)是不是對應(yīng)eg軌道態(tài)密度。 [ Last edited by nkleof on 2010-5-26 at 17:14 ] |
| 20 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|