【答案】应助回帖

11楼2013-12-25 08:43:40

xuejing141

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2478.5
帖子: 57
在线: 283.7小时
虫号: 1979801
注册: 2012-09-06
专业: 凝聚态物性 II ：电子结构

引用回帖:

10楼: Originally posted by jpchou at 2013-12-24 21:00:01
http://www.ncts.ncku.edu.tw/phys/cmr/page/file/tool/band-procar.f

你可以試試這個
將 PROCAR 轉成 band structure

我试了一下这个程序，输出文件band-dn.dat和band-up.dat都是空的为什么呀，谢谢了

12楼2013-12-25 10:09:02

xuejing141

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2478.5
帖子: 57
在线: 283.7小时
虫号: 1979801
注册: 2012-09-06
专业: 凝聚态物性 II ：电子结构

引用回帖:

10楼: Originally posted by jpchou at 2013-12-24 21:00:01
http://www.ncts.ncku.edu.tw/phys/cmr/page/file/tool/band-procar.f

你可以試試這個
將 PROCAR 轉成 band structure

http://muchong.com/bbs/viewthread.php?tid=2158806&fpage=1 就是这儿说的问题，但是我不知道应该怎么样修改这个脚本才能解决数组越界问题，谢谢你啦

13楼2013-12-25 11:17:05

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

parameter (nbd = 200)
   parameter (nkd = 500)
   parameter (nxd = 300)
   parameter (natmd = 20)

第二~五行
這應該就是要修改的地方
nbd 是 number of band
nkd 是 number of kpoints
nxd 是 number of atom
natmd 是 number of type of atom
改成你要的
應該就可以用了

14楼2013-12-25 14:16:28

xuejing141

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2478.5
帖子: 57
在线: 283.7小时
虫号: 1979801
注册: 2012-09-06
专业: 凝聚态物性 II ：电子结构

引用回帖:

14楼: Originally posted by jpchou at 2013-12-25 14:16:28
parameter (nbd = 200)
   parameter (nkd = 500)
   parameter (nxd = 300)
   parameter (natmd = 20)

第二~五行
這應該就是要修改的地方
nbd 是 number of band
nkd 是 number of kpoints
nx ...

按你说的做了，还是不行，我能加你QQ吗？谢谢了

15楼2013-12-25 16:10:56

甲斐之虎

银虫 (正式写手)

应助: 60 (初中生)
金币: 302.3
散金: 220
红花: 7
帖子: 970
在线: 465小时
虫号: 2273045
注册: 2013-02-03
性别: GG
专业: 凝聚态物性 II ：电子结构

围观，学习下

别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。

16楼2013-12-25 19:36:34

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

引用回帖:

15楼: Originally posted by xuejing141 at 2013-12-25 16:10:56
按你说的做了，还是不行，我能加你QQ吗？谢谢了...

我~~~我是台灣人
沒有在用QQ ....
你執行過程應該會有�?息
貼出?砜纯�
這樣比較能知道你需要改什麼

17楼2013-12-25 23:08:52

sxm2012

木虫 (小有名气)

应助: 24 (小学生)
金币: 3052.9
红花: 5
帖子: 248
在线: 496.9小时
虫号: 1765567
注册: 2012-04-19
专业: 理论和计算化学

引用回帖:

4楼: Originally posted by xuejing141 at 2013-12-23 08:40:52
我会把他们分开了，下面的这两行字是上面那个图的说明，我不知道那个0.054e/A怎么来的，谢谢了

111.jpg
...

请问是怎么讲它们分开的啊？

18楼2014-02-19 20:06:34

前田庆次

木虫 (著名写手)

应助: 202 (大学生)
金币: 394.2
散金: 1140
红花: 87
帖子: 1344
在线: 455.2小时
虫号: 1985685
注册: 2012-09-08
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by jpchou at 2013-12-21 20:06:34
寫程式

在CHGCAR內
前面會先寫出結構資�? (就像POSCAR一樣)
接著是 number of grids (NX, NY, NZ)
有了這個 grid 就可以讀接下?淼� charge
在 spin-polarize 計算的 CHGCAR
會先寫出 spin up + spin down 的 ...

慶次郎前?硗針牵�
晚生最近一直在考慮怎麼做自旋極化電荷密度，看了jpchou前輩的回帖，又看了vasp的manmual有這樣一段：
For spinpolarized calculations, two sets of data can be found in the CHGCAR file. The first set contains the total charge density (spin up plus spin down), the second one the magnetization density (spin up minus spin down). For non collinear calculations the CHGCAR file contains the total charge density and the magnetisation density in the x, y and z direction in this order.
明白（阿爾法-spin） - （�?塔-spin）那種自旋電荷密度怎麼做了……
感謝周前輩一百遍，一百遍

长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。

19楼2014-03-26 11:06:57

前田庆次

木虫 (著名写手)

应助: 202 (大学生)
金币: 394.2
散金: 1140
红花: 87
帖子: 1344
在线: 455.2小时
虫号: 1985685
注册: 2012-09-08
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。

20楼2014-03-26 11:07:45