| 查看: 2994 | 回复: 5 | |||||||||||
Gina88版主
|
[求助]
压电计算-收集计算方法！
|
|
（下面以我计算静水压力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计算过，结果也没变多少） 各位高手，您看上面的步骤中我有没有用错参数的地方？ 回帖： 我觉得你应该用脚小体积的方法加静水压，用PSTRESS给出的结果不是很可靠。用E—V 结果再根据状态方程给出的压力应该比较可靠。弹性常数没算过，应该跟体积关系很密切 回帖： 疑问一、以下我对减小体积加压的理解正确吗（四方相）？ 对每一个相，优化得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为单斜相，上面的加压法能压出单斜相吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 仿真建模与计算 | 第一性原理相关文档 | 第一性原理计算-tg | 第一性原理计算经验 |
| 有效质量 | 计算—磁性 | 计算化学 | soft study |
| 第一性原理计算 |
» 猜你喜欢
316求调剂
已经有3人回复
-
342求调剂
已经有11人回复
-
319求调剂
已经有4人回复
-
【求调剂】085601材料工程专硕 | 总分272 |
已经有10人回复
-
26考研调剂
已经有3人回复
-
294分080500材料科学与工程求调剂
已经有16人回复
-
调剂
已经有3人回复
-
一志愿郑大材料工程290求调剂
已经有16人回复
-
273求调剂
已经有16人回复
-
材料调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 电子轨道能量用什么方法计算比较好？
已经有6人回复
- 分数占据的Cu2Se晶体的热电性质模拟（Nat. Mat+JACS+CM）
已经有10人回复
- 请教压电原理的计算（有点类似电容）。
已经有7人回复
- 压电PZT谐振频率
已经有4人回复
- [原创]一种基于能带计算结果的载流子有效质量计算方法
已经有174人回复
- 压电陶瓷同一面的电容计算
已经有7人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- 第一性原理计算硬度时Eg的计算方法
已经有7人回复
- 第一原理软件的两种磁性计算方法
已经有14人回复
- 请计算过加静水压下体系的性质的高手帮忙看看-----加压下的压电系数
已经有10人回复
- 电容器电容计算方法讨论
已经有56人回复
- 【求助】跪求有机物的荧光寿命测量及计算方法！！！
已经有5人回复
- 【求助】求助关于χ2和P值的计算方法
已经有12人回复
- 【求助】关于检测限的计算方法！
已经有15人回复
- 【求助】请教有效质量和载流子浓度计算方法问题
已经有7人回复
- 【分享】压电性、热释电性、铁电性与点群的关系
已经有8人回复
- 【讨论】[有奖活动]各种计算方法的优点
已经有4人回复
Gina88
主管区长
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
2楼2013-01-21 16:33:00
Gina88
主管区长
- 1ST强帖: 1
- 应助: 9 (幼儿园)
- 金币: 3328.8
- 散金: 124
- 红花: 8
- 帖子: 395
- 在线: 113小时
- 虫号: 1047684
- 注册: 2010-06-26
- 专业: 理论和计算化学
|
算硬度： PRL,2003,015502,"Hardness of covalent crystals" 中计算硬度的方法需要计算材料的平均带隙Eg。这一数值远大于第一性原理计算得到的禁带宽度。请问这个数值是如何计算得到的，请给出具体过程，谢谢。 回帖： 第一性原理计算得到的禁带宽度本来就比实验值偏小 方法挺多，加GW修正，或者用杂化密度泛函 google里面搜索GW + band gap,以及hybrid DFT + band gap.很多很多文献 回帖： 带隙偏低的原因在于，基于密度泛函理论的第一性原理方法得到的带隙和实验测得的带隙定义不同；前者忽略了电子间的作用，对于一个有N个电子的体系，其单电子近似的带隙为Eg=E(N,N+1)-E(N,N),而严格的带隙应该是Eg=E(N+1,N+1)-E(N,N)；因此密度泛函方法低估带隙。克服的方法如下几种：准粒子GW方法进行修正；或者采用杂化泛函，用哈特利-福克的自洽场计算，但是这两种方法据说都很耗时间，实际上是设备流，没有好的设备就悲剧了。另外，DFT+U也可以提高带隙，但这中修真仅限于强关联电子，例如f带和d带电子，而且我个人认为有凑数“嫌疑”，因为U的值是和化学环境有关的，某种程度上讲U是一个经验参数。 如果想要准的带隙，可以考虑GW和杂化泛函方法，GW我知道的ABINIT可以实现，PWSCF貌似也可以；ELK里面貌似可以试验杂化泛函方法。 我没用过这两种方法，只是建议，未必对。 |
3楼2013-01-21 16:39:37
youqin5912
专家顾问
- 应助: 13 (小学生)
- 金币: 1183.4
- 散金: 785
- 红花: 15
- 帖子: 397
- 在线: 594.4小时
- 虫号: 1141030
- 注册: 2010-11-07
- 性别: GG
- 专业: 同步辐射技术及其应用
好贴啊，顶一下 4楼2014-08-14 23:29:47
obaica
版主
|
本帖内容被屏蔽 |
5楼2016-11-14 20:35:32
obaica
管理员
|
本帖内容被屏蔽 |
6楼2016-11-14 20:35:41
