 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1612 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
dx0620木蟲 (正式寫手)
|
[求助]
用dos_procar處理dos數(shù)據(jù)不準(zhǔn) 已有1人參與
|
|
|
我的vasp版本是5.3，用的dos_procar處理的PROCAR，想分別得到s，p，d的態(tài)密度，但是最后得到的態(tài)密度曲線不對，總態(tài)密度跟s，p，d的對不上，如下圖 我的dos_procar.f90的代碼如下 !=========================================================================== ! INPUT FILES : OUTCAR , PROCAR !=========================================================================== implicit real*8(a-h,o-z) character*64 tmp dimension ddosP(4),ddosT(4) allocatable wt(  ,oc(:,:,:,,eigen(:,,ee(,na(allocatable gdosP(:, ,gdosT(:,pi=4.0d0*atan(1.0d0) ; netot = 1001 allocate(ee(netot),gdosP(4,netot),gdosT(4,netot)) open(1,file='OUTCAR') !================================================================= nline = 0 ; open(99,status='scratch') 111 read(1,'(a64)',end=222) tmp if(tmp(2:8) .eq. 'E-fermi') nline=nline+1 if(tmp(2:8) .eq. 'E-fermi') write(99,'(a64)') tmp go to 111 222 rewind(99) do i=1,nline ; read(99,*) tmp,tmp,fermi ; end do close(99) !================================================================= write(6,'("fermi energy = ",f10.5)') fermi write(6,'("enter the gaussian "  ') ; read(5,*) gaussianwrite(6,'("ISPIN(1/2) = ? " ') ; read(5,*) ispinwrite(6,'("How many atoms are integrated? " ') ; read(5,*) Nnaallocate(na(Nna)) write(6,*) ' atom number :' ; read(5,*) na open( 7,file='PROCAR') if(ispin .eq. 1) then open(11,file='gdosP-s.dat') ; open(51,file='gdosT-s.dat') open(12,file='gdosP-p.dat') ; open(52,file='gdosT-p.dat') open(13,file='gdosP-d.dat') ; open(53,file='gdosT-d.dat') open(14,file='gdosP-A.dat') ; open(54,file='gdosT-A.dat') else if(ispin .eq. 2) then open(11,file='gdosP-up-s.dat') ; open(51,file='gdosT-up-s.dat') open(12,file='gdosP-up-p.dat') ; open(52,file='gdosT-up-p.dat') open(13,file='gdosP-up-d.dat') ; open(53,file='gdosT-up-d.dat') open(14,file='gdosP-up-A.dat') ; open(54,file='gdosT-up-A.dat') open(21,file='gdosP-dn-s.dat') ; open(61,file='gdosT-dn-s.dat') open(22,file='gdosP-dn-p.dat') ; open(62,file='gdosT-dn-p.dat') open(23,file='gdosP-dn-d.dat') ; open(63,file='gdosT-dn-d.dat') open(24,file='gdosP-dn-A.dat') ; open(64,file='gdosT-dn-A.dat') else write(6,*) 'Input ERROR, you must input 1 or 2 ' ; stop end if !=========================================================================== do is=1,ispin if(is .eq. 1) read(7,'(a64)') tmp read(7,'(16x,i3,20x,i4,19x,i4)') nk,nband,nion if(is .eq. 1) then write(6,'("number of kpoints: ",i4)') nk write(6,'("number of bands: ",i4)') nband write(6,'("number of ions: ",i4)') nion end if emin = 9999.0 ; emax= -9999.0 allocate(wt(nk),eigen(nk,nband),oc(nk,nband,nion+1,4)) do k=1,nk read(7,'(a64)') tmp read(7,'(10x,i3,5x,3f11.8,13x,f11.8)') kp,pt1,pt2,pt3,wt(k) read(7,'(a64)') tmp do nb=1,nband read(7,'(8x,9x,f14.8,7x,f12.8)') eigen(k,nb),occ eigen(k,nb) = eigen(k,nb) - fermi if(eigen(k,nb) .gt. emax) emax=eigen(k,nb) if(eigen(k,nb) .lt. emin) emin=eigen(k,nb) read(7,'(a64)') tmp ; read(7,'(a64)') tmp niont = nion +1 if(nion .eq. 1) niont = 1 do ion = 1,niont read(7,'(3x,4f7.3)') (oc(k,nb,ion,j),j=1,4) end do read(7,'(a64)') tmp end do end do ! ############################################################## weight = 0.0 ; gdosP=0.0d0 ; gdosT=0.0d0 do k=1,nk ; weight = weight + wt(k) ; end do do k=1,nk ; wt(k) = wt(k) / weight ; end do estart_ev = int(emin -5.0) eend_ev = int(emax +5.0) de_ev = (eend_ev - estart_ev)/(netot-1) do ne = 1,netot ; ee(ne) = estart_ev + (ne-1) * de_ev ; end do ascal = 1.0/(gaussian*sqrt(pi)) ! ############################################################### do k=1,nk ; do nb=1,nband do i=1,4 ; ddosT(i) = oc(k,nb,niont,i)*wt(k) ; end do do i=1,4 ; ddosP(i) = sum(oc(k,nb,na( ,i)*wt(k)) ; end dodo ne=1,netot dij = ( eigen(k,nb)-ee(ne) )**2 / (gaussian**2) do i=1,4 gdosT(i,ne) = gdosT(i,ne) + ascal*ddosT(i)*exp(-dij) gdosP(i,ne) = gdosP(i,ne) + ascal*ddosP(i)*exp(-dij) end do end do end do ; end do ! ############################################################### if(is .eq. 1) ns=+1 ; if(is .eq. 2) ns=-1 do i=1,4 ; write(is*10+i ,1) (ee(j),ns*gdosP(i,j),j=1,netot) ; end do do i=1,4 ; write(is*10+i+40,1) (ee(j),ns*gdosT(i,j),j=1,netot) ; end do deallocate(wt,eigen,oc) end do !=========================================================================== 1 format(f10.5,f12.4) end ############################################################################ INCAR 如下： ENCUT = 400 ISTART = 1 ICHARG = 11 ISMEAR = 1 LORBIT = 11 #EMIN = -25 #EMAX = 22.5 NBANDS = 64 #SIGMA = 0.1 #EDIFF = 1E-6 #EDIFFG = -1E-3 #NSW = 200 #IBRION = 2 #POTIM = 0.2 #ISIF = 3 #ALGO = Normal #LREAL = Auto PREC = Accurate PSTRESS = 500.00 |
回復(fù)此樓» 猜你喜歡
一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑
已經(jīng)有6人回復(fù)
-
268求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有11人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué)
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
313求調(diào)劑
已經(jīng)有3人回復(fù)
-
南昌大學(xué)材料專碩311分求調(diào)劑
已經(jīng)有6人回復(fù)
-
316求調(diào)劑
已經(jīng)有6人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- dos-procar.f90
已經(jīng)有94人回復(fù)
- 求解split_dos分割文件該怎么理解！？
已經(jīng)有11人回復(fù)
- vasp計(jì)算，DOS圖的分析，分波態(tài)密度的數(shù)據(jù)分析~·
已經(jīng)有4人回復(fù)
- 請教split_dos處理DOS的相關(guān)問題
已經(jīng)有4人回復(fù)
- 【分享】 最新VASP腳本大全 支持VASP 5.3 (有split_dos; VP; Sumdos; sumup_dos 等)
已經(jīng)有615人回復(fù)
- split_dos分割以后，如何找到對應(yīng)的原子？
已經(jīng)有11人回復(fù)
- 【求助成功】split_dos分割DOSCAR的問題
已經(jīng)有14人回復(fù)
- 【求助】那位強(qiáng)人能在c++程序中調(diào)用DOS命令行��！指點(diǎn)一下唄（懸賞給二十金幣）
已經(jīng)有5人回復(fù)
- 【求助成功】怎么使用DOSCAR畫出DOS圖形
已經(jīng)有23人回復(fù)
- 【求助】murn.f gk.f和pbnd.f split_dos和vp程序
已經(jīng)有17人回復(fù)
mywai520
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 339 (大學(xué)生)
- 金幣: 10425.6
- 紅花: 45
- 帖子: 1664
- 在線: 296.9小時(shí)
- 蟲號(hào): 265142
- 注冊: 2006-07-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2014-05-02 16:50:24
mywai520
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 339 (大學(xué)生)
- 金幣: 10425.6
- 紅花: 45
- 帖子: 1664
- 在線: 296.9小時(shí)
- 蟲號(hào): 265142
- 注冊: 2006-07-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2014-05-02 12:24:18
dx0620
木蟲 (正式寫手)
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 4053.6
- 散金: 5
- 紅花: 6
- 帖子: 393
- 在線: 910小時(shí)
- 蟲號(hào): 1303311
- 注冊: 2011-05-22
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-05-02 15:03:31
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +11 | 呼呼呼。。。。 2026-03-17 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 8/400 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|