24小时热门版块排行榜    

查看: 4803  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 刘仕晨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +5 scalable 2026-08-24 5/250 2026-08-24 17:03 by spiderzhuzhi
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见