我用的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算出来的压电系数 - 第一原理 - 小木虫 - 学术 科研 互动社区...相关推荐

  1. vasp服务器中断,求助VASP能带计算的中断原因 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    我也遇到类似情况,发现服务器内存没有用到多少?是不是并行环境有问题?谢谢! DA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and K ...

  2. vasp 模拟退火_【转】vasp的分子动力学模拟 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    vasp做分子动力学的好处,由于vasp是近些年开发的比较成熟的软件,在做电子scf速度方面有较好的优势. 缺点:可选系综太少. 尽管如此,对于大多数有关分子动力学的任务还是可以胜任的. 主要使用的系 ...

  3. linux根目录下并行计算,【求助】MS在Linux下的并行计算问题 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    先在win下生成计算castep所需的文件,param和cell,导入Linux系统. 在计算的时候,如果不指定cpu的个数,直接"./RunCastep.sh name",就能计 ...

  4. 超级计算机 任务提交,vasp在超算中心的任务提交 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    这是我用来测试晶格常数的脚本,最后提取总能E和晶格常数a,并写出来: #!/bin/sh #SBATCH -x xx #SBATCH -x xx .............. for i in  0. ...

  5. linux使用gaussian提交命令,在linux系统提交gaussian任务失败，求助 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    设置用户环境后,提交任务. 输入文件是 %chk=hcccc.ch %mem=500MB %nproc=1 #p b3lyp/sto-3g opt hcccc 0 1 C C              ...

  6. linux安装软件报失败,【求助】linux下安装MS软件不成功 - 第一原理 - 小木虫 - 学术 科研 互动社区...

    以下拷下来是安装过程,感觉就安装了个LicensePack,不知道步骤对不对,问题出在哪里.除了没有HP-MPI,以下安装过程对不对啊?请大家帮忙. ------------------------- ...

  7. openmpi参数_使用gfortan，openmpi以及BLAS, BLACS, Lapack, Scalapack编译SIESTA - 第一原理 - 小木虫 - 学术 科研 互动社区...

    1.编译openmpi (http://www.open-mpi.org/software/ompi/v1.6/) (1) ./configure --prefix=/home/x_guazh/sof ...

  8. java 红外光谱数据库_免费的谱图数据库20个 - 晶体 - 小木虫 - 学术 科研 互动社区...

    [名称]ChemExper化学品目录CDD (包括MSDS.5000张红外谱图) [资源简介] The database contains currently more than 70 000 che ...

  9. arcgis 同名图层合并_【求助】ArcGIS中怎样合并图层 - 地学 - 小木虫 - 学术 科研 互动社区...

    这个问题 最简单的方式是用editor中的切割多边形功能  可以把一个矢量文件按照你定义的线分割  可直接在属性表中算出面积  要是想直接用一个线文件去切割面文件的话 操作有点复杂  下面是用edit ...

最新文章

  1. 中文版-He Knows My Name(他知道我的名字)-祂认得我-陈熙（音乐河4）
  2. HVM guest配置文件
  3. linux哪些分区用固态硬盘,关于linux:简单聊聊固态硬盘的物理结构
  4. struts2教程（7）--拦截器
  5. easyui from 缓存问题处理
  6. 北斗导航 | 卫星导航基础知识（坐标系）
  7. Flex 幻灯片播放
  8. 网络基础：网络IP 、子网掩码 、路由器 、DNS知识分享
  9. php tp3 操作绑定到类,操作绑定到类 · ThinkPHP3.2.3完全开发手册 · 看云
  10. 人工智能平台 H2O.ai
  11. event mpm php,Ubuntu Apache 切换到php-fpm+mpm_event模式
  12. 基于微信朋友圈实现时间显示样式
  13. 使用fiddler代理，手机无法上网
  14. 软件设计第一步——分离关注点和单一职责原则
  15. html在搜索按钮中加放大镜,点击放大镜搜索图标
  16. 计算机屏幕的显示分辨率与什么有关,计算机屏幕分辨率高低主要跟什么有关？...
  17. MySQL　一起重新认识下count(*) ,count(1),count(id)
  18. 现代经济感悟(4)-为什么一分钱一分货，两分货却要一块钱
  19. 旺旺号userid转换店铺shopid和评分性别
  20. 持久记录精彩时刻，汉印CP4000体验

热门文章

  1. python获取摄像头型号_python opencv设置摄像头分辨率以及各个参数的方法_python
  2. matlab绘图崩溃，重新采用硬件加速绘图
  3. 如何启动Intel VT-x
  4. matplotlib官方中文手册pdf下载
  5. win10黑屏Duilib资源加载失败
  6. python官方文档学习_Python3.5.2官方文档学习备忘录
  7. Maya角色UV展平技巧笔记
  8. 【加拿大签证】加拿大政府指定的签证办理中国体检医院一览【2019官方最新版，加拿大签证体检必看】
  9. xmind 8系列便携版：关联.xmind文件的打开方式后，打开文件会在当前文件夹下产生configuration等子目录的问题解决办法
  10. httprunner（9）运行测试用例的方式总结