24小时热门版块排行榜    

查看: 904  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见