版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

spur

金蟲 (正式寫手)

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

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

請(qǐng)教如何用vasp5.2中的DFPT方法算Born effective charge。
我按照vasp手冊(cè)上說(shuō)的設(shè)置了INCAR，但沒(méi)在OUTCAR中得到Born effective charge。有算過(guò)這個(gè)東西的朋友請(qǐng)指點(diǎn)一下。謝謝！
下面是我的INCAR，有不對(duì)的地方請(qǐng)指教。
#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)自洽計(jì)算得到CHGCAR，然后設(shè)置ICHARG=11來(lái)做，結(jié)果vasp不執(zhí)行這樣的計(jì)算。

回復(fù)此樓

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

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

后天一

木蟲 (小有名氣)

我們的自由！

1ST強(qiáng)帖: 3
應(yīng)助: 9 (幼兒園)
金幣: 2544.6
紅花: 4
帖子: 204
在線: 311小時(shí)
蟲號(hào): 519608
注冊(cè): 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強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
貴賓: 1.21
金幣: 4994.7
散金: 559
紅花: 5
帖子: 866
在線: 374.9小時(shí)
蟲號(hào): 215683
注冊(cè): 2006-03-11
性別: GG
專業(yè): 第一性原理計(jì)算

引用回帖:

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小時(shí)
蟲號(hào): 780288
注冊(cè): 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)該可以的，我試過(guò)的！

贊一下

回復(fù)此樓

科研無(wú)止境

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

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

spur

金蟲 (正式寫手)

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

★
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)該可以的，我試過(guò)的！

這個(gè)是Berry phase的計(jì)算方法，跟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小時(shí)
蟲號(hào): 596846
注冊(cè): 2008-09-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

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小時(shí)
蟲號(hào): 721776
注冊(cè): 2009-03-13
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

贊一下

回復(fù)此樓

不做專家對(duì)不起自己

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小時(shí)
蟲號(hào): 596846
注冊(cè): 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:
我能問(wèn)一下用Berry Phase方法怎么得到自發(fā)極化強(qiáng)度單位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小時(shí)
蟲號(hào): 923748
注冊(cè): 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的回帖

隨便注冊(cè)

銀蟲 (初入文壇)

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

好高端啊~~~

回復(fù)此樓

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 spur 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 311求調(diào)劑 +13	冬十三 2026-03-15	14/700	2026-03-21 22:10 by peike
[考研] 354求調(diào)劑 +6	Tyoumou 2026-03-18	9/450	2026-03-21 20:47 by lbsjt
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 086500 325 求調(diào)劑 +3	領(lǐng)帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu

24小時(shí)熱門版塊排行榜

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

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: