 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (5018)
- >導(dǎo)師招生 (1089)
- >考研 (1087)
- >文獻(xiàn)求助 (270)
- >考博 (112)
- >蟲(chóng)友互識(shí) (112)
- >休閑灌水 (107)
- >碩博家園 (87)
- >招聘信息布告欄 (77)
- >論文道賀祈福 (57)
- >基金申請(qǐng) (43)
- >博后之家 (36)
- >SciFinder/Reaxys (32)
- >論文投稿 (27)
- >綠色求助(高懸賞) (21)
- >有機(jī)資源 (18)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5502 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
donggaomin金蟲(chóng) (正式寫(xiě)手)
|
[求助]
請(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ù)此樓» 收錄本帖的淘帖專輯推薦
 凝聚態(tài)物理 | soft study |
» 猜你喜歡
274求調(diào)劑
已經(jīng)有7人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有5人回復(fù)
-
網(wǎng)絡(luò)空間安全0839招調(diào)劑
已經(jīng)有4人回復(fù)
-
329求調(diào)劑
已經(jīng)有3人回復(fù)
-
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有8人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 急��！求助，vasp計(jì)算中 spin up/spin down 電荷密度圖怎么畫(huà)
已經(jīng)有11人回復(fù)
- 卡方分布、t分布及F分布的性質(zhì)及應(yīng)用
已經(jīng)有71人回復(fù)
- CASTEP計(jì)算出現(xiàn)你們經(jīng)常出現(xiàn)的問(wèn)題了，過(guò)來(lái)瞧瞧。
已經(jīng)有11人回復(fù)
- 求過(guò)柱子的高手幫忙
已經(jīng)有37人回復(fù)
- 請(qǐng)大家?guī)兔纯矗鈱?dǎo)這樣的回復(fù)是婉拒嗎？
已經(jīng)有7人回復(fù)
- 金屬配合物穩(wěn)定常數(shù)分子模擬計(jì)算
已經(jīng)有7人回復(fù)
- 幫忙給看看選取電子衍射圖，給指點(diǎn)指點(diǎn)，是多晶還是單晶，怎么跟文獻(xiàn)不太一樣
已經(jīng)有12人回復(fù)
- 德國(guó)老板回信了，請(qǐng)各位幫忙看看，這種有沒(méi)有戲？
已經(jīng)有7人回復(fù)
- 加溶劑后計(jì)算，不收斂，請(qǐng)高手指點(diǎn)
已經(jīng)有5人回復(fù)
- 做加壓試驗(yàn)時(shí)，反應(yīng)時(shí)間起始點(diǎn)應(yīng)該怎么算？
已經(jīng)有4人回復(fù)
- 冷等靜壓成型時(shí)，升壓階段有必要分段加壓?jiǎn)幔?/a>
已經(jīng)有11人回復(fù)
- 看看別人僅用VASP做DOS計(jì)算發(fā)好文章（兩篇PRL）
已經(jīng)有30人回復(fù)
- 從cif文件直接導(dǎo)入的晶體可以用于直接計(jì)算么
已經(jīng)有4人回復(fù)
- 請(qǐng)幫忙看看這個(gè)服務(wù)器配置單，主要用于Materials-studio計(jì)算
已經(jīng)有10人回復(fù)
donggaomin
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時(shí)
- 蟲(chóng)號(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
youzhizhe(金幣-1): 操作錯(cuò)誤，抱歉。以后回帖時(shí)稍微注意就是了。 2011-05-02 19:58:39
|
不好意思！還沒(méi)有寫(xiě)完就錯(cuò)按了回車鍵。 疑問(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為單斜相，上面的加壓法能壓出單斜相嗎？ 謝謝！ |
4樓2011-05-02 15:28:53
2樓2011-05-02 02:31:14
donggaomin
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時(shí)
- 蟲(chóng)號(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
youzhizhe(金幣-1, 1ST強(qiáng)帖-1): 操作錯(cuò)誤，抱歉。 2011-05-02 19:58:00
 謝謝！有兩個(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)下的總能。 |
3樓2011-05-02 15:10:41
xx2008
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6.3
- 散金: 65
- 帖子: 228
- 在線: 95.2小時(shí)
- 蟲(chóng)號(hào): 646702
- 注冊(cè): 2008-11-05
5樓2011-05-04 09:17:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 312求調(diào)劑 +3 | 上岸吧ZJY 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑 +3 | 鈕恩雪 2026-03-25 | 3/150 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +10 | 碩星赴 2026-03-23 | 10/500 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +4 | 呼呼？~+123456 2026-03-24 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +6 | 匹克i 2026-03-23 | 6/300 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +5 | 谷雨上岸 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？ +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開(kāi)頭均可，過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 315分，誠(chéng)求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|