应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 建立石墨烯（Graphene）的模型
51/1返回列表
查看: 16283  |  回复: 26
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价19次，作者souledge增加金币 14.6 个
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

souledge

专家顾问 (著名写手)

[资源] 建立石墨烯（Graphene）的模型

现在石墨烯很火，虽然本人不是做这个方向的，但是偶尔也尝试下。
但是惊奇的发现，一个建模问题居然难倒了N多人，还讨论来讨论去的，实在有点无语，事实上，只要对晶体学的知识比较熟悉，自己画那么一层六方框框出来，然后把六方坐标轴在六方晶系常作原点的地方试几下，这个问题就解决了。
好吧，我不得不多说一句，晶体学是建模的基本功，虽然Materials Studio或者其他工具可以让建模变得轻松惬意，但是没有这些基础的知识，谁也无法保证建出来的模型是中看还中用的～
直接说建模的结论了，石墨烯中的碳原子位置是与六方二硼化钛中的硼原子完全等价的。
在PWscf中建模时（也包括其他所有的第一性原理软件输入文件），只能针对14种布拉菲格子的原胞基矢来建立，而空间群的对称性都是再软件计算时进行分析的。
所以，如果写为某空间群形式，石墨烯中等价位置其实只有一个，但是作为一个原胞，就必须是两个，因为对称性没有空间群中那么完善，需要自己补充。
所以首先建立六方的坐标：
CODE:
ibrav = 4,
celldm(1) = a*sqrt(3),
celldm(3) = N,

其中a是石墨烯中两个相邻碳原子的距离，N可以设置的足够大。
而如果手动指定基矢，则是：
CODE:
ibrav = 0,
celldm(1) = a*sqrt(3),
CELL_PARAMETERS cubic
1.0  0.0  0.0
-0.5  sqrt(3)/2  0
0.0  0.0  N

第二个方法适用于所有的第一性原理软件输入文件建模。
然后是原子位置：
CODE:
ATOMIC_POSITIONS crystal
  C  2/3  1/3  1/2
  C  1/3  2/3  1/2

把C放再Z轴中间是为了在可视化软件中易于观察和调节角度（当然，Z轴位置在哪里都是同样的计算结果）。
如此，就可以得到简单的石墨烯原胞模型：

不放心可以看沿X和Y轴方向扩展的结构：

扩展一下，如果是Abinit，可以建模为：
CODE:
acell  2*a*sqrt(3)  N
rprim
1.0  0.0  0.0
-0.5  sqrt(3)/2  0
0.0  0.0  1.0
natom  2
ntypat  2
znucl  6
typat  1  1
xred
  1/3  2/3  1/2
  2/3  1/3  1/2

而ELK，则可以是:
CODE:
scale
  a*sqrt(3)
avec
  1.0  0.0  0.0
  -0.5  sqrt(3)/2  0
   0.0  0.0  N
atoms
1
'C.in'
2
0.3333333333  0.6666666667  0.0000000000  0.0000000000  0.0000000000  0.0000000000
0.6666666667  0.3333333333  0.0000000000  0.0000000000  0.0000000000  0.0000000000

VASP的POSCAR：
CODE:
Graphene
  a*sqrt(3)
        1.0  0.0  0.0
        -0.5  sqrt(3)/2  0
         0.0  0.0  N
  2
Direct
   1/3  2/3  1/2
   2/3  1/3  1/2

最后是CASTEP的*.cell文件：
CODE:
%BLOCK LATTICE_CART
  a*sqrt(3)  0.0  0.0
  -a*sqrt(3)/2  a*3/2  0
  0.0  0.0  N
%ENDBLOCK LATTIC_CART
%BLOCK POSITIONS_FRAC
  C  1/3  2/3  1/2
  C  2/3  1/3  1/2
%ENDBLOCK POSITIONS_FRAC

注：除Abinit外，其他软件不具备计算除法和开方的功能，所以要把分数和无理数替换为有效数字稍微多一点的小数形式
回复此楼

» 收录本帖的淘帖专辑推荐

锂电收藏 石墨烯小百科 第一性原理 第一性原理相关文档
晶体模型计算分析 第一性原理计算经验 VASP资料 石墨烯模拟计算
科研软件小神器 量化软件 第一性 计算化学
科研 学习贴 田凯 爪机党
第一性原理计算 纳米材料 石墨烯 第一性原理
画图 每日新知识

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-01-25 00:47:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)

★★★★★ 五星级,优秀推荐

souledge: 取消置顶 2012-04-12 22:57:59
可视化的软件名是什么？
赞一下 回复此楼
10楼2012-04-12 21:07:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 27 个回答

yibabamsini

金虫 (小有名气)

★★★ 三星级,支持鼓励

虽然看不太懂，但觉得楼主的思维方式很正确，有道理，赞一个
赞一下 回复此楼
2楼2012-01-25 11:18:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独唱团

金虫 (正式写手)

★★★★★ 五星级,优秀推荐

学习了。。。。。
回复此楼
6楼2012-02-20 22:59:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flymice

金虫 (著名写手)

★★★★★ 五星级,优秀推荐

晶体基础不好的人伤不起呀。
赞一下 回复此楼
7楼2012-03-06 15:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 27 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:28 by G6APbkg8SA6w
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 09:17 by G6APbkg8SA6w
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