| 9 | 1/1 | 返回列表 |
| 查看: 6070 | 回復(fù): 8 | ||||||
nkleof銀蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】表面靜電勢(shì)怎么計(jì)算?
|
看到有的文獻(xiàn)中做出了從表面深入到體材料的各層的靜電勢(shì),如下圖所示:![]() 請(qǐng)問(wèn)這個(gè)靜電勢(shì)數(shù)據(jù)從vasp計(jì)算的結(jié)果中如何提。 |
第一性原理和電化學(xué) | 第一性原理 | 不錯(cuò)的帖子 |
木蟲(chóng) (小有名氣)
銀蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (職業(yè)作家)
|
網(wǎng)上可以找到相關(guān)的工具,vasp組弄的workshop的材料里面好像就有類(lèi)似的工具,臺(tái)灣學(xué)術(shù)圈辦的workshop的材料里面也有處理LOCPOT的工具。自己寫(xiě)一個(gè)也很簡(jiǎn)單。 讀入LOCPOT,然后寫(xiě)一個(gè)求和的子過(guò)程就可以得到了。 pot (i,jk) 這是local potential的數(shù)組,它的維數(shù)在LOPOT文件里面有給出,也一般是NGX, NGY, NGZ。假設(shè)表面或界面的法向方向是在z方向上,那么 do k =1, NGZ summ = 0.0 do j= 1, NGY do i = 1, NGX summ = pot(i,j,k) +summ end do end do summ = 表面或界面的面積 * summ/ (NGX * NGY) write(*,*) k*原胞c軸的長(zhǎng)度, summ end do 要注意的就是各個(gè)量的單位。 |
銅蟲(chóng) (小有名氣)
|
看帖http://www.gaoyang168.com/bbs/viewthread.php?tid=2357737&fpage=1 復(fù)制過(guò)來(lái)的時(shí)候記。 WRITE(*,*) 'Error opening or reading file LOCPOT.' 之前少了個(gè)空格 |
銅蟲(chóng) (小有名氣)
銀蟲(chóng) (正式寫(xiě)手)
銅蟲(chóng) (小有名氣)
銅蟲(chóng) (小有名氣)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|