版塊導(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í)熱門版塊排行榜

返回列表

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 請(qǐng)計(jì)算過加靜水壓下體系的性質(zhì)的高手幫忙看看-----加壓下的壓電系數(shù)

祝各位五一節(jié)快樂！
下面是我計(jì)算壓電系數(shù)的步驟，請(qǐng)各位大俠看一下（用下面的步驟重復(fù)吳志剛的那篇PRL時(shí)，得出的結(jié)果和文獻(xiàn)的變化趨勢(shì)是一樣的，但是e15的大小卻約等于文獻(xiàn)結(jié)果的一半，我分別用LDA和GGA算過，但結(jié)果差不多，也改變過精度，還是不行，所以我懷疑我的計(jì)算過程有問題，所以請(qǐng)各位幫忙看一下）：
文獻(xiàn)中的結(jié)果：

我計(jì)算的結(jié)果（縱坐標(biāo)的負(fù)值不用管，只看大小，因?yàn)檫@只是與相對(duì)的方向有關(guān)）：

（下面以我計(jì)算靜水壓力P=8GPa時(shí)的步驟為例子說明）
一、優(yōu)化原子坐標(biāo)及晶格常數(shù)
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 3
POTIM = 0.2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.

KPOINTS：(后面采用同一個(gè)KPOINTS）
Kpoints
0
Monkhorst-pack
11 11 11
0 0 0

POSCAR：
PbTiO3_p4mm 3
1.000000000000000
   3.9020000000000000 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.9020000000000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 4.1550000000000000
3 1 1
Direct
  0.5000000000000000  0.5000000000000000  0.1131297125826391
  0.5000000000000000  0.0000000000000000  0.6196488104654303
  0.0000000000000000  0.5000000000000000  0.6196488104654303
  0.5000000000000000  0.5000000000000000  0.5397451459032098
  0.0000000000000000  0.0000000000000000 -0.0047183237261610

POTCAR：LDA （O Ti Pb）（后面采用同一個(gè)POTCAR）

結(jié)果為：
CONTCAR：
PbTiO3_p4mm 3
1.000000000000000
   3.8646149601818020 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.8646149601818020 0.0000000000000000
   0.0000000000000000 0.0000000000000000 4.0560945114966480
O Ti Pb
3 1 1
Direct
  0.5000000000000000  0.5000000000000000  0.1035077853086349
  0.5000000000000000  0.0000000000000000  0.6148831221823548
  0.0000000000000000  0.5000000000000000  0.6148831221823548
  0.5000000000000000  0.5000000000000000  0.5451943700111713
  0.0000000000000000  0.0000000000000000  0.0089857560060325

二、加靜水壓并優(yōu)化原子位置和晶格常數(shù)
cp CONTCAR POSCAR
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 3
PSTRESS = 80
POTIM = 0.2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.

KPOINTS,POTCAR 同上
結(jié)果：
CONTCAR：
PbTiO3_p4mm
1.000000000000000
   3.8347769794747630 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.8347769794747630 0.0000000000000000
   0.0000000000000000 0.0000000000000000 3.8969610358476200
O Ti Pb
3 1 1
Direct
  0.5000000000000000  0.5000000000000000  0.0882145213088690
  0.5000000000000000  0.0000000000000000  0.6050703395668161
  0.0000000000000000  0.5000000000000000  0.6050703395668161
  0.5000000000000000  0.5000000000000000  0.5550015891094247
  0.0000000000000000  0.0000000000000000  0.0340973661386298

三、計(jì)算壓電常數(shù)
先做一次ISIF=2的原子坐標(biāo)馳豫
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.2
NSW = 60; IBRION = 2
EDIFF = 1E-5
EDIFFG = -1E-2
ISIF = 2
PSTRESS = 80
POTIM = 0.2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.
再計(jì)算壓電常數(shù)
CP CONTCAR POSCAR
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
LEPSILON=.TRUE.
IBRION= 8
PSTRESS= 80
ISTART = 0
ICHARG = 2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.

KPOINTS，POTCAR同上
結(jié)果：
e33 = 7.4546  e15 = 9.88125
文獻(xiàn)結(jié)果：e33約等于8，  e15約等于50
也就是說e33的值和文獻(xiàn)差不多，但e15數(shù)值差很大

（另外，我去掉上一個(gè)INCAR中的PSTRESS= 80，計(jì)算結(jié)果一樣）
（去掉PREC = High，用ENCUT =1200計(jì)算，結(jié)果也差不多）
（我改用GGA計(jì)算過，結(jié)果也沒變多少）

各位高手，您看上面的步驟中我有沒有用錯(cuò)參數(shù)的地方？

[ Last edited by donggaomin on 2011-5-1 at 10:52 ]

回復(fù)此樓

好好學(xué)習(xí)

1樓 2011-05-01 10:49:16

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

回帖支持 ( 顯示支持度最高的前 50 名 )

youqin5912

金蟲 (正式寫手)

應(yīng)助: 13 (小學(xué)生)
金幣: 1183.4
散金: 785
紅花: 15
帖子: 397
在線: 594.4小時(shí)
蟲號(hào): 1141030
注冊(cè): 2010-11-07
性別: GG
專業(yè): 同步輻射技術(shù)及其應(yīng)用

樓主你好，請(qǐng)問你添加LEPSILON=.TRUE.后，OUTCAR中的 PIEZOELECTRIC TENSOR for field in x, y, z (C/m^2)和PIEZOELECTRIC TENSOR IONIC CONTR for field in x, y, z (C/m^2)分別代表什么？是e33么？非常感謝

贊一下(1人)

回復(fù)此樓

8樓2014-08-21 15:27:40

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

普通回帖

jghe

金蟲 (著名寫手)

應(yīng)助: 10 (幼兒園)
貴賓: 0.225
金幣: 112.2
散金: 5
紅花: 3
帖子: 1481
在線: 123小時(shí)
蟲號(hào): 322338
注冊(cè): 2007-03-11
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 量子化學(xué)

★ ★
youzhizhe(金幣+2): 謝謝交流。 2011-05-02 09:15:02

我覺得你應(yīng)該用腳小體積的方法加靜水壓，用PSTRESS給出的結(jié)果不是很可靠。用E—V 結(jié)果再根據(jù)狀態(tài)方程給出的壓力應(yīng)該比較可靠。彈性常數(shù)沒算過，應(yīng)該跟體積關(guān)系很密切
如果結(jié)果已經(jīng)收斂了，你用再大的ENCUT沒有任何意義

贊一下(1人)

回復(fù)此樓

2樓2011-05-02 02:31:14

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

youzhizhe(金幣+1, 1ST強(qiáng)帖+1): 謝謝交流。 2011-05-02 19:57:20
youzhizhe(金幣-1, 1ST強(qiáng)帖-1): 操作錯(cuò)誤，抱歉。 2011-05-02 19:58:00

引用回帖:

Originally posted by jghe at 2011-05-02 02:31:14:
我覺得你應(yīng)該用腳小體積的方法加靜水壓，用PSTRESS給出的結(jié)果不是很可靠。用E—V 結(jié)果再根據(jù)狀態(tài)方程給出的壓力應(yīng)該比較可靠。彈性常數(shù)沒算過，應(yīng)該跟體積關(guān)系很密切
如果結(jié)果已經(jīng)收斂了，你用再大的ENCUT沒有任 ...

謝謝！
有兩個(gè)小疑問，還請(qǐng)您解答一下：

疑問一、對(duì)每一個(gè)相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化(ISIF=4)，將優(yōu)化得到再進(jìn)行一次靜態(tài)的總能計(jì)算，得到這10來個(gè)體積點(diǎn)下的總能。

贊一下(2人)

回復(fù)此樓

好好學(xué)習(xí)

3樓2011-05-02 15:10:41

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

youzhizhe(金幣+1): 謝謝交流。以后回帖時(shí)稍微注意就是了，呵呵。 2011-05-02 19:57:39
youzhizhe(金幣-1): 操作錯(cuò)誤，抱歉。以后回帖時(shí)稍微注意就是了。 2011-05-02 19:58:39

引用回帖:

Originally posted by donggaomin at 2011-05-02 15:10:41:

謝謝！
有兩個(gè)小疑問，還請(qǐng)您解答一下：

疑問一、對(duì)每一個(gè)相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化(ISIF=4)，將優(yōu)化得到再進(jìn)行一次靜態(tài)的總能計(jì)算，得 ...

不好意思！還沒有寫完就錯(cuò)按了回車鍵。

  疑問一、以下我對(duì)減小體積加壓的理解正確嗎（四方相）？
            對(duì)每一個(gè)相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化(ISIF=4)，將優(yōu)化得到再進(jìn)行一次靜態(tài)的總能計(jì)算，得到這10來個(gè)體積點(diǎn)下的總能。
         用EOS擬合出P-V后，對(duì)于每一個(gè)P，有一個(gè)V。這個(gè)P就是所加壓強(qiáng)。要計(jì)算其它性質(zhì)時(shí)，就直接用這個(gè)V根據(jù)a/c比例不變計(jì)算出a和c，然后做固定體積的優(yōu)化(ISIF=4)得到CONTCAR，用這個(gè)CONTCAR 直接變成POSCAR，計(jì)算其它性質(zhì)。

疑問二，有些物質(zhì)的相，是在加壓下才出現(xiàn)，比如PbTiO3在零壓下為四方相，加壓到9GPa為單斜相，上面的加壓法能壓出單斜相嗎？

謝謝！

贊一下(2人)

回復(fù)此樓

好好學(xué)習(xí)

4樓2011-05-02 15:28:53

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

xx2008

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 6.3
散金: 65
帖子: 228
在線: 95.2小時(shí)
蟲號(hào): 646702
注冊(cè): 2008-11-05

【答案】應(yīng)助回帖

★
donggaomin(金幣+2): 咱們可以多交流，計(jì)算過程詳見上面說明 2011-05-06 22:45:10
youzhizhe(金幣+1): 謝謝交流。 2011-05-12 19:51:43

引用回帖:

Originally posted by donggaomin at 2011-05-02 15:28:53:
不好意思！還沒有寫完就錯(cuò)按了回車鍵。

疑問一、以下我對(duì)減小體積加壓的理解正確嗎（四方相）？
對(duì)每一個(gè)相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化( ...

我目前也正在計(jì)算壓力下的性質(zhì)，我覺得您的理解很正確，我也正在試試計(jì)算。另外，您能說說您的壓電常數(shù)怎么計(jì)算的嗎？謝謝

贊一下(1人)

回復(fù)此樓

5樓2011-05-04 09:17:29

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

donggaomin

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by xx2008 at 2011-05-04 09:17:29:
我目前也正在計(jì)算壓力下的性質(zhì)，我覺得您的理解很正確，我也正在試試計(jì)算。另外，您能說說您的壓電常數(shù)怎么計(jì)算的嗎？謝謝

再計(jì)算壓電常數(shù)
CP CONTCAR POSCAR
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
LEPSILON=.TRUE.
IBRION= 8
PSTRESS= 80
ISTART = 0
ICHARG = 2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.

我現(xiàn)在是這樣算的，但是在加壓下，的壓電系數(shù)不準(zhǔn)，所以我才出了這一個(gè)貼子，如果你算加壓下的壓電系數(shù)算對(duì)了，請(qǐng)貼出來讓我也學(xué)習(xí)一下，共同進(jìn)步，交流

贊一下

回復(fù)此樓

好好學(xué)習(xí)

6樓2011-05-04 11:35:03

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

future_wl

木蟲 (著名寫手)

應(yīng)助: 87 (初中生)
金幣: 3295.4
散金: 1881
紅花: 5
帖子: 1404
在線: 676.6小時(shí)
蟲號(hào): 950891
注冊(cè): 2010-01-30
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

507204樓: Originally posted by donggaomin at 2011-05-04 11:35:03
再計(jì)算壓電常數(shù)
CP CONTCAR POSCAR
INCAR：
SYSTEM = PbTiO3_p4mm
ENCUT = 800
LEPSILON=.TRUE.
IBRION= 8
PSTRESS= 80
ISTART = 0
ICHARG = 2
PREC = High
LWAVE = .FALSE.
LCHARG =.FALSE.

我現(xiàn)在 ...

你好，壓電系數(shù)是否分電子貢獻(xiàn)和離子貢獻(xiàn)兩部分，你的這個(gè)參數(shù)能將兩個(gè)部分都算出來嗎？
謝謝！

贊一下

回復(fù)此樓

未來就是現(xiàn)在

7樓2012-09-14 16:03:54

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

vasp001

金蟲 (正式寫手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線: 326小時(shí)
蟲號(hào): 727265
注冊(cè): 2009-03-20
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

內(nèi)容已刪除

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

9樓2014-11-27 17:13:18

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

lusi_1993

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5.5
帖子: 12
在線: 14.9小時(shí)
蟲號(hào): 4119842
注冊(cè): 2015-10-05
專業(yè): 理論和計(jì)算化學(xué)

樓主，我最近也在計(jì)算壓電，和你遇到的情況類似，請(qǐng)問樓主解決了嗎？

贊一下

回復(fù)此樓

10樓2016-03-24 09:51:50

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

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

返回列表

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

[考研] 281求調(diào)劑（0805） +3	煙汐憶海 2026-03-16	8/400	2026-03-18 00:47 by 煙汐憶海
[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4	share_joy 2026-03-16	4/200	2026-03-17 21:32 by hmn_wj
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky

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

[求助] 請(qǐng)計(jì)算過加靜水壓下體系的性質(zhì)的高手幫忙看看-----加壓下的壓電系數(shù)

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

» 猜你喜歡

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

【答案】應(yīng)助回帖

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