应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp5.2光学性质问题
192/2返回列表上一页12
查看: 4455  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

laplacetalor

银虫 (小有名气)
引用回帖:
2415845楼: Originally posted by 707462558 at 2012-12-19 11:54:22
楼主你好  我想问一下  计算光学性质的流程是什么样的呀  也是 优化---i静态自洽---非自洽吗？  输入文件依然是那四个吗  有什么改变呀...

流程就是那样。
静态计算的时候增加NBANDS到优化时的3倍， 并且保留算出来的chgcar和wavecar
非自洽的时候使用刚才提到的2个文件，修改incar如下
ISTART=1 ICHARG=11 LOPTICS=.T. NEDOS=3000 NBANDS和静态时一样
赞一下 回复此楼
11楼2012-12-19 13:07:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

707462558

银虫 (著名写手)
引用回帖:
2415854楼: Originally posted by laplacetalor at 2012-12-19 13:07:39
流程就是那样。
静态计算的时候增加NBANDS到优化时的3倍， 并且保留算出来的chgcar和wavecar
非自洽的时候使用刚才提到的2个文件，修改incar如下
ISTART=1 ICHARG=11 LOPTICS=.T. NEDOS=3000 NBANDS和静态时一样...

我每一步都没有设置NBAND呀  只有dos时候设NEDOS=2000  那这个设多少呀
赞一下 回复此楼
无论生活和你开了多大的玩笑我们都得继续生活
12楼2012-12-19 15:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)
★
laplacetalor(金币+3): 谢谢参与
引用回帖:
2415854楼: Originally posted by laplacetalor at 2012-12-19 13:07:39
流程就是那样。
静态计算的时候增加NBANDS到优化时的3倍， 并且保留算出来的chgcar和wavecar
非自洽的时候使用刚才提到的2个文件，修改incar如下
ISTART=1 ICHARG=11 LOPTICS=.T. NEDOS=3000 NBANDS和静态时一样...

非自洽时只要wavecar吧（在有的帖子里看到过），所以ICHARG=11或2都可以吧？NBANDS在静态计算要加到优化时的2-3倍？优化和静态都取默认（不设置NBANDS），而在非自洽时再加上NBANDS（等于scf时的2-3倍）就行了吧？有经验的同学给确认一下，谢谢

[ Last edited by bjwang on 2012-12-20 at 11:06 ]
赞一下 回复此楼
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
13楼2012-12-20 10:56:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laplacetalor

银虫 (小有名气)
引用回帖:
2415863楼: Originally posted by 707462558 at 2012-12-19 15:39:24
我每一步都没有设置NBAND呀  只有dos时候设NEDOS=2000  那这个设多少呀...

vasp manual说是默认值的3倍，一般优化就用默认值，去outcar里面找找有的。
赞一下 回复此楼
14楼2012-12-21 08:58:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laplacetalor

银虫 (小有名气)
引用回帖:
2415911楼: Originally posted by bjwang at 2012-12-20 10:56:17
非自洽时只要wavecar吧（在有的帖子里看到过），所以ICHARG=11或2都可以吧？NBANDS在静态计算要加到优化时的2-3倍？优化和静态都取默认（不设置NBANDS），而在非自洽时再加上NBANDS（等于scf时的2-3倍）就行了 ...

非自洽要求保持电荷密度不变，应该要chgcar吧，否则怎么做到不变。其实wavecar是不需要的，只是计算时间长一点而已。
静态是如果用默认NBANDS，而非自洽用3倍NBANDS会导致波函数不一样，这样非自洽时就可能读取波函数错误，虽然也可以算，但浪费时间要随机生成波函数，对于有些体系会不收敛。
赞一下 回复此楼
15楼2012-12-21 09:02:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjwang

铁杆木虫 (正式写手)
引用回帖:
2415966楼: Originally posted by laplacetalor at 2012-12-21 09:02:43
非自洽要求保持电荷密度不变，应该要chgcar吧，否则怎么做到不变。其实wavecar是不需要的，只是计算时间长一点而已。
静态是如果用默认NBANDS，而非自洽用3倍NBANDS会导致波函数不一样，这样非自洽时就可能读取波 ...

我是看到这个杂化泛函计算光学性质的帖子http://muchong.com/bbs/viewthread.php?tid=3753879  说：后面的一直是用前面的wavecar (始终ICHARG = 2)；并且，那里DFT时没有加NBANDS(应该用的默认)，在HSE06计算时才加NBANDS

       而在这个帖子了：http://cms.mpi.univie.ac.at/vasp ... topic.php?4.6633.20说“WAVECAR only. CHGCAR is not relevant for restart of HSE calculations.”
上面这些都是对杂化泛函计算光学有效，对DFT计算光学性质无效吗？
在这个回帖中（http://muchong.com/bbs/viewthread.php?tid=3753879&page=6）强调：“ICHARG一定不要设置成11”。
   
     我验证了一个体系，利用在SCF计算时不改变NBANDS值，在非自洽（optic）计算时增加NABANDS的数值为SCF时的2-3倍,且保持ICHARGE=2，结果显示和文献符合的好些（利用楼上那个方法的epsilon2多了小峰）。

    不过感觉楼上的说的也有道理，哪位提供一下权威说法，谢谢！

附：
Dear professor:
>         I am a VASP user. I want to calculate band structure and optical property with HSE06.I have read the post about "Band Structure Calculation With HSE06" in VASP forum.But now I also have several question about "admin's answer".http://cms.mpi.univie.ac.at/vasp ... topic.php?4.6633.10
>        Now my questions are as follows:
>        1.Dose "-1) First perform " for HSE scf need the WAVECAR of DFT calculation? In this step,how to set the parameter ICHARG?

Doing DFT calculations first is always a good idea for HSE.
ICHARG does not matter when you start a HSE calculation (it is disregarded).

>        2. The "-5) Perform a second VASP run: " is  A Hartee-Fock calculation ? Dose this step for band need  the  WAVECAR and CHGCAR file of prevous HSE06 calculation ? How to set the
> parameters ISTARTand ICHARG?
>WAVECAR only. CHGCAR is not relevant for restart of HSE calculations.

Georg Kresse

[ Last edited by bjwang on 2012-12-22 at 10:40 ]
赞一下 回复此楼
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
16楼2012-12-21 17:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfeng513

铜虫 (初入文坛)
★
laplacetalor(金币+3): 谢谢参与
请问DFT光学性质计算时，究竟是自洽还是非自洽计算啊？求求大侠提点提点
赞一下 回复此楼
17楼2013-09-28 19:35:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苦行客

木虫 (正式写手)
★ ★ ★
laplacetalor(金币+3): 谢谢参与
franch: 金币+2, 谢谢回帖交流 2013-10-15 09:37:16
引用回帖:
17楼: Originally posted by yangfeng513 at 2013-09-28 19:35:05
请问DFT光学性质计算时，究竟是自洽还是非自洽计算啊？求求大侠提点提点

非自洽
回复此楼
行有不得，反求诸己～
18楼2013-10-14 10:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaiyong2013

新虫 (初入文坛)
★
laplacetalor(金币+3): 谢谢参与
引用回帖:
10楼: Originally posted by 707462558 at 2012-12-19 11:54:22
楼主你好  我想问一下  计算光学性质的流程是什么样的呀  也是 优化---i静态自洽---非自洽吗？  输入文件依然是那四个吗  有什么改变呀...

把你的计算经验分享一下好吗，我是新手，谢谢
赞一下 回复此楼
学无止境
19楼2014-12-29 10:52:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 laplacetalor 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 05:00 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:57 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
关闭
请填处理意见
关闭