24小时热门版块排行榜    

查看: 892  |  回复: 0
【悬赏金币】回答本帖问题,作者小茂茂.将赠送您 20 个金币

小茂茂.

新虫 (初入文坛)

[求助] BaTiO3扩胞以后掺杂怎么计算压电性

下面是我的计算步骤:
(1)优化
INCAR:
SYSTEM = BTZO
ISTART =  1            
ISPIN  =  1            
LREAL  = .FALSE.      
ENCUT  =  500        
PREC   =  Normal      
LWAVE  = .FALSE.        
LCHARG = .FALSE.        
ADDGRID= .TRUE.      
ISMEAR =  0            
SIGMA  =  0.05         
EDIFF  =  1E-05        
NSW    =  100         
IBRION =  2            
ISIF   =  2            
EDIFFG = -2E-02     
(部分)POSCAR(3x3x2):
BaTiZrO3
1.0
       12.1066999435         0.0000000000         0.0000000000
        0.0000000000        12.1066999435         0.0000000000
        0.0000000000         0.0000000000         8.0712003708
   Ba   Ti    O   Zr
   18   17   54    1
Direct
     0.000000000         0.000000000         0.000000000
     0.333330005         0.000000000         0.000000000
     0.666670024         0.000000000         0.000000000
KPOINTS:
K-Spacing Value to Generate K-Mesh: 0.020
0
Monkhorst-Pack
   4   4   6
0.0  0.0  0.0
POTCAR(略)
(2)自洽(略)
(3)压电性(0压力0电场)
INCAR:
SYSTEM = BTZO
ISTART =  1            
PREC   =  Normal           
ENCUT  =  500        
PREC   =  Normal      
LWAVE  = .FALSE.        
LCHARG = .FALSE.        
ADDGRID= .TRUE.        
ISMEAR =  0            
SIGMA  =  0.05        
LORBIT =  10
NELM   =  60         
EDIFF  =  1E-05        
LEPSILON=.TRUE.
IBRION= 8
最后结果如下:
ELASTIC MODULI IONIC CONTR (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX          -8.4988      0.1400      0.1400      0.0000     -0.0000     -0.0000
YY           0.1400     -8.4988      0.1400      0.0000      0.0000     -0.0000
ZZ           0.1400      0.1400     -8.4988     -0.0000      0.0000      0.0000
XY           0.0000      0.0000     -0.0000     -4.7557      0.0000      0.0000
YZ          -0.0000      0.0000      0.0000      0.0000     -4.7557      0.0000
ZX          -0.0000     -0.0000      0.0000      0.0000      0.0000     -4.7557
--------------------------------------------------------------------------------


PIEZOELECTRIC TENSOR IONIC CONTR  for field in x, y, z        (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.00000    -0.00000    -0.00000    -0.00000    -0.00000    -0.00000
  y     0.00000    -0.00000     0.00000     0.00000     0.00000     0.00000
  z     0.00000     0.00000    -0.00000     0.00000     0.00000    -0.00000
问题1:为什么要扩胞为3x3x2,因为参考了李强博士的毕业论文(无铅压电材料压电特性的第一性原理研究)中说到在掺杂Zr在10%以下才有压电性,那么为什么我的压电张量全部为0?还有为什么我计算的没有离子部和电子部?
问题2:压电性d33、d31、d15是怎么计算出来的?文献中单位是C/N。怎么转换的?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小茂茂. 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:44 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:43 by 3VGOm6paSgqv
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:32 by 3VGOm6paSgqv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 20:32 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 20:24 by 3VGOm6paSgqv
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:21 by 3VGOm6paSgqv
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 12/600 2026-08-02 20:21 by drbart
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:11 by 3VGOm6paSgqv
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 cBGhfJFuNMca 2026-08-02 17/850 2026-08-02 20:03 by 3VGOm6paSgqv
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cBGhfJFuNMca 2026-08-02 17/850 2026-08-02 19:54 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 18:34 by ca3DVAx4OB19
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见