slow_winds

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
帖子: 9
在線: 40.7小時(shí)
蟲號: 2865052
注冊: 2013-12-11

[交流] VASP計(jì)算壓電張量不對稱

請教各位一下，在下剛剛轉(zhuǎn)來用VASP做鈣鈦礦結(jié)構(gòu)壓電材料的第一性原理計(jì)算，但是現(xiàn)在計(jì)算中出現(xiàn)了一些問題，我計(jì)算P4mm的四方相結(jié)構(gòu)時(shí)，壓電離子貢獻(xiàn)24與15方向不相等，導(dǎo)致計(jì)算出來的e24也不等于e15，請教下是什么原因。
POSCAR
KNbO3
1.00000000000000
   3.9780718153912340 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.9780718153912340 0.0000000000000000
   0.0000000000000000 0.0000000000000000 4.0697524641330913
K Nb O
1 1 3
Direct
  0.0000000000000000  0.0000000000000000  0.0208122828936224
  0.5000000000000000  0.5000000000000000  0.5370745926342191
  0.5000000000000000  0.5000000000000000  0.9956543127557665
  0.0000000000000000  0.5000000000000000  0.4982294058581956
  0.5000000000000000  0.0000000000000000  0.4982294058581956

INCAR
SYSTEM = KNbO3

  ISTART =  0
  ICHARG =  2
  PREC =  Accurate
  EDIFF  =  1E-06
  EDIFFG =  -0.005
  NSW =  1
  POTIM  =  0.2
  NFREE  =  2
  IBRION =  8
  LEPSILON = Ture
  ISIF =  3
  ISPIN  =  1
  NELM =  60
ISMEAR = 0
SIGMA  = 0.2
LREAL = Auto
LCHARGE = False
LWAVE = False

計(jì)算結(jié)果
PIEZOELECTRIC TENSOR  for field in x, y, z       (C/m^2)
      XX       YY       ZZ       XY       YZ       ZX
  --------------------------------------------------------------------------------
  x    0.00000    0.00000    0.00000    0.00000    0.00000    0.19952
  y    0.00000    0.00000    0.00000    0.00000    0.19957    0.00000
  z    0.16714    0.16714    0.06456    0.00000    0.00000    0.00000

PIEZOELECTRIC TENSOR IONIC CONTR  for field in x, y, z       (C/m^2)
      XX       YY       ZZ       XY       YZ       ZX
  --------------------------------------------------------------------------------
  x -0.00020    0.00055 -0.00176    0.00000    0.00000 -5.47889
  y -0.00076    0.00097 -0.00008    0.00000    1.73507    0.00000
  z -0.46524 -0.46524    4.83952    0.00000    0.00000    0.00000

d33值利用d=e*s計(jì)算出來只有20pC/N多點(diǎn)。煩請各位告知哪里有問題

回復(fù)此樓

» 猜你喜歡

342求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)專業(yè)279調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有20人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有9人回復(fù)
一志愿a區(qū)211，085601-307分求調(diào)劑已經(jīng)有23人回復(fù)
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑已經(jīng)有8人回復(fù)
材料調(diào)劑已經(jīng)有10人回復(fù)
288求調(diào)劑一志愿哈工大材料與化工已經(jīng)有13人回復(fù)
材料科學(xué)與工程339求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)308分調(diào)劑已經(jīng)有11人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2018-02-20 01:55:29

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

tsglss

鐵桿木蟲 (著名寫手)

應(yīng)助: 41 (小學(xué)生)
金幣: 6014.2
紅花: 19
帖子: 1875
在線: 940小時(shí)
蟲號: 321641
注冊: 2007-03-11
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

請問，為什么要e24和e15相等呢？根據(jù)哪個(gè)文獻(xiàn)？

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

做好自己

2樓2018-02-20 15:40:41

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

tsglss

鐵桿木蟲 (著名寫手)

應(yīng)助: 41 (小學(xué)生)
金幣: 6014.2
紅花: 19
帖子: 1875
在線: 940小時(shí)
蟲號: 321641
注冊: 2007-03-11
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

在你的結(jié)果中在哪里能看出e24和e15呢？另外，d33和d等于es這個(gè)公式是什么意思？謝謝

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

做好自己

3樓2018-02-20 21:03:38

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

slow_winds

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
帖子: 9
在線: 40.7小時(shí)
蟲號: 2865052
注冊: 2013-12-11

引用回帖:

2樓: Originally posted by tsglss at 2018-02-20 15:40:41
請問，為什么要e24和e15相等呢？根據(jù)哪個(gè)文獻(xiàn)？

4mm點(diǎn)群的x，y是相同的，15和24分量可以通過對稱變換得出相等。

贊一下

回復(fù)此樓

4樓2018-02-21 01:36:57

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

slow_winds

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
帖子: 9
在線: 40.7小時(shí)
蟲號: 2865052
注冊: 2013-12-11

引用回帖:

3樓: Originally posted by tsglss at 2018-02-20 21:03:38
在你的結(jié)果中在哪里能看出e24和e15呢？另外，d33和d等于es這個(gè)公式是什么意思？謝謝

xZX是15，yYZ是24，離子貢獻(xiàn)上一個(gè)-5.4，一個(gè)1.7.
d=e*s，d和e是3x6矩陣，s是彈性柔順系數(shù)，是彈性勁度系數(shù)c的6x6逆矩陣。

贊一下

