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)動健身，喜歡讀書，沉迷于三教九流，陶醉于藝術(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ā)表

[考研] 生醫(yī)工0831調(diào)劑求推薦 +4	小熊睿睿_s 2026-03-27	6/300	2026-03-31 17:12 by 記事本2026
[考研] 085600 材料與化工 329分求調(diào)劑 +20	Mr. Z 2026-03-25	21/1050	2026-03-31 16:53 by Zzxxxs
[考研] 求收留 +8	1943443204 2026-03-28	8/400	2026-03-31 15:00 by -迷了路啊路
[考研] 調(diào)劑申請 +6	張張張張zy 2026-03-31	6/300	2026-03-31 14:20 by 記事本2026
[考研] 322求調(diào)劑 +10	舊吢 2026-03-24	10/500	2026-03-31 13:31 by 無際的草原
[考研] 272求調(diào)劑，接受跨專業(yè)調(diào)劑！ +3	閑魚盧 2026-03-31	3/150	2026-03-31 13:00 by 替代品000
[考研] 323分食品與營養(yǎng)調(diào)劑 +3	嘿ooo 2026-03-31	3/150	2026-03-31 09:38 by longlotian
[考研] 08工科求調(diào)劑286 +5	tgs_001 2026-03-28	5/250	2026-03-31 08:18 by 一只好果子?
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 285求調(diào)劑 +6	AZMK 2026-03-29	9/450	2026-03-30 21:02 by dophin1985
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +7	YYYYYNNNNN 2026-03-26	9/450	2026-03-30 20:29 by YYYYYNNNNN
[考研] 085601材料工程找調(diào)劑 +17	oatmealR 2026-03-29	18/900	2026-03-30 19:21 by Wang200018
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i