| 查看: 3809 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
dashuai1133新虫 (小有名气)
|
[求助]
phonopy生成声子谱的PRIMITIVE_AXIS和PRIMITIVE_AXIS设置 已有1人参与
|
||
|
用phonopy将所有的vasprun.xml结果汇总得到力常数矩阵已做完，生成FORCE_SETS，看手册上说运行phonopy -f band.conf就能出结果。但对于band.conf的具体参数不太理解，不会设置。请各位老师指点：PRIMITIVE_AXIS=？ PRIMITIVE_AXIS=?  在线等结果。 |
回复此楼» 猜你喜欢
309求调剂
已经有9人回复
-
环境科学与工程334分求调剂
已经有4人回复
-
085404求调剂，总分309，本科经历较为丰富
已经有6人回复
-
285求调剂
已经有5人回复
-
求化学调剂
已经有7人回复
-
材料工程专硕求调剂
已经有4人回复
-
317求调剂
已经有10人回复
-
318求调剂
已经有4人回复
-
化学308分调剂
已经有3人回复
-
一志愿河北工业大学0817化工278分求调剂
已经有18人回复
fyuewen
金虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 1620.4
- 散金: 151
- 红花: 10
- 帖子: 480
- 在线: 125.2小时
- 虫号: 1506790
- 注册: 2011-11-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
原题引用“用phonopy将所有的vasprun.xml结果汇总得到力常数矩阵已做完，生成FORCE_SETS，看手册上说运行phonopy -f band.conf就能出结果。但对于band.conf的具体参数不太理解，不会设置。请各位老师指点：PRIMITIVE_AXIS=？ PRIMITIVE_AXIS=?” PRIMITIVE_AXIS的定义，其实在官网上的settings中已经给出（虽然有一处作者有点笔误）。而且楼上有一位其实已经说的很清楚了。M_p是用来对unit cell转化成 primitive cell用的。 具体定义为（a_p, b_p, c_p) = (a_u, b_u, c_u)*M_p 此处p代表primitive cell, u 代表unit cell. |
21楼2016-06-09 16:46:09
dashuai1133
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 210.4
- 散金: 60
- 红花: 1
- 帖子: 96
- 在线: 95.2小时
- 虫号: 2200052
- 注册: 2012-12-22
- 专业: 理论和计算化学
2楼2014-12-24 15:16:45
dashuai1133
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 210.4
- 散金: 60
- 红花: 1
- 帖子: 96
- 在线: 95.2小时
- 虫号: 2200052
- 注册: 2012-12-22
- 专业: 理论和计算化学
3楼2014-12-24 16:19:27
dx0620
木虫 (正式写手)
- 应助: 14 (小学生)
- 金币: 4053.6
- 散金: 5
- 红花: 6
- 帖子: 393
- 在线: 910小时
- 虫号: 1303311
- 注册: 2011-05-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢交流！ 2014-12-25 20:02:09
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢交流！ 2014-12-25 20:02:09
|
你可以把这一行去掉，因为一般你计算的时候是不会选取原胞的，这个是原胞基矢的选项，所以不用加这一个参量。band.conf可以这样写： ATOM_NAME= DIM = 2 2 2 BAND = 0 0 0 0 0 0.5 -0.5 0 0.5 -0.5 0 0 -0.5 0.5 0 0 0.5 0 0 0.5 0.5 -0.5 0.5 0.5 BAND_POINTS = 101 BAND_LABELS = \Gamma Z T Y S X U R BAND确实可以选取跟castep一样的，或者自己选点也行 |
4楼2014-12-25 09:21:31
