| 查看: 5300 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
donggaomin金虫 (正式写手)
|
[求助]
请计算过加静水压下体系的性质的高手帮忙看看-----加压下的压电系数
|
|||
|
祝各位五一节快乐！ 下面是我计算压电系数的步骤，请各位大侠看一下（用下面的步骤重复吴志刚的那篇PRL时，得出的结果和文献的变化趋势是一样的，但是e15的大小却约等于文献结果的一半，我分别用LDA和GGA算过，但结果差不多，也改变过精度，还是不行，所以我怀疑我的计算过程有问题，所以请各位帮忙看一下）： 文献中的结果：  我计算的结果（纵坐标的负值不用管，只看大小，因为这只是与相对的方向有关）：  （下面以我计算静水压力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计算过，结果也没变多少） 各位高手，您看上面的步骤中我有没有用错参数的地方？ [ Last edited by donggaomin on 2011-5-1 at 10:52 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 凝聚态物理 | soft study |
» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
急！！求助，vasp计算中 spin up/spin down 电荷密度图怎么画
已经有11人回复- 卡方分布、t分布及F分布的性质及应用
已经有71人回复
- CASTEP计算出现你们经常出现的问题了，过来瞧瞧。
已经有11人回复
- 求过柱子的高手帮忙
已经有37人回复
- 请大家帮忙看看，外导这样的回复是婉拒吗？
已经有7人回复
- 金属配合物稳定常数分子模拟计算
已经有7人回复
- 帮忙给看看选取电子衍射图，给指点指点，是多晶还是单晶，怎么跟文献不太一样
已经有12人回复
- 德国老板回信了，请各位帮忙看看，这种有没有戏？
已经有7人回复
- 加溶剂后计算，不收敛，请高手指点
已经有5人回复
- 做加压试验时，反应时间起始点应该怎么算？
已经有4人回复
- 冷等静压成型时，升压阶段有必要分段加压吗？
已经有11人回复
- 看看别人仅用VASP做DOS计算发好文章（两篇PRL）
已经有30人回复
- 从cif文件直接导入的晶体可以用于直接计算么
已经有4人回复
- 请帮忙看看这个服务器配置单，主要用于Materials-studio计算
已经有10人回复
2楼2011-05-02 02:31:14