| 查看: 2576 | 回复: 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一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算晶体的能带谱图发现了一个很奇怪的现象
已经有9人回复
- 【专家在线答疑贴之三】杂化泛函计算FAQ
已经有28人回复
- 【求助】各位大侠！小女子求助钙钛矿结构计算磁性的问题？
已经有36人回复
- 【求助】VASP计算体系性能没有f能带电子
已经有24人回复
- 【求助】castep 优化铁时 initial spin该如何设置
已经有26人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【原创】DFT计算体会
已经有20人回复
- 【其他】在线答疑:VASP的问题
已经有314人回复
- 【求助】vasp noncollinear计算能带的时候只能算spin-unpolarized吗
已经有12人回复
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
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
wutianhuan
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 6697.4
- 散金: 12
- 红花: 3
- 帖子: 1723
- 在线: 216.2小时
- 虫号: 284387
- 注册: 2006-10-09
- 性别: MM
- 专业: 金属材料的合金相、相变及
5楼2012-07-12 11:25:28
20