小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
vasp能带计算求助
请问为什么在计算能带的时候会出现这样的结果？这是我计算能带时候的INCAR，能带图也在下面。各位好心人，谢谢啦
SYSTEM=FeSe2
ISTART=1;ICHARG=11
PREC=Accurate
#LORBIT=11
EDIFF=1E-4
NELM=200;NELMIN=4
# GGA=PBE0;VOSKOWN=1
#GGA=B3 GGA type
#ISMEAR=-5
#SIGMA=0.1
# ALGO = ALL ; TIME = 0.4
# AEXX = 0.5150 exact exchange contribution
# LHFCALC = .TRUE. ; HFSCREEN = 0.2
# ALGO = TIME = 0.4
# ALGO =D TIME = 0.4
# ENCUTFOCK = 0
# NKRED =2
LDAU=.TRUE.
LDAUTYPE=1
LDAUL=-1 -1
LDAUU=3.4 0
LDAUJ=0.5 0
LDAUPRINT=1
Ionic relaxation
EDIFFG = -0.01 ! stopping-criterion for IOM
NSW = 0 ! number of steps for IOM
IBRION = -1 ! ionic relax: 0-MD 1-quasi-New 2-CG
ISIF = 3 ! stress and relaxation
POTIM = 0.2 !
优化的时候是这样的INCAR
SYSTEM=FeSe2
ISTART=0;ICHARG=2
PREC=Accurate
#NBANDS=24
#LORBIT=11
EDIFF=1E-4
NELM=200;NELMIN=4
# GGA=PBE0;VOSKOWN=1
#GGA=B3 GGA type
#SIGMA=0.1
# ALGO = ALL ; TIME = 0.4
# AEXX = 0.5150 exact exchange contribution
# LHFCALC = .TRUE. ; HFSCREEN = 0.2
# ALGO = TIME = 0.4
# ALGO =D TIME = 0.4
# ENCUTFOCK = 0
# NKRED =2
LDAU=.TRUE.
LDAUTYPE=1
LDAUL=-1 -1
LDAUU=3.4 0
LDAUJ=0.5 0
LDAUPRINT=1
Ionic relaxation
EDIFFG = -0.01 ! stopping-criterion for IOM
NSW = 200 ! number of steps for IOM
IBRION = 2 ! ionic relax: 0-MD 1-quasi-New 2-CG
ISIF = 3 ! stress and relaxation
POTIM = 0.2 !
KPOINTS如下
计算能带的KPOINTS如下
k-points along high symmetry
0.000 0.000 0.000
0.000 0.000 0.500
0.000 0.000 0.500
-0.500 0.000 0.500
-0.500 0.000 0.500
-0.500 0.000 0.000
-0.500 0.000 0.000
-0.500 0.500 0.000
-0.500 0.500 0.000
0.000 0.500 0.000
0.000 0.500 0.000
0.000 0.500 0.500
有一点很奇怪，当把IALGO改为38时就没有问题，但是我觉得这个不具有说服力。&&Thank you very much！
自洽计算的INCAR
ISTART=0;ICHARG=2
PREC=Accurate
#NBANDS=24
#LORBIT=11
EDIFF=1E-4
NELM=200;NELMIN=4
# GGA=PBE0;VOSKOWN=1
#GGA=B3 GGA type
#SIGMA=0.1
# ALGO = ALL ; TIME = 0.4
# AEXX = 0.5150 exact exchange contribution
# LHFCALC = .TRUE. ; HFSCREEN = 0.2
# ALGO = TIME = 0.4
# ALGO =D TIME = 0.4
# ENCUTFOCK = 0
# NKRED =2
LDAU=.TRUE.
LDAUTYPE=1
LDAUL=-1 -1
LDAUU=3.4 0
LDAUJ=0.5 0
LDAUPRINT=1
Ionic relaxation
EDIFFG = -0.01 ! stopping-criterion for IOM
NSW = 0 ! number of steps for IOM
IBRION = -1 ! ionic relax: 0-MD 1-quasi-New 2-CG
ISIF = 3 ! stress and relaxation
POTIM = 0.2 !
优化的过程不就是自洽的过程吗？为什么还要跑一遍自洽。。。
进行态密度和能带计算的时候必须要用到前面一步自洽计算得到的CHG啊，是这样的吧。
是这样，但优化过程得到的CHG不能用吗？
我用的倒空间，所以是R，参考的MS。是错的吗？ 谢谢。
不能 优化的CHGCAR是不精确的，静态计算是为得到更好CHGCAR
自旋不是只有磁性材料才用到的吗？对这块我不怎么理解。还请指点！
我听别人说的 ，NBANDS是设置的OUTCAR里面的NBANDS就可以了，请问您是怎样根据不同情况设置的？
我看到你的能带图上面有锯齿形状，我以前也出现过.所以你看看静态中OUTCAR中NBANDS是多少，在INCAR中在设置大一些 就可以避免这个问题出现。你试试看。
这个我还真没考虑过 ，只是用OUTCAR中的值 ，也没改&&，谢谢 。
我试试 ，不过我觉得这个问题不大 ，因为他会按照参数要求来计算 ，既然是静态计算 ，而且NSW=0了已经 ，即使我设置ISIF为3，晶格参数也不会动的才是 。不知道我的想法是不是对的 。谢谢！
你试一下吧，我看你别的参数也没什么问题。我也用+U的算过，只有这个不同，呵呵
您好！我是一个刚刚接触VASP的菜鸟。刚刚拜读了您的应助。感觉到您对计算模拟这块儿挺有见解的。因为最近要算GGA+U，但是我对他又是几乎一窍不通。现在想麻烦您解释一下一些INCAR文件的设置问题，特别是
LDAU=.TRUE.
LDAUTYPE=1
LDAUL=-1 -1
LDAUU=3.4 0
LDAUJ=0.5 0
LDAUPRINT=1
这些具体指什么意思。如果有幸得到您的帮助，建不胜感激。您如果方便的话能加一下您的QQ吗？以方便以后向您请教。我的QQ号是.
这些参数的意思可以从vasp手册中查到，虽然都是英文的，但是多看几遍就可以看懂了，我对这些也只是一知半解的:cry:。
KPOINTS确实很让我头疼，我一般就是在MS里面导出来的，没有细想。大神，请问该怎么确定KPOINTS？
我以为大神会有什么秘籍呢:D:D是该好好看手册:work::work:
我不是大神，也没什么秘籍，就手册吧……
很详细，说的很有道理，很感谢。
最近没上小木虫，这个DFT+U的所有参数在手册都有介绍，你只要对应看一看就行。在2009版英文手册6.63节，96页。
首先感谢大侠能在百忙之中抽时间回复小弟的帖子。谢谢！DFT+U小弟已经通过虫友的帮忙和您说的VASP手册弄得稍微懂点儿皮毛了。不过老板又让搞LDA+DMFT，不知道大侠您对这个有研究吗？希望您能指点指点。由于发现2012年的VASP手册没有介绍有关的LDA+DMFT内容。您如果方便的话能不能给推荐几个参考手册。再次谢谢您！！
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
MS 里面的所有模块都是基于第一性原理来编写的么？ GULP模块是不是？
如题！& &大神们& &给点建议哈！
粗略的应用gulp拟合过势参数，当时的思路就是运用第一性原理计算的能量和gulp经典的势参数计算的能量比较
红花送上！&&大神！& &如果我用GULP模块来模拟离子在晶体中的扩散作用，可行么？&&望赐教！
你这离子扩散属于MD的范畴吧&&dl_poly，lammps可能比gulp更合适点
鲜花送上！&&大神！&&我想要得到离子在晶体中的扩散轨迹，&&用MS的那个模块计算出来的数据&&社会认可度高点哈！& & 望赐教！&&感谢先！
MS里面好像没有适合离子晶体扩散计算的模块，Gulp也可以，但好像J.D. Gale不太建议GULP来做MD计算，我的建议还会你用dl_poly来做离子晶体相关扩散的研究，用这款软件包来做的话，大家都是认可的。dl_poly也很好学
还有个问题&&希望大神赐教：
我在MS中用一个模块去论述一个问题& &再用另外一个模块去辅助证明一下这个问题
&&比如& &用CASTEP去算一下迁移势垒&&得到一个最佳路径&&再用Dmol3 去算一下MD证明最佳迁移路径& && &这样用两个模块在论文中可能使用？&&求大神指点+++
& & 鄙人金币虽少& &但必不吝啬！
如果算离子晶体中的迁移能垒的话，CASTEP算的不准确，文献中常见的CASTEP做过渡态一般都采用constrained minimization approach（很多文献你可以自己搜一下）。Dmol做MD可以，
你的想法不错，我目前没见过有你这样想法的论文。
看到过用VASP 计算能垒和AIMD来证明的文章，可能有异曲同工之妙。
但也可能会存在问题一定的问题，另外，不知道你为什么对MS这么情有独钟。
鲜花送上！&&非常感谢！&&科研经费有限& &只买了MS& &没有买VASP！
资源限制+++++
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研第一性计算原理_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
第一性计算原理
上传于||暂无简介
阅读已结束，如果下载本文需要使用0下载券
想免费下载更多文档？
定制HR最喜欢的简历
下载文档到电脑，查找使用更方便
还剩5页未读，继续阅读
定制HR最喜欢的简历
你可能喜欢小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
【原创】第一原理计算的一些心得
把自己关于第一原理的一些理解写出：
（1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质；
评述：K-S方程的计算水平达到了H-F水平，同时还考虑了电子间的交换关联作用。
（2）关于DFT中密度泛函的Functional，其实是交换关联泛函
包括LDA，GGA，杂化泛函等等
一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案；
GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BLYP泛函也属于GGA；
此外还有一些杂化泛函，B3LYP等。
（3）关于赝势
在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法，一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。
赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。
& & 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。
（4）关于收敛测试
（a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。
（b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。
（5）关于磁性
一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。
（6）关于几何优化
包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。
在PRL一篇文章中见到过只优化原子位置，晶格常数用实验值的例子（PRL 100, 08)）；也见到过晶格常数先优化，之后固定晶格常数优化原子位置的情况；更多的情况则是Full geometry optimization。
一般情况下，也有不优化几何结构直接计算电子结构的，但是对于缺陷形成能的计算则往往要优化。
（7）关于软件
软件大致分为基于平面波的软件，如CASTEP、PWSCF和ABINIT等等，计算量大概和体系原子数目的三次方相关；还有基于原子轨道线性组合的软件(LCAO)，比如openmx，siesta，dmol等，计算量和体系原子数目相关，一般可模拟较多原子数目的体系。
参考文献：
Designing meaningful density functional theory calculations in materials science—a primer
理解错误，呵呵，那想请教一下，PWSCF的USPP可以计算红外和拉曼光谱不？
这个要顶！有一处疑问和大家讨论：
侯老师的书上指出，这种全部优化，即晶格常数和原子位置同时优化，主要是用来估计的，计算得到较合理，但精度不高。不知怎样去理解？我想同时放开各种自由度应该得到更为可靠的结果。请各位指教！谢谢！
书在哪？可否传我一份？感激不尽啊
文章是侯老师推荐的，嘿嘿
NC：IR（frequency and intensity）+Raman（frequency and intensity）
US：IR（frequency and intensity）+Raman（frequency）
呵呵，，，顶一个！呵呵，，，顶一个！呵呵，，，顶一个！
想问一下优化的时候晶格常数和原子位置怎么设置才是一起优化，怎样才算单独优化啊，在优化以前需要怎样设置。我原来认为都是一起优化出来了，不可分开，谢谢高手指教。
wo ding a :D:D:D:D:D:D:D:D:D:D
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研}