vasp计算压电系数_求助DFTP算出来的压电系数
我用的VASP5.2 下的LEPSILON-tag算的压电系数和BEC等东西,之前发过类似的贴,不过这次想一起贴上来在好好求助一下:
我的INCAR
SYSTEM=PLPK
ENCUT=400
PREC=Normal
ISMEAR=0
SIGMA=0.05
EDIFF=1E-3
ISTART=0
ICHARG=2
IBRION=8
POTIM=0.1
LEPSILON=.TRUE.
LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 -1 3 -1 -1 -1 -1 -1 -1
LDAUU=0 0 6 0 0 0 0 0 0
LDAUJ=0 0 0 0 0 0 0 0 0
LDAUPRINT=2
算的是一个20个原子的四方相结构,含极性不连续的。主要算它的极化强度和压电特性在应变条件下的变化,现在只是OUTCAR中的结果还差那么点儿分析,想多请教下算过相关内容的前辈和同行,希望不吝赐教。
比如0应变的结构,我从OUTCAR中找到的有用内容如下:(目前提取出来了弹性模量,BEC和压电张量矩阵,还有个介电张量)
MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT)
---------------------------------------------
7.562 0.000 0.000
0.000 7.562 0.000
0.000 0.000 6.072
---------------------------------------------
PIEZOELECTRIC TENSOR for field in x, y, z (e Angst)
XX YY ZZ XY YZ ZX
--------------------------------------------------------------------------------
x -0.57655 0.53337 1.41549 0.00000 0.00000 -2.48765
y 0.74595 -0.80028 0.02920 0.00000 -2.29471 0.00000
z -2.13003 -2.13003 -1.71919 0.00000 0.00000 0.00000
PIEZOELECTRIC TENSOR 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.16001
y 0.00000 0.00000 0.00000 0.00000 -0.16005 0.00000
z -0.14256 -0.14256 -0.11506 0.00000 0.00000 0.00000
BORN EFFECTIVE CHARGES (in e, cummulative output)
-------------------------------------------------
ion 1
1 4.16337 0.00000 -0.35555
2 0.00000 4.25218 -0.00273
3 0.00000 0.00000 3.08453
ion 2
1 5.65192 0.00000 -1.45143
2 0.00000 5.62322 0.02908
3 0.00000 0.00000 6.09140
ion 3
1 4.59437 0.00000 -0.77293
2 0.00000 4.57861 0.02654
3 0.00000 0.00000 4.69899
ion 4
1 2.97250 0.00000 -0.46660
2 0.00000 2.96029 0.01771
3 0.00000 0.00000 4.07500
ion 5
1 4.17827 0.00000 -0.05032
2 0.00000 4.15128 0.00691
3 0.00000 0.00000 3.68609
ion 6
1 6.14756 0.00000 0.98037
2 0.00000 6.12975 -0.01380
3 0.00000 0.00000 4.99445
ion 7
1 1.18633 0.00000 0.29818
2 0.00000 1.18496 -0.00823
3 0.00000 0.00000 1.34329
ion 8
1 8.15684 0.00000 0.91563
2 0.00000 7.88065 -0.02060
3 0.00000 0.00000 6.62523
ion 9
1 -2.67478 0.00000 0.42944
2 0.00000 -4.21338 0.03375
3 0.00000 0.00000 -2.01709
ion 10
1 -4.10024 0.00000 0.31148
2 0.00000 -2.70871 -0.04404
3 0.00000 0.00000 -2.01709
ion 11
1 -3.29422 0.00000 0.52821
2 0.00000 -3.19993 -0.00621
3 0.00000 0.00000 -4.46866
ion 12
1 -3.03210 0.00000 0.35551
2 0.00000 -2.55061 0.00564
3 0.00000 0.00000 -3.02006
ion 13
1 -2.56321 0.00000 0.32517
2 0.00000 -2.93166 -0.02865
3 0.00000 0.00000 -3.02006
ion 14
1 -2.59899 0.00000 0.66267
2 0.00000 -2.64079 -0.01770
3 0.00000 0.00000 -3.98642
ion 15
1 -2.45927 0.00000 -0.27148
2 0.00000 -4.36881 0.03943
3 0.00000 0.00000 -1.74718
ion 16
1 -4.39566 0.00000 -0.19074
2 0.00000 -2.44264 -0.03256
3 0.00000 0.00000 -1.74718
ion 17
1 -3.00654 0.00000 -0.27400
2 0.00000 -3.01064 -0.00666
3 0.00000 0.00000 -3.64142
ion 18
1 -2.20096 0.00000 0.01679
2 0.00000 -5.89128 0.02448
3 0.00000 0.00000 -1.75444
ion 19
1 -6.05680 0.00000 -0.01867
2 0.00000 -2.16964 -0.03114
3 0.00000 0.00000 -1.75444
ion 20
1 -1.57589 0.00000 -0.92688
2 0.00000 -1.56198 0.02343
3 0.00000 0.00000 -4.92454
另外还有个离子贡献的介电张量:
MACROSCOPIC STATIC DIELECTRIC TENSOR IONIC CONTRIBUTION
---------------------------------------------
0.155 -0.004 -0.114
-0.004 0.000 0.001
-0.114 0.001 7.979
---------------------------------------------
那么问题来了:
1、首先是两个压电张量,为什么除了单位不一样,它们输出的格式也不一样,后者好多项是0,请问那是eij 的矩阵吗?可不可以换算成dij 的矩阵?
2、离子的介电张量贡献有用吗?或者说这部分贡献分析什么东西有用?
目前就这俩问题,请大神们帮看看吧!!!
vasp计算压电系数_求助DFTP算出来的压电系数相关推荐
- vasp计算压电系数_求助DFTP算出来的压电系数 - 第一原理 - 小木虫 - 学术 科研 互动社区...
我用的VASP5.2 下的LEPSILON-tag算的压电系数和BEC等东西,之前发过类似的贴,不过这次想一起贴上来在好好求助一下: 我的INCAR SYSTEM=PLPK ENCUT=400 PRE ...
- VASP计算笔记_声子谱计算
声子谱计算 主要使用phonopy软件计算声子谱.版本:phonopy-2.1.3 密度泛函微软理论/线性响应方法 (DFPT) 必要的输入文件: INCAR KPOINTS POSCAR-unitc ...
- VASP计算笔记_弹性常数计算
二阶弹性常数计算 主要采用能量–应变方法 AELAS软件计算 必要的文件 elas_strain.sh elas_cal.sh POTCAR INPOS 优化后的POSCAR.可在脚本设置或者单独给出 ...
- java语言编程三角形图形_编程题:编写程序输入三角形的3条边长,计算并输出... 求助一道JAVA编程题:编写一个类似记事本的图形用......
导航:网站首页 > 编程题:编写程序输入三角形的3条边长,计算并输出... 求助一道JAVA编程题:编写一个类似记事本的图形用... 编程题:编写程序输入三角形的3条边长,计算并输出... 求助 ...
- Java计算时间差_传统的SimpleDateFormat类
Java计算时间差_传统的SimpleDateFormat类 SimpleDateFormat simpleFormat = new SimpleDateFormat("yyyy-MM-dd ...
- 匿名飞控TI版_姿态解算
匿名飞控TI版_姿态解算 准备电赛 准备大创 先看看匿名姿态解算的代码 文章目录 匿名飞控TI版_姿态解算 一,姿态解算原理 1.介绍 2.方向余弦矩阵 (1)方向余弦 (2)DCM矩阵 3.欧拉角 ...
- 使用VASP计算材料能带和态密度
前言 一.前期准备 二.计算流程 1.静态自洽(scf) 2.能带计算 3.态密度计算 总结 前言 能带结构是目前采用第一性原理计算所得到的常用信息,可用来结合解释金属.半导体和绝缘体的区别.能带用来 ...
- VASP计算AgGaS2能带及态密度及光学性质
第一步:构型优化 1.准备四个输入文件 POSCAR INCAR POTCAR KPOINTS POSCAR: 从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包.这样, ...
- VASP计算光学性质
PS:请见文末的打赏选项. VASP计算光学性质的INCAR是根据HSE06计算的自洽INCAR来改的 SCF INCAR: SYSTEM = System ISTART = 1 ICHARG = 1 ...
最新文章
- llist对象两个属性相乘在相加_Java8使用stream实现list中对象属性的合并(去重并求和)...
- 进程线程001 进程线程结构体和KPCR
- JVM监控-命令行篇
- 错误记录​:使用CLLocationManager不走代理方法
- QT每日一练day9:设计师界面
- Redis学习总结(23)——Redis集群化方案对比:Codis、Twemproxy、Redis Cluster
- Quartz 触发器、过期触发策略 、排它日历、数据持久化
- 好听的歌曲单片机c语言程序,单片机音乐播放C语言程序
- 动态链接库(.dll) 动态导入库(.lib) 静态链接库(.lib)
- 拿它们练Python爬虫,是在法律边缘试探吗?爬虫圈香饽饽之视频网站的评论区采集
- NFT 推荐|Funguys 王国作品集第一弹
- 假如生产环境出现CPU占用过高,如何排查与定位
- 计算机二级程序设计提交,程序设计方法与风格(计算机二级复习指导)
- 《Android软件安全权威指南》原生程序分析 阅读笔记
- 设计模式之(九)--包装模式(decorator)
- 【研讨会现场】广州天嵌科技受邀参加2014年TI9月深圳研讨会
- hbase命令使用笔记
- 英语单词:through; fall through; write through
- c语言的文件打开模式,c语言文件打开模式 - osc_b2jf5efr的个人空间 - OSCHINA - 中文开源技术交流社区...
- 动态规划之钱币兑换问题
热门文章
- 2k的地址范围 计算机组成原理,计算机组成原理课后习题
- 难道主键除了自增就是GUID?支持k8s等分布式场景下的id生成器了解下
- linux centos安装scp,CentOS使用不了scp命令怎么办?
- CabloyJS微信模块、企业微信模块已出齐
- three.js obj模型的mtl材质贴图不显示
- c++串口配置及DCB结构体
- c++智能指针简单使用
- android ip查看工具,安卓手机查看IP地址的两种方法
- 如何快速梳理领域文献
- 【案例分享】让新时代教育发展与“数”俱进