回復(fù)此樓

5樓2018-02-21 01:41:08

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

tsglss

鐵桿木蟲 (著名寫手)

應(yīng)助: 41 (小學(xué)生)
金幣: 6014.2
紅花: 19
帖子: 1875
在線: 940小時(shí)
蟲號: 321641
注冊: 2007-03-11
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

感謝耐心回答，你計(jì)算的彈性常數(shù)是否與他人計(jì)算文獻(xiàn)一致呢？若一樣，不應(yīng)該出現(xiàn)問題吧。若不一樣，是不是優(yōu)化的過程中isif的參數(shù)設(shè)置是否合理，換句話說優(yōu)化的時(shí)候結(jié)構(gòu)的對稱性發(fā)生相變了呢。

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

做好自己

6樓2018-02-21 08:13:49

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

slow_winds

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
帖子: 9
在線: 40.7小時(shí)
蟲號: 2865052
注冊: 2013-12-11

引用回帖:

6樓: Originally posted by tsglss at 2018-02-21 08:13:49
感謝耐心回答，你計(jì)算的彈性常數(shù)是否與他人計(jì)算文獻(xiàn)一致呢？若一樣，不應(yīng)該出現(xiàn)問題吧。若不一樣，是不是優(yōu)化的過程中isif的參數(shù)設(shè)置是否合理，換句話說優(yōu)化的時(shí)候結(jié)構(gòu)的對稱性發(fā)生相變了呢。
...

彈性常數(shù)是用IBRION=6另算的，我對比實(shí)驗(yàn)值還比較符合。優(yōu)化是直接ISIF=3，結(jié)構(gòu)沒有變化，POSCAR這個(gè)數(shù)據(jù)就是我優(yōu)化后的四方相P4mm的晶格。

贊一下

回復(fù)此樓

7樓2018-02-21 08:56:00

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

ak81

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12.3
帖子: 15
在線: 40.8小時(shí)
蟲號: 5486965
注冊: 2017-01-15
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

樓主，你好，請問一下，你用vasp計(jì)算壓電時(shí)，是直接用PIEZOELECTRIC TENSOR 做壓電張量矩陣，還是考慮離子PIEZOELECTRIC TENSOR IONIC CONTR這部分，如果考慮離子，那么，你是直接計(jì)算電子與離子部分的和，還是其他？求解答

贊一下

回復(fù)此樓

8樓2018-05-31 11:53:28

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

honoured

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 181.1
散金: 42
紅花: 1
帖子: 86
在線: 49.7小時(shí)
蟲號: 4102301
注冊: 2015-09-25
性別: GG

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

可能是k點(diǎn)小，或者精度低造成的

贊一下

回復(fù)此樓

愛好運(yùn)動(dòng)健身，喜歡讀書，沉迷于三教九流，陶醉于藝術(shù)，致力于學(xué)術(shù)研究

9樓2018-10-06 11:39:44

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

不二3366

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 300.5
帖子: 22
在線: 75.7小時(shí)
蟲號: 6673869
注冊: 2017-06-03
性別: MM
專業(yè): 原子和分子物理

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

7樓: Originally posted by slow_winds at 2018-02-21 08:56:00
彈性常數(shù)是用IBRION=6另算的，我對比實(shí)驗(yàn)值還比較符合。優(yōu)化是直接ISIF=3，結(jié)構(gòu)沒有變化，POSCAR這個(gè)數(shù)據(jù)就是我優(yōu)化后的四方相P4mm的晶格。...

我感覺應(yīng)該是您的模型建立可能有問題，麻煩問一下那您有沒有把理論計(jì)算值和實(shí)驗(yàn)值進(jìn)行對比？

贊一下

回復(fù)此樓

10樓2018-12-25 11:50:19

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

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

[考研] 284求調(diào)劑 +9	小熊～～ 2026-03-31	9/450	2026-03-31 18:22 by 253863592
[考研] 085602 307分求調(diào)劑 +9	不知道叫什么！ 2026-03-26	9/450	2026-03-31 18:10 by 544594351
[考研] 085600，321分求調(diào)劑 +7	大饞小子 2026-03-31	7/350	2026-03-31 18:02 by 無際的草原
[考研] 290分調(diào)劑求助 +10	吉祥止止陳 2026-03-25	10/500	2026-03-31 17:54 by 544594351
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +6	12.27 2026-03-27	6/300	2026-03-31 16:36 by vendetta_GG
[考研] 調(diào)劑310 +13	溫柔的晚安 2026-03-25	14/700	2026-03-31 13:03 by 記事本2026
[考研] 土木304求調(diào)劑 +4	兔突突突， 2026-03-31	5/250	2026-03-31 11:29 by 北風(fēng)之神.
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6	星沉uu 2026-03-26	7/350	2026-03-31 10:19 by GdShizy
[考研] 274求調(diào)劑 +6	xiao愛同學(xué) 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 282求調(diào)劑 +3	呼吸都是減肥 2026-03-31	3/150	2026-03-31 10:00 by wangjy2002
[考研] 一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-30	5/250	2026-03-31 09:41 by zhshch
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4	橙子（勝意） 2026-03-30	4/200	2026-03-31 00:44 by Linzejun
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 285求調(diào)劑 +6	AZMK 2026-03-29	9/450	2026-03-30 21:02 by dophin1985
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 343求調(diào)劑 +6	愛羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無際的草原
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG