| 查看: 1270 | 回复: 8 | |||
tfl03木虫 (职业作家)
|
[求助]
Lammps： 请帮忙建立一个 石墨 graphite 的输入晶格模型，发短息亦可
|
|
用下面的语句建立模型，老出问题，主要是开始原子个数都不对，谢谢！ units metal boundary p p p dimension 3 atom_style atomic #想建立石墨graphite的晶格，根据MS的数据，六方结构，a=b=2.46和c=6.8 Angstroms，a和b程120°，其余90° lattice custom 2.46 a1 1 0 0 a2 -0.5 0.866 0 a3 0 0 2.74 & basis 0 0 0 basis 0.333 0.667 0 & basis 0 0 0.5 basis 0.667 0.333 0.5 #region box block 0 1 0 1 0 1 #create_box 1 box #create_atoms 1 box #说有6个原子，丢了4个 region myr1 prism 0 1 0 1 0 1 0 0 0 create_box 1 myr1 create_atoms 1 region myr1 #同样说有6个原子，丢了4个 mass 1 12.0 velocity all create 5.0 87287 loop geom pair_style rebo pair_coeff * * CH.airebo C pair_modify shift yes neighbor 0.3 bin neigh_modify delay 0 every 1 check yes fix 1 all nve #fix 2 all enforce dump 1 all atom 100 dump.min thermo 100 run 1000 neigh_modify delay 0 every 1 check yes dump_modify 1 every 25 thermo 50 minimize 1.0e-26 0.00001 1000 10000 [ Last edited by tfl03 on 2012-6-24 at 14:59 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 吡啶型氮与石墨型氮有什么区别
已经有24人回复
- 谁能给我一些石墨烯/陶瓷块体复合材料的相关论文，谢谢！！！
已经有6人回复
- 怎样测量氧化石墨以及RGO的碳氧比C/O
已经有9人回复
- 一次型NiAl结构优化，晶格矢量变化非常大，求助~
已经有5人回复
- Accounts of Chemical Research最新综述：二维软材料---氧化石墨烯的另一面
已经有314人回复
- 氧化石墨烯在马弗炉里面还原，结果被烧成灰了
已经有32人回复
- 普通石墨粉加到铝熔体中会氧化变成二氧化碳？
已经有8人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
tfl03: 金币+10, ★有帮助, 我想做的是石墨，不是石墨烯，能给点建议吗？ 2012-06-24 15:13:36
chaizhm: 金币+4, 谢谢~ 2012-06-25 08:57:44
感谢参与，应助指数 +1
tfl03: 金币+10, ★有帮助, 我想做的是石墨，不是石墨烯，能给点建议吗？ 2012-06-24 15:13:36
chaizhm: 金币+4, 谢谢~ 2012-06-25 08:57:44
|
好久不用lammps了，懒得细看in文件， 我给你推荐个更简单的方法建立石墨烯模型 VMD，里面有插件直接可以生成你所需要的模型大小和性质。 比如我建立的一个正方形，中间有个孔的石墨烯单层，要转换成lammps数据，只需要载入模型： mol load psf g_hole.psf pdb g_hole.pdb 然后输入命令： topo writelammpsdata ok.data 就得到了lammps的data文件 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : g_hole.psf
- 附件 2 : g_hole.pdb
- 附件 3 : ok.data
2012-06-24 15:12:18, 580.92 K
2012-06-24 15:12:19, 264.61 K
2012-06-24 15:12:29, 489.86 K
2楼2012-06-24 15:11:56
3楼2012-06-24 15:19:02
tfl03
木虫 (职业作家)
- 应助: 24 (小学生)
- 金币: 3128.5
- 散金: 1126
- 红花: 42
- 帖子: 4453
- 在线: 1812.2小时
- 虫号: 772313
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-06-24 15:20:40
5楼2012-06-24 15:22:36
tfl03
木虫 (职业作家)
- 应助: 24 (小学生)
- 金币: 3128.5
- 散金: 1126
- 红花: 42
- 帖子: 4453
- 在线: 1812.2小时
- 虫号: 772313
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2012-06-24 15:24:16
tfl03
木虫 (职业作家)
- 应助: 24 (小学生)
- 金币: 3128.5
- 散金: 1126
- 红花: 42
- 帖子: 4453
- 在线: 1812.2小时
- 虫号: 772313
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
今天又实验了好几次，还是不对。lammps的graphite输入文件错在哪里吗？ units metal boundary p p p dimension 3 atom_style atomic #lattice custom 2.46 a1 1 0 0 a2 -0.5 0.866 0 a3 0 0 2.74 & #basis 0 0 0 basis 0.333 0.667 0 & #basis 0 0 0.5 basis 0.667 0.333 0.5 #region box block 0 1 0 1 0 1 #create_box 1 box #create_atoms 1 box #lattice custom 1.42 a1 1.5 -0.866 0 a2 1.5 0.866 0 a3 0 0 5.46 & # basis 0 0 0 basis 0.5 0.866 0 basis 0 0 2.73 #lattice custom 1.42 a1 1.732 0 0 a2 0.866 1.5 0 a3 0.866 0.5 2.73 & # basis 0 0 0 basis 0 1 0 #lattice custom 1.42 a1 1.5 -0.866 0 a2 1.5 0.866 0 a3 0 0 4.7887 & # basis 0 0 0 basis 0.5 0.866 0 basis 0 0 0.5 lattice custom 2.46 a1 0.5 -0.866 0 a2 0.5 0.866 0 a3 0 0 2.73 & basis 0 0 0 basis 0.5 0.433 0 & basis 0 0 0.5 basis 0.5 -0.433 0.5 region myr1 prism 0 1 0 1 0 1 0 0 0 create_box 1 myr1 create_atoms 1 region myr1 |
7楼2012-06-25 07:40:05
| 好 |
8楼2012-06-25 08:09:48
 |
9楼2012-06-25 08:41:18
10