| 查看: 2223 | 回复: 7 | ||
stevencong新虫 (初入文坛)
|
[求助]
~~有关磁性体系计算能带时遇到的问题~~
|
|
把全部金币贴进来了，求大牛们指导啊~~ 现在在做磁性纳米管的能带计算，但是不太清楚计算方法。我的计算过程是这样的： 1. 计算体系的Charge Density，输入文件如下（其中MAGMOM是从弛豫后的OUTCAR文件中抄下来的）： SYSTEM = Nanotube ENCUT = 800 ISTART = 0; ICHARG = 2 ISMEAR = -5 ISPIN = 2 MAGMOM = 0.03 0.03 -0.01 0.007 0.03 0.03 0.038 0.038 0.006 -0.01 0.039 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.039 0.039 0.007 -0.01 0.038 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.038 0.038 0.006 -0.01 0.038 0.039 PREC = Accurate 2. 再做BS计算，输入文件为： SYSTEM = Nanotube ENCUT = 800 ISTART = 1; ICHARG = 11 NBANDS = 120 ISPIN = 2 MAGMOM = 0.03 0.03 -0.01 0.007 0.03 0.03 0.038 0.038 0.006 -0.01 0.039 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.039 0.039 0.007 -0.01 0.038 0.038 0.03 0.03 -0.01 0.006 0.03 0.03 0.038 0.038 0.006 -0.01 0.038 0.039 PREC = Accurate 但是每次在计算BS时VASP刚开始运行就挂掉了，out文件输出： running on 16 nodes distr: one band on 1 nodes, 16 groups vasp.5.2.2 15Apr09 complex POSCAR found type information on POSCAR C POSCAR found : 1 types and 36 ions ----------------------------------------------------------------------------- | | | ADVICE TO THIS USER RUNNING 'VASP/VAMP' (HEAR YOUR MASTER'S VOICE ...): | | | | You have a (more or less) 'large supercell' and for larger cells | | it might be more efficient to use real space projection operators | | So try LREAL= Auto in the INCAR file. | | Mind: If you want to do a very accurate calculations keep the | | reciprocal projection scheme (i.e. LREAL=.FALSE.) | | | ----------------------------------------------------------------------------- LDA part: xc-table for Pade appr. of Perdew WARNING: stress and forces are not correct POSCAR, INCAR and KPOINTS ok, starting setup FFT: planning ...( 1 ) 然后，然后就什么都没有了。。。不知道哪位能解决这个问题啊？？？在此先谢过了~ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算晶体的能带谱图发现了一个很奇怪的现象
已经有9人回复
- 【专家在线答疑贴之三】杂化泛函计算FAQ
已经有28人回复
- 【求助】各位大侠！小女子求助钙钛矿结构计算磁性的问题？
已经有36人回复
- 【求助】VASP计算体系性能没有f能带电子
已经有24人回复
- 【求助】castep 优化铁时 initial spin该如何设置
已经有26人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【原创】DFT计算体会
已经有20人回复
- 【其他】在线答疑:VASP的问题
已经有314人回复
- 【求助】vasp noncollinear计算能带的时候只能算spin-unpolarized吗
已经有12人回复
mousekingadv
金虫 (小有名气)
- 应助: 51 (初中生)
- 金币: 1141.8
- 红花: 1
- 帖子: 140
- 在线: 28.7小时
- 虫号: 1721002
- 注册: 2012-03-27
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指教！ 2012-07-11 16:26:13
stevencong: 金币+20, ★★★很有帮助 2012-07-12 14:34:45
感谢参与，应助指数 +1
fzx2008: 金币+3, 谢谢指教！ 2012-07-11 16:26:13
stevencong: 金币+20, ★★★很有帮助 2012-07-12 14:34:45
|
这不是挂掉，这是条件太苛刻以致于一步都跑不出来…… ENCUT=800...手册上第一个例子就是讲的测试ENCUT吧...无论如何800太高了。 NBANDS并不需要手动提高，如果不关心高能空带的结构、也不需要计算光学性质的话。 另外KPOINTS不能过多。楼主是不是沿着高对称路径“老板，每条路径给我来50个K点”？先用较低的标准测试成功了，再在可接受的世界范围内提高精度。 LREAL=AUTO可以考虑写上。 |
2楼2012-07-11 15:48:06
stevencong
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 8
- 帖子: 18
- 在线: 8.4小时
- 虫号: 1732947
- 注册: 2012-04-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-07-11 16:23:42
stevencong
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 8
- 帖子: 18
- 在线: 8.4小时
- 虫号: 1732947
- 注册: 2012-04-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
结果又变成这样了： running on 8 nodes distr: one band on 1 nodes, 8 groups vasp.5.2.2 15Apr09 complex POSCAR found type information on POSCAR C POSCAR found : 1 types and 36 ions LDA part: xc-table for Pade appr. of Perdew WARNING: stress and forces are not correct POSCAR, INCAR and KPOINTS ok, starting setup WARNING: small aliasing (wrap around) errors must be expected FFT: planning ...( 1 ) reading WAVECAR WARNING: dimensions on CHGCAR file are different ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 ERROR: charge density could not be read from file CHGCAR for ICHARG>10 我两次计算所用的POSCAR应该一样啊，为什么Dimension会不一样呢？ |
4楼2012-07-12 00:19:55
wutianhuan
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 6697.4
- 散金: 12
- 红花: 3
- 帖子: 1723
- 在线: 216.2小时
- 虫号: 284387
- 注册: 2006-10-09
- 性别: MM
- 专业: 金属材料的合金相、相变及
5楼2012-07-12 11:25:28
stevencong
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 8
- 帖子: 18
- 在线: 8.4小时
- 虫号: 1732947
- 注册: 2012-04-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2012-07-12 14:36:18
wutianhuan
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 6697.4
- 散金: 12
- 红花: 3
- 帖子: 1723
- 在线: 216.2小时
- 虫号: 284387
- 注册: 2006-10-09
- 性别: MM
- 专业: 金属材料的合金相、相变及
7楼2012-07-13 08:33:05
【答案】应助回帖
★ ★
franch: 金币+2, 谢谢回帖交流。 2012-07-29 21:46:14
franch: 金币+2, 谢谢回帖交流。 2012-07-29 21:46:14
|
计算能带用ICHARG=11没有问题 能带的计算就是在自洽计算之后的到的电荷密度分布即为CHGCAR的前提下，进行费自洽计算的得到的 也就是说自洽计算得到CHGCAR，用做下一步计算 你是不是只是ICHGCAR=11，但是没有给CHGCAR，或是在计算上改变了什么参数，使得计算Bands时候，无法正确读取电荷分布，所以程序会报错“ERROR: charge density could not be read from file CHGCAR for ICHARG>10”，也就是说参数设置大于10，要求读取CHGCAR但是没有读到相应数据 |
8楼2012-07-29 21:13:46
