| 查看: 2630 | 回复: 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 ) 然后，然后就什么都没有了。。。不知道哪位能解决这个问题啊？？？在此先谢过了~ |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算晶体的能带谱图发现了一个很奇怪的现象
已经有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
10