| 9 | 1/1 | 返回列表 |
| 查看: 1315 | 回復(fù): 8 | ||
tfl03木蟲 (職業(yè)作家)
|
[求助]
Lammps: 請(qǐng)幫忙建立一個(gè) 石墨 graphite 的輸入晶格模型,發(fā)短息亦可
|
|
用下面的語(yǔ)句建立模型,老出問題,主要是開始原子個(gè)數(shù)都不對(duì),謝謝! units metal boundary p p p dimension 3 atom_style atomic #想建立石墨graphite的晶格,根據(jù)MS的數(shù)據(jù),六方結(jié)構(gòu),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個(gè)原子,丟了4個(gè) region myr1 prism 0 1 0 1 0 1 0 0 0 create_box 1 myr1 create_atoms 1 region myr1 #同樣說有6個(gè)原子,丟了4個(gè) 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 ] |

|
好久不用lammps了,懶得細(xì)看in文件, 我給你推薦個(gè)更簡(jiǎn)單的方法建立石墨烯模型 VMD,里面有插件直接可以生成你所需要的模型大小和性質(zhì)。 比如我建立的一個(gè)正方形,中間有個(gè)孔的石墨烯單層,要轉(zhuǎn)換成lammps數(shù)據(jù),只需要載入模型: mol load psf g_hole.psf pdb g_hole.pdb 然后輸入命令: topo writelammpsdata ok.data 就得到了lammps的data文件 |
木蟲 (職業(yè)作家)

木蟲 (職業(yè)作家)

木蟲 (職業(yè)作家)
|
今天又實(shí)驗(yàn)了好幾次,還是不對(duì)。lammps的graphite輸入文件錯(cuò)在哪里嗎? 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 |

| 好 |
|
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3 | 小張zxy張 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +4 | 辰沐5211314 2026-03-26 | 4/200 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +3 | Micky11223 2026-03-25 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
[考研] 調(diào)劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 奶油草莓. 2026-03-22 | 7/350 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|