| 5 | 1/1 | 返回列表 |
| 查看: 1319 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
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 ] |

木蟲 (職業(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 |

|
好久不用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è)作家)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +3 | 阿姨復(fù)古不過 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +5 | 一定上岸aaa 2026-03-27 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +3 | hanamiko 2026-03-27 | 3/150 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績(jī)343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|