| 查看: 5532 | 回复: 10 | ||||
donggaomin金虫 (正式写手)
|
[求助]
请计算过加静水压下体系的性质的高手帮忙看看-----加压下的压电系数
|
|||
|
祝各位五一节快乐！ 下面是我计算压电系数的步骤，请各位大侠看一下（用下面的步骤重复吴志刚的那篇PRL时，得出的结果和文献的变化趋势是一样的，但是e15的大小却约等于文献结果的一半，我分别用LDA和GGA算过，但结果差不多，也改变过精度，还是不行，所以我怀疑我的计算过程有问题，所以请各位帮忙看一下）： 文献中的结果：  我计算的结果（纵坐标的负值不用管，只看大小，因为这只是与相对的方向有关）：  （下面以我计算静水压力P=8GPa时的步骤为例子说明） 一、优化原子坐标及晶格常数 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） 结果为： 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 二、加静水压并优化原子位置和晶格常数 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 同上 结果： 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 三、计算压电常数 先做一次ISIF=2的原子坐标驰豫 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. 再计算压电常数 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同上 结果： e33 = 7.4546 e15 = 9.88125 文献结果：e33约等于8， e15约等于50 也就是说e33的值和文献差不多，但e15数值差很大 （另外，我去掉上一个INCAR中的PSTRESS= 80，计算结果一样） （去掉PREC = High，用ENCUT =1200计算，结果也差不多） （我改用GGA计算过，结果也没变多少） 各位高手，您看上面的步骤中我有没有用错参数的地方？ [ Last edited by donggaomin on 2011-5-1 at 10:52 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 凝聚态物理 | soft study |
» 猜你喜欢
080500求调剂
已经有5人回复
-
269电子信息求调剂，可转专业
已经有5人回复
-
285求调剂
已经有15人回复
-
265求调剂
已经有8人回复
-
274求调剂求调剂
已经有5人回复
-
284求调剂
已经有4人回复
-
土木水利专硕276分求调剂
已经有4人回复
-
材料考研求调剂总分280
已经有9人回复
-
285求调剂
已经有6人回复
-
材料调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急！！求助，vasp计算中 spin up/spin down 电荷密度图怎么画
已经有11人回复
- 卡方分布、t分布及F分布的性质及应用
已经有71人回复
- CASTEP计算出现你们经常出现的问题了，过来瞧瞧。
已经有11人回复
- 求过柱子的高手帮忙
已经有37人回复
- 请大家帮忙看看，外导这样的回复是婉拒吗？
已经有7人回复
- 金属配合物稳定常数分子模拟计算
已经有7人回复
- 帮忙给看看选取电子衍射图，给指点指点，是多晶还是单晶，怎么跟文献不太一样
已经有12人回复
- 德国老板回信了，请各位帮忙看看，这种有没有戏？
已经有7人回复
- 加溶剂后计算，不收敛，请高手指点
已经有5人回复
- 做加压试验时，反应时间起始点应该怎么算？
已经有4人回复
- 冷等静压成型时，升压阶段有必要分段加压吗？
已经有11人回复
- 看看别人仅用VASP做DOS计算发好文章（两篇PRL）
已经有30人回复
- 从cif文件直接导入的晶体可以用于直接计算么
已经有4人回复
- 请帮忙看看这个服务器配置单，主要用于Materials-studio计算
已经有10人回复
youqin5912
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1183.4
- 散金: 785
- 红花: 15
- 帖子: 397
- 在线: 594.4小时
- 虫号: 1141030
- 注册: 2010-11-07
- 性别: GG
- 专业: 同步辐射技术及其应用
楼主你好，请问你添加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么？非常感谢
8楼2014-08-21 15:27:40
2楼2011-05-02 02:31:14
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 673.8
- 散金: 820
- 红花: 9
- 帖子: 566
- 在线: 244.3小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-05-02 15:10:41
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 673.8
- 散金: 820
- 红花: 9
- 帖子: 566
- 在线: 244.3小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
youzhizhe(金币+1): 谢谢交流。以后回帖时稍微注意就是了，呵呵。 2011-05-02 19:57:39
youzhizhe(金币-1): 操作错误，抱歉。以后回帖时稍微注意就是了。 2011-05-02 19:58:39
youzhizhe(金币-1): 操作错误，抱歉。以后回帖时稍微注意就是了。 2011-05-02 19:58:39
|
不好意思！还没有写完就错按了回车键。 疑问一、以下我对减小体积加压的理解正确吗（四方相）？ 对每一个相，优化得V0后，改变V，手算出a,c（要保持a/c的值不变）；然后做固定体积的优化(ISIF=4)，将优化得到再进行一次静态的总能计算，得到这10来个体积点下的总能。 用EOS拟合出P-V后，对于每一个P，有一个V。这个P就是所加压强。要计算其它性质时，就直接用这个V根据a/c比例不变计算出a和c，然后做固定体积的优化(ISIF=4)得到CONTCAR，用这个CONTCAR 直接变成POSCAR，计算其它性质。 疑问二，有些物质的相，是在加压下才出现，比如PbTiO3在零压下为四方相，加压到9GPa为单斜相，上面的加压法能压出单斜相吗？ 谢谢！ |
4楼2011-05-02 15:28:53
5楼2011-05-04 09:17:29
donggaomin
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 673.8
- 散金: 820
- 红花: 9
- 帖子: 566
- 在线: 244.3小时
- 虫号: 853446
- 注册: 2009-09-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2011-05-04 11:35:03
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2012-09-14 16:03:54
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2014-11-27 17:13:18
10楼2016-03-24 09:51:50
