 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1754 | 回復(fù): 4 | |||
那年的小艾金蟲 (初入文壇)
|
[求助]
關(guān)于態(tài)密度計算的一些問題
|
|
本人先做了晶體的結(jié)構(gòu)優(yōu)化后，將INCAR CONTCAR POTCAR KPOINTS考入新建的DOS文件夾中并將CONTCAR改成了POSCAR，準備進行態(tài)密度的計算。請問INCAR中的設(shè)置哪些需要修改的（體系中有Eu,要考慮自旋） 新手一枚，望各位多加指點。 |
回復(fù)此樓» 猜你喜歡
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有6人回復(fù)
-
318求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有6人回復(fù)
-
一志愿武理材料工程348求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有6人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有7人回復(fù)
-
303求調(diào)劑
已經(jīng)有6人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 計算態(tài)密度的時候設(shè)置的NPAR是什么意思
已經(jīng)有17人回復(fù)
- vasp計算態(tài)密度結(jié)果出問題
已經(jīng)有11人回復(fù)
- VASP 計算DOS，態(tài)密度中費米能的作用
已經(jīng)有3人回復(fù)
- 能帶結(jié)構(gòu)和態(tài)密度問題
已經(jīng)有19人回復(fù)
- 怎么根據(jù)態(tài)密度圖積分算區(qū)間的電子數(shù)
已經(jīng)有8人回復(fù)
- 關(guān)于太陽能電池態(tài)密度和載流子復(fù)合問題
已經(jīng)有8人回復(fù)
- 用vasp如何計算能帶圖和態(tài)密度圖呢？
已經(jīng)有11人回復(fù)
- VASP計算態(tài)密度出現(xiàn)問題
已經(jīng)有12人回復(fù)
- abinit能否處理一個位置被兩種元素占據(jù)的固溶體的能帶和態(tài)密度
已經(jīng)有10人回復(fù)
- 【求助】關(guān)于MS軟件中能帶和態(tài)密度的計算問題
已經(jīng)有8人回復(fù)
- 【求助】計算AlN的態(tài)密度請教
已經(jīng)有4人回復(fù)
- 用MS的castep模塊計算能帶結(jié)構(gòu)和態(tài)密度，為什么很長時間計算還不結(jié)束
已經(jīng)有6人回復(fù)
- 【求助】呼叫zhangguangping...關(guān)于態(tài)密度的計算
已經(jīng)有11人回復(fù)
- 【求助】MS中局域態(tài)密度的計算
已經(jīng)有3人回復(fù)
- 【求助】siesta 計算DOS應(yīng)該注意什么問題
已經(jīng)有3人回復(fù)
- 【求助】關(guān)于vasp態(tài)密度計算
已經(jīng)有10人回復(fù)
- 【求助】關(guān)于態(tài)密度的計算的基礎(chǔ)問題
已經(jīng)有10人回復(fù)
- 【求助】VASP計算電子態(tài)密度
已經(jīng)有34人回復(fù)
- 【求助】MS５.０ Castep 計算不出部分態(tài)密度（PDOS）怎么辦？（站內(nèi)有人提問MS4.3，?/a>
已經(jīng)有19人回復(fù)
- 【求助】怎樣由能帶結(jié)構(gòu)和態(tài)密度看計算的結(jié)構(gòu)是不是可能存在超導(dǎo)現(xiàn)象
已經(jīng)有12人回復(fù)
- 【求助】態(tài)密度計算究竟是自洽計算，還是非自洽計算？
已經(jīng)有24人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
identation: 金幣+5, 謝謝指教 2013-10-28 17:43:26
那年的小艾: 金幣+30, ★★★很有幫助, 非常感謝，很詳細的解答 2013-10-28 20:33:50
感謝參與，應(yīng)助指數(shù) +1
identation: 金幣+5, 謝謝指教 2013-10-28 17:43:26
那年的小艾: 金幣+30, ★★★很有幫助, 非常感謝，很詳細的解答 2013-10-28 20:33:50
|
VASP中電子態(tài)密度計算的流程 主要分成三步：一、結(jié)構(gòu)優(yōu)化；二、靜態(tài)自洽計算；三、非自洽計算 以Al－FCC為例子 第一步 結(jié)構(gòu)優(yōu)化 輸入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步 靜態(tài)自洽計算 INCAR： PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5 輸入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步計算是在結(jié)構(gòu)優(yōu)化的結(jié)果上進行的，所以開始第二步的時候，將第一步中的輸入文件INCAR, POTCAR, POSCAR, KPOINT 以及 C* 文件放入靜態(tài)自洽計算中去，并且將CONTCAR 拷貝到 POSCAR中，然后運行VASP。計算結(jié)果中的Fermi能是準確的，需要記下。 第三步 非自洽計算 INCAR： PREC = Medium，ICHARG = 11，ISMEAR = -5，LORBIT = 10 或者 11（這時可不設(shè)RWIGS)，ISTART = 1 在第二步自洽計算的基礎(chǔ)上進行，修改輸入文件INCAR, POTCAR, POSCAR, KPOINT。 INCAR文件 System=Al ISTART=1 ISMEAR=-5 SIGMA=0.2 ICHARG=11 RWIGS=1.402 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 1.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 21 21 21 0.0 0.0 0.0 VASP中文手冊：在進行能帶和DOS計算時，ISMEAR 不能使用N階MP方法。因為MP方法在空軌道上有負的占據(jù)，所以求得的能帶和DOS是不正確的。但是從其它地方看到---“提示：在計算能帶結(jié)構(gòu)時，采ISMEAR = 0或1對結(jié)果的影響非常小，可以認為是一樣的。但是不能采用ISMEAR = -5 或-4。”ISMEAR到底多少？ 計算能帶：ICHARG = 11 導(dǎo)體的話，用ISMEAR=1； 半導(dǎo)體或絕緣體，用ISMEAR=0 。 計算 DOS： ICHARG = 11 ISMEAR = -5 計算的時候，金屬可用0、1，非金屬不要大過0，體材料可用－4、－5（面的話就用－1、0吧？） 參考http://210.45.66.120/doc/vasp5.0/node224.html 設(shè)置完成后進行計算，計算完后，得到包含了態(tài)密度值的DOSCAR文件，采用split_dos對態(tài)密度文件DOSCAR進行分割，得到總態(tài)密度DOS0，各個原子的分波態(tài)密度DOS1，DOS2……。另外在運行split_dos程序?qū)OSCAR文件分割時，要保證當前目錄下有對應(yīng)的OUTCAR和POSCAR文件。 分割后的DOS0，DOS1…等文件的能量值是以費米能級作為能量參考零點。DOS0的第一列數(shù)據(jù)是能量值，單位為eV；第二列數(shù)據(jù)是總態(tài)密度的值，單位 State/eV.unit cell；第三列數(shù)據(jù)是總態(tài)密度的積分值，也就是電子數(shù)，單位為electrons。DOS1是第一個原子的分波態(tài)密度值，其中的第一列數(shù)據(jù)是能量值，單位為eV；第二、三、四列數(shù)據(jù)分別對應(yīng)于s、p、d態(tài)的分波態(tài)密度值，單位為State/eV.atom。其他的DOS文件與DOS1類似。 |
2樓2013-10-28 16:20:37
3樓2013-11-13 12:46:03
4樓2013-11-13 15:17:36
5樓2013-11-13 17:01:35
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 那年的小艾 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|