| 查看: 1865 | 回复: 4 | ||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 fzuvivi 的 13 个金币 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
|
||||
|
最近跑vasp，遇到了一些问题， （1）关于Si 体弹性模量的计算： 根据体模量的定义，我通过单体积下充分弛豫后得到一系列的不同晶格常数和对应能量值的关系（E-A），然后通过一个脚本murn.x利用Murnaghan方程来拟合得到最低能量对应的晶格常数大小（alat）以及体模量值（B0）；可是发现，虽然晶格常数可以与实验值吻合的很好（误差在1%左右）；但体模量却出现较的偏差（误差达10%以上）。后来也通过计算弹性常数来得到体模量，可是准确性没什么改观。。。请问，是不是有什么原因会很大程度的影响体模量的准确性。还是我拟合体模量的方法有问题，是否有其他的方法能得到较精确的模量值？ 附上Si的原子优化脚本： #!/bin/sh rm WAVECAR for i in 5.15 5.20 5.25 5.30 5.35 5.40 5.43 5.455.50 5.53 5.55 5.60 5.65 5.70 do cat >POSCAR< Si-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 Direct 0.0 0.0 0.0 0.25 0.25 0.25 ! echo "a = $i angstrom"; time vasp E=$(grep "TOTEN" OUTCAR |tail -1| awk '{printf "%12.6f \n", $5}') V=$(grep "volume" OUTCAR |tail -1| awk '{printf "%12.4f \n", $5}') echo $i $V $E >>EtVo.dat done INCAR： SYSTEM = Si-fcc ENCUT = 450 ISTART = 0; ICHARG = 2 ISMEAR = -5 IBRION = 2； NSW = 60 ISIF = 2 EDIFF = 1E-4； EDIFFG = -1E-3 PREC = Accurate （2）关于体材料的自旋极化设置问题： 计算Co基磁性材料时，在INCAR中设置了参数ISPIN=2；可是静态自洽scf计算、能带结构band计算等得到的OUTCAR中并未记录magnetization信息；只有在态密度dos计算的OUTCAR文件里有s、p、d 以及Total的magnetization？而且在同种Co基磁性体系的不同结构中，有的结构似乎没有自旋极化（ Total magnetization=0.0）；有的结构却有明显的自旋极化现象（Total magnetization=2.0）； 附上其中一个结构的算DOS的INCAR： SYSTEM = Co-work ENCUT = 450 ISTART = 1; ICHARG = 11 ISMEAR = 0; SIGMA = 0.1 GGA = 91; VOSKOWN =1 ISPIN = 2 RWIGS = 2.163 1.512 PREC = Accurate （3）关于LDA+U的设置： 在计算Co基材料时，由于有3d电子的存在，所以想对GGA势进行“+U”计算。可是不知道其中的U值和J值该怎么设置？（LDAU、LDAUTYPE、LDAUL、LDAUU、LDAUJ）。。vasp手册说要参考实验文献；具体要参考什么实验值？论坛上有人说带隙值，那如果是金属，怎么办？还有，如果没有实验文献值，那以什么来作为判据，选择正确的U和J值。 附：上次一位高手的解答，但是还是很疑惑 [ Last edited by fzuvivi on 2010-12-29 at 08:54 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 科研相关资料 |
» 猜你喜欢
有没有要一起买VASP版权的，可以6个人
已经有7人回复
-
我感觉当老师好累啊。。
已经有10人回复
-
求助！请问海外博后依托内地单位申请青基的优劣？
已经有9人回复
-
这种人痛苦吗
已经有14人回复
-
评正教授需要两个国家级项目，有人用子课题糊弄，结果在评审前资格公示时被举报拿下了
已经有31人回复
-
chemical science和advanced science哪个好一点
已经有6人回复
-
国内期刊审稿人数量
已经有5人回复
-
中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知
已经有18人回复
-
angew三个小修，返修过去已经20多天了没动静，有朋友遇到这种情况吗？
已经有5人回复
-
随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 还是决定了13年考研，希望得到点大家的经验，谢谢！
已经有229人回复
- 希望文章能顺利接受
已经有29人回复
- 求助：origin高斯拟合
已经有10人回复
- 拟合过过弹性模量的没？
已经有43人回复
- 提供一个 Birch-Murnaghan 拟合的算法（更新1）
已经有95人回复
- 丙烯酸酯弹性体
已经有23人回复
- 金属电沉积的极化曲线
已经有33人回复
- 为什么说，单个高分子链已具备了高弹性？
已经有5人回复
- 对数据进行曲线拟合后如何得到拟合误差
已经有30人回复
- 钢为什么会有弹性？
已经有5人回复
- 弹性常数计算出现诡异。。。。来看看喽
已经有5人回复
- 如何做出拟合曲线求值
已经有24人回复
- 求助椭偏拟合
已经有3人回复
- 求哪里可以查得到一些材料的杨氏模量、泊松比、热膨胀系数等参数呀？
已经有12人回复
- 石墨烯中的自旋轨道耦合
已经有7人回复
- 超精细相互作用下的自旋翻转
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
坐标川渝，非诚勿扰
+1/254
-
重庆大学化学化工学院，招聘电催化方向科研助理
+1/98
-
分离纯化和制备载样量的交流
+1/87
-
湖州第一人民医院/湖州师范学院医学院招聘合作博士后，助理研究员等优秀人才（全职）
+1/83
-
善准旋转型等离子体低温灰化仪两种装样方式煤粉低温灰化更均匀更充分
+2/52
-
Texas Tech University HSC 诚聘生物信息学、计算基因组学数据分析软件开发博士后数名
+1/37
-
放下难，拾起更难，知难而进，方得始终
+1/33
-
江苏大学程明教授课题组——招收钙钛矿太阳能电池领域博士生
+1/33
-
西北农林科技大学 生物质功能高分子材料课题组 招收2025申请考核制博士生
+1/31
-
中国人民大学化学与生命资源学院博士后招聘
+2/22
-
求职帖子，各位虫友
+1/17
-
伦敦大学学院（UCL）诚招手术机器人方向博士生（机械/电子/自动化/控制）
+1/7
-
硅碳软包电极拆解
+1/6
-
燕山大学高压科学中心李鹏教授课题组招收博士后及博士生
+3/6
-
北京化工大学程元徽教授团队博士研究生招生
+1/4
-
俄亥俄大学化学与生物化学系 招收全奖博士研究生
+1/3
-
韩国科学技术院KAIST招聘博后一名
+1/3
-
山东大学集成电路学院太赫兹团队招聘青年教师和博士后（长期有效）
+1/3
-
西班牙ICMAB-UAB招收2025博士生
+1/3
-
中国科学技术大学国家级青年人才招聘材料方向博士/博士后/特聘副研究员
+1/1
+1/254
2楼2010-12-29 08:56:47
