| 10 | 1/1 | 返回列表 |
| 查看: 3268 | 回復(fù): 9 | ||
風(fēng)吹葉123金蟲 (小有名氣)
|
[求助]
【求助】關(guān)于LAMMPS建模
|
|
各位版主、大神: 你們好。我想問幾個關(guān)于LAMMPS建模的基本問題。我知道可能有點(diǎn)簡單,但我真的是認(rèn)真查看過Manual中相關(guān)的命令說明之后,仍有一些疑問。我自己想了很久,也問過很多人,但是還是沒有解決。 命令來自一個很基礎(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)行任何計算,然后想在LAMMPS中導(dǎo)出該模型,最后在VMD中查看該模型。那么在整個代碼的最后一行需要添加一個怎樣的dump語句?(我查過手冊中的dump語句,也試過很多次自己寫,但是都不對。) 跪求各路大神幫助。小弟在此謝過! |
金蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
|
對于問題1:好好看看手冊,說的蠻詳細(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. 兩條命令有差別; 對于問題2:我覺得可以去掉; 對于問題3:dump 1 all atom 1 dump.lammpstrj 就可以了,其中的1不可以和其他的dump重復(fù);1000為每多少步輸出一次,可以自己設(shè)定,看手冊就知道啥意思了。關(guān)鍵在于文件名后綴為lammpstrj。 PS:我沒仔細(xì)看你的程序,只是針對你的問題。你可以試下! |

金蟲 (小有名氣)
金蟲 (著名寫手)

禁蟲 (初入文壇)
|
本帖內(nèi)容被屏蔽 |
鐵蟲 (小有名氣)
新蟲 (正式寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 298求調(diào)劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|