 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5501 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
donggaomin金蟲 (正式寫手)
|
[求助]
請計(jì)算過加靜水壓下體系的性質(zhì)的高手幫忙看看-----加壓下的壓電系數(shù)
|
|||
|
祝各位五一節(jié)快樂！ 下面是我計(jì)算壓電系數(shù)的步驟，請各位大俠看一下（用下面的步驟重復(fù)吳志剛的那篇PRL時，得出的結(jié)果和文獻(xiàn)的變化趨勢是一樣的，但是e15的大小卻約等于文獻(xiàn)結(jié)果的一半，我分別用LDA和GGA算過，但結(jié)果差不多，也改變過精度，還是不行，所以我懷疑我的計(jì)算過程有問題，所以請各位幫忙看一下）： 文獻(xiàn)中的結(jié)果： 我計(jì)算的結(jié)果（縱坐標(biāo)的負(fù)值不用管，只看大小，因?yàn)檫@只是與相對的方向有關(guān)）： （下面以我計(jì)算靜水壓力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 三、計(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ù)值差很大 （另外，我去掉上一個INCAR中的PSTRESS= 80，計(jì)算結(jié)果一樣） （去掉PREC = High，用ENCUT =1200計(jì)算，結(jié)果也差不多） （我改用GGA計(jì)算過，結(jié)果也沒變多少） 各位高手，您看上面的步驟中我有沒有用錯參數(shù)的地方？ [ Last edited by donggaomin on 2011-5-1 at 10:52 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 凝聚態(tài)物理 | soft study |
» 猜你喜歡
材料與化工考研調(diào)劑
已經(jīng)有5人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料292調(diào)劑
已經(jīng)有5人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有4人回復(fù)
-
收08調(diào)劑生
已經(jīng)有3人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
086003食品工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 急！！求助，vasp計(jì)算中 spin up/spin down 電荷密度圖怎么畫
已經(jīng)有11人回復(fù)
- 卡方分布、t分布及F分布的性質(zhì)及應(yīng)用
已經(jīng)有71人回復(fù)
- CASTEP計(jì)算出現(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ì)算
已經(jīng)有7人回復(fù)
- 幫忙給看看選取電子衍射圖，給指點(diǎn)指點(diǎn)，是多晶還是單晶，怎么跟文獻(xiàn)不太一樣
已經(jīng)有12人回復(fù)
- 德國老板回信了，請各位幫忙看看，這種有沒有戲？
已經(jīng)有7人回復(fù)
- 加溶劑后計(jì)算，不收斂，請高手指點(diǎn)
已經(jīng)有5人回復(fù)
- 做加壓試驗(yàn)時，反應(yīng)時間起始點(diǎn)應(yīng)該怎么算？
已經(jīng)有4人回復(fù)
- 冷等靜壓成型時，升壓階段有必要分段加壓嗎？
已經(jīng)有11人回復(fù)
- 看看別人僅用VASP做DOS計(jì)算發(fā)好文章（兩篇PRL）
已經(jīng)有30人回復(fù)
- 從cif文件直接導(dǎo)入的晶體可以用于直接計(jì)算么
已經(jīng)有4人回復(fù)
- 請幫忙看看這個服務(wù)器配置單，主要用于Materials-studio計(jì)算
已經(jīng)有10人回復(fù)
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
謝謝！
3樓2011-05-02 15:10:41
2樓2011-05-02 02:31:14
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 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)的總能計(jì)算，得到這10來個體積點(diǎn)下的總能。 用EOS擬合出P-V后，對于每一個P，有一個V。這個P就是所加壓強(qiáng)。要計(jì)算其它性質(zhì)時，就直接用這個V根據(jù)a/c比例不變計(jì)算出a和c，然后做固定體積的優(yōu)化(ISIF=4)得到CONTCAR，用這個CONTCAR 直接變成POSCAR，計(jì)算其它性質(zhì)。 疑問二，有些物質(zhì)的相，是在加壓下才出現(xiàn)，比如PbTiO3在零壓下為四方相，加壓到9GPa為單斜相，上面的加壓法能壓出單斜相嗎？ 謝謝！ |
4樓2011-05-02 15:28:53
5樓2011-05-04 09:17:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +9 | 迷人的哈哈 2026-03-23 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 12/600 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +18 | chaojifeixia 2026-03-19 | 20/1000 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +3 | WW.' 2026-03-23 | 5/250 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|