| 查看: 4956 | 回復(fù): 18 | ||
future_wl木蟲(chóng) (著名寫(xiě)手)
|
[求助]
VASP計(jì)算介電張量的結(jié)果文件
|
|
計(jì)算完介電張量后(LPEAD=.TRUE LEPSILON=.TRUE), 在結(jié)果文件中搜索DIELECTRIC 發(fā)現(xiàn)兩部分矩陣: HEAD OF MICROSCOPIC STATIC DIELECTRIC TENSOR (independent particle, excluding Hartree and local field effects) ------------------------------------------------------ 7.531256 0.000000 0.000000 0.000000 7.531256 0.000000 0.000000 0.000000 6.922124 ------------------------------------------------------ MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT) ------------------------------------------------------ 7.314593 0.000000 0.000000 0.000000 7.314593 0.000000 0.000000 0.000000 6.607187 ------------------------------------------------------ 請(qǐng)問(wèn)這兩部分分別表示什么?我如果要取介電張量的話應(yīng)該用哪部分? 謝謝大家! |

木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
|
體系微觀介電常數(shù)epsilon(r,r')是一個(gè)二變量的函數(shù),對(duì)其作傅立葉變換后,傅立葉系數(shù)可以寫(xiě)成一個(gè)矩陣的形式epsilon( G, G' ),G和G'是倒格矢。 宏觀介電常數(shù)由微觀介電常數(shù)在原胞內(nèi)積分取均值獲得,可以證明最終結(jié)果等于其傅立葉變換系數(shù)的(0,0)項(xiàng),也就是epsilon(0,0)。 但是一般程序在計(jì)算時(shí)算的不是epsilon(G,G'),而是epsilon^-1(G,G')。這樣必須先取整個(gè)epsilon^-1(G,G')矩陣的逆,再取其(0,0)元,這樣做時(shí)就包含了局域場(chǎng)效應(yīng)。 如果為了計(jì)算簡(jiǎn)便,僅對(duì)epsilon^-1(G,G')的對(duì)角元取逆,就忽略了局域場(chǎng)效應(yīng)。 一般來(lái)說(shuō),局域場(chǎng)效應(yīng)對(duì)inhomogeneous系統(tǒng)影響較大。 |
木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
|
RPA是一種計(jì)算節(jié)點(diǎn)函數(shù)的理論,局域場(chǎng)效應(yīng)是計(jì)算介電函數(shù)過(guò)程中所涉及的一種偏數(shù)學(xué)的近似手段。 這里VASP手冊(cè)里寫(xiě)得我也看得不是太明白,按我理解應(yīng)該是這樣的: 算介電函數(shù)時(shí)先算的是極化函數(shù)X,然后由X計(jì)算eps^-1。由前面DFT計(jì)算得到的波函數(shù)和本征值可以根據(jù)RPA算出X0,更精確的應(yīng)該是由X0算出的X?疾豢紤]RPA的區(qū)別主要在X上。 你的理解應(yīng)該是對(duì)的。 推薦去看下DP和exciting這兩個(gè)TDDFT程序的主頁(yè),上面有講解TDDFT理論的PPT,里面對(duì)極化率節(jié)點(diǎn)函數(shù)講得很清楚。 |
木蟲(chóng) (著名寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|