 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5517 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
donggaomin金蟲 (正式寫手)
|
[求助]
請計算過加靜水壓下體系的性質(zhì)的高手幫忙看看-----加壓下的壓電系數(shù)
|
|||
|
祝各位五一節(jié)快樂！ 下面是我計算壓電系數(shù)的步驟，請各位大俠看一下（用下面的步驟重復(fù)吳志剛的那篇PRL時，得出的結(jié)果和文獻(xiàn)的變化趨勢是一樣的，但是e15的大小卻約等于文獻(xiàn)結(jié)果的一半，我分別用LDA和GGA算過，但結(jié)果差不多，也改變過精度，還是不行，所以我懷疑我的計算過程有問題，所以請各位幫忙看一下）： 文獻(xiàn)中的結(jié)果： 我計算的結(jié)果（縱坐標(biāo)的負(fù)值不用管，只看大小，因?yàn)檫@只是與相對的方向有關(guān)）： （下面以我計算靜水壓力P=8GPa時的步驟為例子說明） 一、優(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：(后面采用同一個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）（后面采用同一個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 三、計算壓電常數(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. 再計算壓電常數(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ù)值差很大 （另外，我去掉上一個INCAR中的PSTRESS= 80，計算結(jié)果一樣） （去掉PREC = High，用ENCUT =1200計算，結(jié)果也差不多） （我改用GGA計算過，結(jié)果也沒變多少） 各位高手，您看上面的步驟中我有沒有用錯參數(shù)的地方？ [ Last edited by donggaomin on 2011-5-1 at 10:52 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 凝聚態(tài)物理 | soft study |
» 猜你喜歡
347求調(diào)劑
已經(jīng)有9人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有8人回復(fù)
-
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有8人回復(fù)
-
考研調(diào)劑求助
已經(jīng)有3人回復(fù)
-
362求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有11人回復(fù)
-
290分調(diào)劑求助
已經(jīng)有5人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有8人回復(fù)
-
085602 307分 求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 急！！求助，vasp計算中 spin up/spin down 電荷密度圖怎么畫
已經(jīng)有11人回復(fù)
- 卡方分布、t分布及F分布的性質(zhì)及應(yīng)用
已經(jīng)有71人回復(fù)
- CASTEP計算出現(xiàn)你們經(jīng)常出現(xiàn)的問題了，過來瞧瞧。
已經(jīng)有11人回復(fù)
- 求過柱子的高手幫忙
已經(jīng)有37人回復(fù)
- 請大家?guī)兔纯矗鈱?dǎo)這樣的回復(fù)是婉拒嗎？
已經(jīng)有7人回復(fù)
- 金屬配合物穩(wěn)定常數(shù)分子模擬計算
已經(jīng)有7人回復(fù)
- 幫忙給看看選取電子衍射圖，給指點(diǎn)指點(diǎn)，是多晶還是單晶，怎么跟文獻(xiàn)不太一樣
已經(jīng)有12人回復(fù)
- 德國老板回信了，請各位幫忙看看，這種有沒有戲？
已經(jīng)有7人回復(fù)
- 加溶劑后計算，不收斂，請高手指點(diǎn)
已經(jīng)有5人回復(fù)
- 做加壓試驗(yàn)時，反應(yīng)時間起始點(diǎn)應(yīng)該怎么算？
已經(jīng)有4人回復(fù)
- 冷等靜壓成型時，升壓階段有必要分段加壓嗎？
已經(jīng)有11人回復(fù)
- 看看別人僅用VASP做DOS計算發(fā)好文章（兩篇PRL）
已經(jīng)有30人回復(fù)
- 從cif文件直接導(dǎo)入的晶體可以用于直接計算么
已經(jīng)有4人回復(fù)
- 請幫忙看看這個服務(wù)器配置單，主要用于Materials-studio計算
已經(jīng)有10人回復(fù)
5樓2011-05-04 09:17:29
2樓2011-05-02 02:31:14
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
謝謝！3樓2011-05-02 15:10:41
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
youzhizhe(金幣+1): 謝謝交流。以后回帖時稍微注意就是了，呵呵。 2011-05-02 19:57:39
youzhizhe(金幣-1): 操作錯誤，抱歉。以后回帖時稍微注意就是了。 2011-05-02 19:58:39
youzhizhe(金幣-1): 操作錯誤，抱歉。以后回帖時稍微注意就是了。 2011-05-02 19:58:39
|
不好意思！還沒有寫完就錯按了回車鍵。 疑問一、以下我對減小體積加壓的理解正確嗎（四方相）？ 對每一個相，優(yōu)化得V0后，改變V，手算出a,c（要保持a/c的值不變）；然后做固定體積的優(yōu)化(ISIF=4)，將優(yōu)化得到再進(jìn)行一次靜態(tài)的總能計算，得到這10來個體積點(diǎn)下的總能。 用EOS擬合出P-V后，對于每一個P，有一個V。這個P就是所加壓強(qiáng)。要計算其它性質(zhì)時，就直接用這個V根據(jù)a/c比例不變計算出a和c，然后做固定體積的優(yōu)化(ISIF=4)得到CONTCAR，用這個CONTCAR 直接變成POSCAR，計算其它性質(zhì)。 疑問二，有些物質(zhì)的相，是在加壓下才出現(xiàn)，比如PbTiO3在零壓下為四方相，加壓到9GPa為單斜相，上面的加壓法能壓出單斜相嗎？ 謝謝！ |
4樓2011-05-02 15:28:53
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 土木304求調(diào)劑 +4 | 兔突突突， 2026-03-31 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
|
[考研] 求調(diào)劑 +4 | 圖鑒212 2026-03-30 | 4/200 |
|
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +4 | zzll03 2026-03-31 | 4/200 |
|
|
|
[考研] 276求調(diào)劑 +3 | 趙久華 2026-03-29 | 3/150 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|