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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

donggaomin

金蟲(chóng) (正式寫(xiě)手)

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

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

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

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

（下面以我計(jì)算靜水壓力P=8GPa時(shí)的步驟為例子說(shuō)明）
一、優(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
也就是說(shuō)e33的值和文獻(xiàn)差不多，但e15數(shù)值差很大

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

各位高手，您看上面的步驟中我有沒(méi)有用錯(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的回帖

future_wl

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 87 (初中生)
金幣: 3295.4
散金: 1881
紅花: 5
帖子: 1404
在線(xiàn): 676.6小時(shí)
蟲(chóng)號(hào): 950891
注冊(cè): 2010-01-30
性別: GG
專(zhuān)業(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è)部分都算出來(lái)嗎？
謝謝！

贊一下

回復(fù)此樓

未來(lái)就是現(xiàn)在

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

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

查看全部 11 個(gè)回答

jghe

金蟲(chóng) (著名寫(xiě)手)

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 820
紅花: 9
帖子: 565
在線(xiàn): 244.3小時(shí)
蟲(chóng)號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專(zhuān)業(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:
我覺(jué)得你應(yīng)該用腳小體積的方法加靜水壓，用PSTRESS給出的結(jié)果不是很可靠。用E—V 結(jié)果再根據(jù)狀態(tài)方程給出的壓力應(yīng)該比較可靠。彈性常數(shù)沒(méi)算過(guò)，應(yīng)該跟體積關(guān)系很密切
如果結(jié)果已經(jīng)收斂了，你用再大的ENCUT沒(méi)有任 ...

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

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

贊一下(2人)

回復(fù)此樓

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

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

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

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 820
紅花: 9
帖子: 565
在線(xiàn): 244.3小時(shí)
蟲(chóng)號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專(zhuān)業(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:

不好意思！還沒(méi)有寫(xiě)完就錯(cuò)按了回車(chē)鍵。

  疑問(wèn)一、以下我對(duì)減小體積加壓的理解正確嗎（四方相）？
            對(duì)每一個(gè)相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化(ISIF=4)，將優(yōu)化得到再進(jìn)行一次靜態(tài)的總能計(jì)算，得到這10來(lái)個(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ì)。

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

謝謝！

贊一下(2人)

回復(fù)此樓

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

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

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

查看全部 11 個(gè)回答

普通表情龍兔虎貓

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

[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 0856求調(diào)劑 +7	楒桉 2026-03-28	7/350	2026-03-29 08:28 by fmesaito
[考研] 329求調(diào)劑 +7	星野? 2026-03-26	7/350	2026-03-29 06:43 by 544594351
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚(yú) 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來(lái)啊
[考研] 311（085601）求調(diào)劑 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專(zhuān)碩-298 +4	手機(jī)用戶(hù) 2026-03-26	4/200	2026-03-28 07:17 by mmm just
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +5	dongfang59 2026-03-27	5/250	2026-03-28 04:53 by wxiongid
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚(yú)的貓
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣(mài)魚(yú) 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4	我愛(ài)學(xué)電池 2026-03-23	4/200	2026-03-25 00:59 by 1027_324
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月