版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (2590)
>虫友互识 (120)
>基金申请 (93)
>论文投稿 (49)
>导师招生 (31)
>考博 (25)
>教师之家 (12)
>硕博家园 (12)
>找工作 (12)
>论文道贺祈福 (11)
>公派出国 (10)
>博后之家 (8)
>文献求助 (8)
>考研 (8)
>休闲灌水 (7)
>留学DIY (1)

返回列表

【奖励】本帖被评价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外，其他软件不具备计算除法和开方的功能，所以要把分数和无理数替换为有效数字稍微多一点的小数形式

回复此楼

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

» 猜你喜欢

2026年叶企孙基金已经有7人回复
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？已经有6人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有5人回复
有多少人是今天查系统知道结果的？已经有18人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有4人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有5人回复
国自然评审意见已经有15人回复
为什么资助数各大高校都创新高，自己申请怎么就这么难已经有13人回复
国自然面上复盘~欢迎讨论已经有14人回复

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

» 本主题相关价值贴推荐，对您同样有帮助:

吸附计算可否人为添加化学键？已经有10人回复
VASP 计算石墨烯与石墨烯带的区别已经有11人回复
在石墨烯表面吸附水分子MS该如何建模已经有7人回复
ACS Nano最新研究进展：氧化石墨烯的酸性起源、水中的不稳定性以及动力学结构模型已经有305人回复
分享：超强吸附材料--氧化石墨烯的文章（Journal of hazardous materials）已经有225人回复
关于Graphene, graphene oxide (GO) and reduced graphene oxide (rGO) 的区别？已经有7人回复
石墨烯中反π键轨道已经有4人回复
求助石墨烯纳米带的建模~ 已经有19人回复
【求助】请教MS建模，怎样在一种晶体表面上加上一层石墨烯？已经有12人回复
【求助】石墨烯的MS建模已经有13人回复
【求助】怎样用MS构建石墨烯模型？已经有12人回复
【求助成功】求助如何建立石墨烯模型及酰胺化学键模型已经有5人回复
【求助】掺氮石墨烯的紫外可见吸收已经有8人回复
【求助】石墨烯吸附气体分子的优化问题已经有23人回复
【交流】石墨与氢气在高温下是否会发生反应已经有3人回复
【讨论】求助：：NOTEEXPRESS 直接连线维普数据库下载全文已经有3人回复

1楼 2012-01-25 00:47:12

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

yolanda1234

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 3115.3
帖子: 169
在线: 179小时
虫号: 1669070

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

晶体学没有一点基础的人，伤不起，该学习什么能够入门的？楼主建议一个~~~

赞一下

回复此楼

15楼2013-01-14 17:02:06

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

查看全部 27 个回答

yibabamsini

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1269
帖子: 94
在线: 53.5小时
虫号: 1134428

★★★ 三星级,支持鼓励

虽然看不太懂，但觉得楼主的思维方式很正确，有道理，赞一个

赞一下

回复此楼

2楼2012-01-25 11:18:14

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

独唱团

金虫 (正式写手)

应助: 14 (小学生)
金币: 483
帖子: 729
在线: 142小时
虫号: 1308772

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

学习了。。。。。

回复此楼

6楼2012-02-20 22:59:29

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

flymice

金虫 (著名写手)

应助: 20 (小学生)
金币: 1166.7
帖子: 2592
在线: 498.9小时
虫号: 1053017

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

晶体基础不好的人伤不起呀。

赞一下

回复此楼

7楼2012-03-06 15:51:51

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

查看全部 27 个回答

☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级

普通表情龙兔虎猫高级回复 (可上传附件)

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	ASdOkHsho7FD 2026-08-28	5/250	2026-08-29 06:21 by gy1nBQXYQJqL
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13	晴天加油 2026-08-26	14/700	2026-08-28 18:36 by huagongfeihu
[基金申请] 基金不中，共勉 +11	eulota 2026-08-26	11/550	2026-08-28 14:22 by 火星超人xi
[基金申请] 怎么查啊 +6	huang1991js 2026-08-26	6/300	2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5	yuleib84 2026-08-26	6/300	2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5	ssyjh 2026-08-27	5/250	2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3	myrtle 2026-08-26	3/150	2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7	process2012 2026-08-23	10/500	2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9	炎黄贵胄 2026-08-22	10/500	2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25	木水思豆 2026-08-25	28/1400	2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9	trojank 2026-08-26	9/450	2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4	bjdxyxy 2026-08-26	4/200	2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7	一怀馨秋 2026-08-26	7/350	2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +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	beefly 2026-08-26	8/400	2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6	chengyan1220 2026-08-23	6/300	2026-08-25 19:45 by zfd97