24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qingeraaa

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 10.7
散金: 18
红花: 5
帖子: 219
在线: 195.9小时
虫号: 1503611
注册: 2011-11-21
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] ibrav=0找不到全部对称性已有1人参与

我用QE计算石墨烯的声子谱，需要找到Gamma点的全部对称性。但是用ibrav=0时，不论怎么调整晶格参数，都没有成功，还希望高人指点一二。
&SYSTEM
ibrav = 0,
celldm(1)=4.6478
nat = 2,
ntyp = 1,
ecutwfc = 100.D0 ,
CELL_PARAMETERS
1.00000000    0.00000000       0.00000000
-0.5000000    0.866025000    0.00000000
0.00000000    0.00000000       4.0000000
ATOMIC_POSITIONS crystal
C    0.333333333 0.66666667 0.5000000000
C    0.666666667 0.33333333 0.5000000000
Result： 8 Sym.Ops. (with inversion)  --问题所在

&SYSTEM
   ibrav = 4,
   celldm(1) = 4.6478,
   celldm(3) = 4.0,
   nat = 2,
   ntyp = 1,
   ecutwfc = 100.D0 ,
ATOMIC_POSITIONS crystal
C    0.333333333 0.66666667 0.5000000000
  C    0.666666667 0.333333333 0.5000000000
Result：24 Sym.Ops. (with inversion)  ---ibrav=4时可以找到，而且对应点群为D6h。

因为有需要要用ibrav=0，所以还希望大家帮帮忙。我把ibrav=0的输入文件附上

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : ene.rx.in

2014-11-21 22:25:13, 1.01 K

» 本帖已获得的红花（最新10朵）

qingeraaa

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

对于未来，我一点也不不担心，因为我相信，时光会把我变得更好！

1楼 2014-11-21 22:22:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Dirac12345

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 12
散金: 19
红花: 8
帖子: 88
在线: 64.4小时
虫号: 2638624
注册: 2013-09-05

引用回帖:

12楼: Originally posted by boya_tech at 2014-11-22 18:33:27
1 Bohr=0.529177A,这个是hexgonal，直接将a换算即可。
alat是默认的，可以不设置。当然也可以设置为 cell_para alat...

你好。。
这两天一直在跟一个虫子讨论 ibrav=0 celldm(1)的取值问题
据我现在了解： ibrav=0时，若不设置celldm(1)，则cell_para 选项卡的晶格参数可以由 A来表示
若设置 celldm(1), 则它由缩放系数来决定。 VASP里缩放系数一般是1，所以 celldm(1)=1/0.529. cell_para 直接用VASP里 POSCAR 的分数坐标来表示

但是这里我又不太明白了，究竟 ibrav=0 时， celldm(1) 取值或不取值应该分别如何决定呢？以及后面的 cell_para 参数决定等问题。

你提到 hexgonal 可以直接将a换算，与cubic有何不同吗？
但是我看说明书上写 cell_para 后面得加 { hexgonal 或 cubic}，也又说要加 {alat，等等的}。

我现在已经为 ibrav=0时 celldm（1）以及 cell_para 的取值快研究一周了。

能麻烦你再给我解释下这块的关系吗？我在网上也查了好多也还是没弄明白。
如果可以的话能加下你Q吗？谢谢
我的事 838839057