spur

金蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
貴賓: 1.21
金幣: 4994.7
散金: 559
紅花: 5
帖子: 866
在線: 374.9小時
蟲號: 215683
注冊: 2006-03-11
性別: GG
專業(yè): 第一性原理計算

[交流] 【求助】如何用vasp5.2中的DFPT方法算Born effective charge？已有5人參與

請教如何用vasp5.2中的DFPT方法算Born effective charge。
我按照vasp手冊上說的設(shè)置了INCAR，但沒在OUTCAR中得到Born effective charge。有算過這個東西的朋友請指點一下。謝謝！
下面是我的INCAR，有不對的地方請指教。
#Startparameter for this run
ISTART=0
ISPIN=2
PREC=Accurate
ENCUT=400
GGA=91
VOSKOWN =1
MAGMOM=6*0  5  -5  2*0
LNABLA=.TURE.

IALGO=48
LDIAG=.TURE.
NPAR=2
LPLANE = .TRUE.

#Electronic relaxation
NELM=80; NELMIN=6;
NELMDL=3
LREAL=AUTO
EDIFF=0.1E-3
ISMEAR=0; SIGMA=0.002

#LDA+U
LDAU=.TURE.
LDAUTYPE= 2
LDAUL= -1  2  -1
LDAUU= 0.00  3.00  0.0
LDAUJ= 0.00  1.00  0.0
LDAUPRINT= 0

LMAXMIX=4
AMIX    = 0.2
BMIX    = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001

#Ionic relaxation
#EDIFFG=-0.02
#NSW=400
IBRION=8
#ISIF=2
ISYM=0
POTIM=0.1
LEPSILON= .TURE.
LRPA= .FALSE.

#Write flag
LWAVE=.F.
LCHARG=.F.
之后我又做了靜態(tài)自洽計算得到CHGCAR，然后設(shè)置ICHARG=11來做，結(jié)果vasp不執(zhí)行這樣的計算。

回復(fù)此樓

1樓 2010-04-14 22:26:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

后天一

木蟲 (小有名氣)

我們的自由！

1ST強帖: 3
應(yīng)助: 9 (幼兒園)
金幣: 2544.6
紅花: 4
帖子: 204
在線: 311小時
蟲號: 519608
注冊: 2008-03-06
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
qasd(金幣+2):謝謝 2010-04-19 23:51

NB: LDAUL, LDAUU, and LDAUJ must be specified for all atomic species!

贊一下(1人)

回復(fù)此樓

2樓2010-04-18 14:58:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

spur

金蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
貴賓: 1.21
金幣: 4994.7
散金: 559
紅花: 5
帖子: 866
在線: 374.9小時
蟲號: 215683
注冊: 2006-03-11
性別: GG
專業(yè): 第一性原理計算

引用回帖:

Originally posted by 后天一 at 2010-04-18 14:58:43:
NB: LDAUL, LDAUU, and LDAUJ must be specified for all atomic species!

Thank you for your response.
That is not the point. There are three elements in my calculation, I believe the LDA+U parameters are corresponding to each element, respectively.
Actually, I found the answer several days ago. It is because of the wrong compilation of vasp5.2.

贊一下

回復(fù)此樓

3樓2010-04-19 22:25:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dhc198611

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2488.4
散金: 20
紅花: 2
帖子: 280
在線: 573.9小時
蟲號: 780288
注冊: 2009-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

ICHARG= 1

LBERRY= .T.
IGPAR=3
NPPSTR = 8
DIPOL = 0.15 0.15 0.15

這樣設(shè)置應(yīng)該可以的，我試過的！

贊一下

回復(fù)此樓

科研無止境

4樓2010-04-20 08:20:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

spur

金蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
貴賓: 1.21
金幣: 4994.7
散金: 559
紅花: 5
帖子: 866
在線: 374.9小時
蟲號: 215683
注冊: 2006-03-11
性別: GG
專業(yè): 第一性原理計算

★
xiaohunhun(金幣+1):xiexie 2010-04-20 12:14

引用回帖:

Originally posted by dhc198611 at 2010-04-20 08:20:00:
ICHARG= 1

LBERRY= .T.
IGPAR=3
NPPSTR = 8
DIPOL = 0.15 0.15 0.15

這樣設(shè)置應(yīng)該可以的，我試過的！

這個是Berry phase的計算方法，跟linear response不同，呵呵。Berry phase 用vasp4.6就可以做，而linear response方法只有5.2以上的版本才行。

贊一下(1人)

回復(fù)此樓

5樓2010-04-20 09:48:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunray55

銀蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 370.8
散金: 92
紅花: 5
帖子: 226
在線: 109.4小時
蟲號: 596846
注冊: 2008-09-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

版主能否詳細(xì)講一下你的經(jīng)驗（比如，編譯要注意什么，incar如何設(shè)置，結(jié)果如何分析等等），對我們這些也要計算鐵電材料的人是個很好的教材。非常感謝。
宏劍那個教材太簡單了。

引用回帖:

Originally posted by spur at 2010-04-19 22:25:43:

Thank you for your response.
That is not the point. There are three elements in my calculation, I believe the LDA+U parameters are corresponding to each element, respectively.
Actually, I foun ...

贊一下

回復(fù)此樓

6樓2010-06-17 16:18:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ydeng2639

木蟲 (小有名氣)

半導(dǎo)體納米材料

應(yīng)助: 3 (幼兒園)
金幣: 1687.6
帖子: 203
在線: 241.5小時
蟲號: 721776
注冊: 2009-03-13
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

我能問一下用Berry Phase方法怎么得到自發(fā)極化強度單位C/m2？

贊一下

回復(fù)此樓

不做專家對不起自己

7樓2010-06-18 15:09:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunray55

銀蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 370.8
散金: 92
紅花: 5
帖子: 226
在線: 109.4小時
蟲號: 596846
注冊: 2008-09-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
cenwanglai(金幣+1):謝謝回復(fù) 2010-06-21 16:16:36

VASP的結(jié)果是e/A,除以體積后，換算一下就得到C/m2了
e~1.6XXC
A~10-10m

引用回帖:

Originally posted by ydeng2639 at 2010-06-18 15:09:34:
我能問一下用Berry Phase方法怎么得到自發(fā)極化強度單位C/m2？

贊一下(1人)

回復(fù)此樓

8樓2010-06-21 16:15:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lizhentc

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2625.6
帖子: 93
在線: 22.9小時
蟲號: 923748
注冊: 2009-12-09
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

????

回復(fù)此樓

Where there is a will,there is a way.Take time while time is,for time will be away.

9樓2010-12-06 09:39:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

隨便注冊

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 445.2
紅花: 2
帖子: 48
在線: 99.8小時
蟲號: 2000882
注冊: 2012-09-14
專業(yè): 光學(xué)

好高端啊~~~

回復(fù)此樓

10樓2013-04-02 15:23:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 化學(xué)工程321分求調(diào)劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領(lǐng)袖
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學(xué)員8dgXkO
[考研] 298求調(diào)劑 +4	上岸6666@ 2026-03-20	4/200	2026-03-21 17:14 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan