| 5 | 1/1 | 返回列表 |
| 查看: 3280 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
風(fēng)吹葉123金蟲 (小有名氣)
|
[求助]
【求助】關(guān)于LAMMPS建模
|
|
|
各位版主、大神: 你們好。我想問幾個(gè)關(guān)于LAMMPS建模的基本問題。我知道可能有點(diǎn)簡(jiǎn)單,但我真的是認(rèn)真查看過Manual中相關(guān)的命令說明之后,仍有一些疑問。我自己想了很久,也問過很多人,但是還是沒有解決。 命令來自一個(gè)很基礎(chǔ)的LAMMPS例子,該例子的全部命令如下: # Find minimum energy fcc configuration # Mark Tschopp, 2010 # ---------- Initialize Simulation --------------------- clear units metal dimension 3 boundary p p p atom_style atomic atom_modify map array # ---------- Create Atoms --------------------- lattice fcc 4 region box block 0 1 0 1 0 1 units lattice create_box 1 box lattice fcc 4 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 create_atoms 1 box replicate 1 1 1 # ---------- Define Interatomic Potential --------------------- pair_style eam/alloy pair_coeff * * Al99.eam.alloy Al neighbor 2.0 bin neigh_modify delay 10 check yes # ---------- Define Settings --------------------- compute eng all pe/atom compute eatoms all reduce sum c_eng # ---------- Run Minimization --------------------- reset_timestep 0 fix 1 all box/relax iso 0.0 vmax 0.001 thermo 10 thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms min_style cg minimize 1e-25 1e-25 5000 10000 variable natoms equal "count(all)" variable teng equal "c_eatoms" variable length equal "lx" variable ecoh equal "v_teng/v_natoms" print "Total energy (eV) = ${teng};" print "Number of atoms = ${natoms};" print "Lattice constant (Angstoms) = ${length};" print "Cohesive energy (eV) = ${ecoh};" print "All done!" 我的問題是: (1)lattice fcc 4 和 lattice fcc 4 orient x 100 orient y 0 1 0 orient z 0 0 1這兩句有什么區(qū)別?效果不是一樣的么?為什么在lattice fcc 4建立格子之后,還要有l(wèi)attice fcc 4 orient x 100 orient y 0 1 0 orient z 0 0 1這句?意義何在? (2)replicate 1 1 1這句不相當(dāng)于沒有復(fù)制么?為什么要寫這條命令呢?意義何在? (3)如果我只想建模,不進(jìn)行任何計(jì)算,然后想在LAMMPS中導(dǎo)出該模型,最后在VMD中查看該模型。那么在整個(gè)代碼的最后一行需要添加一個(gè)怎樣的dump語(yǔ)句?(我查過手冊(cè)中的dump語(yǔ)句,也試過很多次自己寫,但是都不對(duì)。) 跪求各路大神幫助。小弟在此謝過! |
金蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
|
對(duì)于問題1:好好看看手冊(cè),說的蠻詳細(xì)。The other 3 arguments are the crystallographic direction in the lattice that you want to orient along that axis, specified as integers. E.g. "orient x 2 1 0" means the x-axis in the simulation box will be the [210] lattice direction. 兩條命令有差別; 對(duì)于問題2:我覺得可以去掉; 對(duì)于問題3:dump 1 all atom 1 dump.lammpstrj 就可以了,其中的1不可以和其他的dump重復(fù);1000為每多少步輸出一次,可以自己設(shè)定,看手冊(cè)就知道啥意思了。關(guān)鍵在于文件名后綴為lammpstrj。 PS:我沒仔細(xì)看你的程序,只是針對(duì)你的問題。你可以試下! |

金蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 080500-315分復(fù)試調(diào)劑 +4 | 上岸3821 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
130+3
|
Yishengeryi 2026-03-30 | 3/150 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|