| 查看: 2354 | 回复: 4 | ||||
[交流]
【求助】关于拟合得到体弹性模量&自旋极化&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 ] |
» 收录本帖的淘帖专辑推荐
科研相关资料 |
» 猜你喜欢
一志愿西北工业大学289 085602
已经有27人回复
281求调剂
已经有7人回复
284求调剂
已经有11人回复
0858求调剂
已经有8人回复
287求调剂
已经有21人回复
求调剂
已经有5人回复
人工智能320调剂08工类还有机会吗
已经有8人回复
290求调剂
已经有3人回复
272分材料子求调剂
已经有25人回复
332求调剂
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
还是决定了13年考研,希望得到点大家的经验,谢谢!
已经有229人回复
希望文章能顺利接受
已经有29人回复
求助:origin高斯拟合
已经有10人回复
拟合过过弹性模量的没?
已经有43人回复
提供一个 Birch-Murnaghan 拟合的算法(更新1)
已经有95人回复
丙烯酸酯弹性体
已经有23人回复
金属电沉积的极化曲线
已经有33人回复
为什么说,单个高分子链已具备了高弹性?
已经有5人回复
对数据进行曲线拟合后如何得到拟合误差
已经有30人回复
钢为什么会有弹性?
已经有5人回复
弹性常数计算出现诡异。。。。来看看喽
已经有5人回复
如何做出拟合曲线求值
已经有24人回复
求助椭偏拟合
已经有3人回复
石墨烯中的自旋轨道耦合
已经有7人回复
超精细相互作用下的自旋翻转
已经有7人回复
» 抢金币啦!回帖就可以得到:
江西科师大招收与湖南大学联合培养学硕1名
+1/488
华北理工大学 矿业工程学院 矿业工程招收调剂研究生4名(学硕)
+5/25
化学、环境工程,上车就走
+1/24
杭州医学院张逢质教授课题组招收药学硕士研究生若干
+1/19
2026广东医科大学生物医学工程学硕、电子信息专硕接收调剂
+1/14
大连大学-贵州省煤炭洁净利用重点实验室联合培养研究生 化学5人+环境工程7人
+1/14
大连大学 化学专业 招收26年硕士研究生调剂
+1/11
大湾区大学李红庚课题组招聘研究助理
+1/7
水产方向招收调剂,090800和095134的可以看过来。
+1/6
江苏海洋大学+环境与化学工程学院+课题组招调剂生+第一批调剂
+2/6
齐齐哈尔大学李莉课题组诚招2026级考研调剂生(学硕和专硕)
+1/6
江西水利电力大学土木水利第二轮调剂系统已开
+1/6
26申博,ai+材料方向,求课题组推荐
+1/5
武汉工程大学化学与环境工程学院招收调剂学生
+1/5
济南大学天然药物化学生物学省重点实验室招收生物学、制药工程学专硕调剂
+1/4
徐工-环境工程学院-招收调剂硕士
+1/4
河北大学分析化学招收多名调剂考生
+1/3
上海应用技术大学化工学部招生调剂
+1/3
温州大学招2026年博士研究生(电催化领域)
+1/1
光学工程考生求调剂
+1/1
2楼2010-12-29 08:56:47
3楼2010-12-29 09:44:18
|
本帖内容被屏蔽 |
4楼2011-01-05 09:14:18
5楼2011-01-05 14:15:49













回复此楼
