24小时热门版块排行榜    

查看: 4005  |  回复: 35

gemucai

木虫 (正式写手)

引用回帖:
20楼: Originally posted by 刘仕晨 at 2013-04-22 13:27:13
您好,我想问下,pwscf计算band时KPOINTS的高对称点的选择,一定要一个一个把所有的点都写出来么,能不能像vasp里面的那样,只给出高对称点和告诉它之间产生的点数,自动生成呢,比如像vasp里这样的:
k-points a ...

貌似有这种程序,但我是自己写的,这个程序自己写也不难
21楼2013-04-23 07:13:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
21楼: Originally posted by gemucai at 2013-04-23 07:13:53
貌似有这种程序,但我是自己写的,这个程序自己写也不难...

你是写出高对称点,然后按中间分多少点,然后自己算了,再写到in文件的么,我看到侯老师的博客中有这样一个fortran程序,http://valenhou.blog.edu.cn/2006/133292.html,我用的是这个,只是我不明白的是,那个计算的fcc的Cu的KPOINTS,那28个K点是怎么产生的,fcc的Brillouin Zones的高对称点好像不是那样的阿?还有为什么我用fortran脚本的时候都必须把里面注释的部分去掉才能成功编译出来,这就是为什么之前您给我发的band的脚本我处理不鸟的原因,而且我的后缀名必须是.f90,奇怪?谢谢。
22楼2013-04-23 15:18:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
22楼: Originally posted by 刘仕晨 at 2013-04-23 15:18:20
你是写出高对称点,然后按中间分多少点,然后自己算了,再写到in文件的么,我看到侯老师的博客中有这样一个fortran程序,http://valenhou.blog.edu.cn/2006/133292.html,我用的是这个,只是我不明白的是,那个计算 ...

我就是看不懂他写的代码,所以才自己编的,我觉得自己编比看别人写的代码要容易很多。不过,fcc的高对称点的确就是那些,可以参见http://en.wikipedia.org/wiki/Bri ... 81.29.2C_BCC.281.29

你提到的fortran有注释就无法编译的问题,可能是你用的编译器与我用的不一样,换个编译器就OK了。
23楼2013-04-23 16:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
23楼: Originally posted by gemucai at 2013-04-23 16:51:49
我就是看不懂他写的代码,所以才自己编的,我觉得自己编比看别人写的代码要容易很多。不过,fcc的高对称点的确就是那些,可以参见http://en.wikipedia.org/wiki/Brillouin_zone#Cubic_lattice_system_CUB.281.29.2 ...

不是,我的意思是,fcc的高对称点不应该是:
G 0 0 0
L 0.5 0.5 0.5
X 0.5 0 0.5
W 0.5 0.25 0.75
K 0.375 0.375 0.75
而那个计算的in文件里的28个点,有些是0 0 1fcc里没有这个高对称点吧,准确的说我就是 不明白,那28个点中的前11个点是怎么来的,是G点到那个点?谢谢。
24楼2013-04-23 17:55:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
24楼: Originally posted by 刘仕晨 at 2013-04-23 17:55:36
不是,我的意思是,fcc的高对称点不应该是:
G 0 0 0
L 0.5 0.5 0.5
X 0.5 0 0.5
W 0.5 0.25 0.75
K 0.375 0.375 0.75
而那个计算的in文件里的28个点,有些是0 0 1fcc里没有这个高对称点吧,准确的说我就是  ...

你如果感兴趣的话,可以做一个坐标变换,http://en.wikipedia.org/wiki/Brillouin_zone里面的fcc坐标轴变换到黄昆《固体物理》第179页的坐标轴上,这样你就知道001是哪个点了。
坐标轴选取得不同,坐标也会不一样
25楼2013-04-23 18:21:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
25楼: Originally posted by gemucai at 2013-04-23 18:21:25
你如果感兴趣的话,可以做一个坐标变换,http://en.wikipedia.org/wiki/Brillouin_zone里面的fcc坐标轴变换到黄昆《固体物理》第179页的坐标轴上,这样你就知道001是哪个点了。
坐标轴选取得不同,坐标也会不一样...

我不知道我理解的对不对阿,您说的坐标轴选取不同,是不是就是最后的单位不同阿,选择的pi/a或者2pi/a阿,如果是这样的话,那pwscf这个程序默认的到底是那个单位呢,或者我要怎么告诉它我用的是那个呢,不然计算不就不对了呢?谢谢!总是麻烦您,真的太感谢了。
26楼2013-04-26 12:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
26楼: Originally posted by 刘仕晨 at 2013-04-26 12:31:46
我不知道我理解的对不对阿,您说的坐标轴选取不同,是不是就是最后的单位不同阿,选择的pi/a或者2pi/a阿,如果是这样的话,那pwscf这个程序默认的到底是那个单位呢,或者我要怎么告诉它我用的是那个呢,不然计算不 ...

坐标轴不同,主要是指它们的方向不同,所以得出高对称点的坐标也就不一样了
27楼2013-04-26 14:45:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
27楼: Originally posted by gemucai at 2013-04-26 14:45:25
坐标轴不同,主要是指它们的方向不同,所以得出高对称点的坐标也就不一样了...

晕了,我明天再看看,那我计算的时候到底选那个坐标轴呢,两个都可以么,计算结果没有影响么?谢谢
28楼2013-04-26 22:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
28楼: Originally posted by 刘仕晨 at 2013-04-26 22:17:59
晕了,我明天再看看,那我计算的时候到底选那个坐标轴呢,两个都可以么,计算结果没有影响么?谢谢...

当然两个都可以,就看你喜欢用什么样的坐标系了
29楼2013-04-27 09:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
29楼: Originally posted by gemucai at 2013-04-27 09:04:46
当然两个都可以,就看你喜欢用什么样的坐标系了...

谢谢!
30楼2013-04-27 10:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘仕晨 的主题更新
信息提示
请填处理意见