| 5 | 1/1 | 返回列表 |
| 查看: 3275 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
風(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ì)看你的程序,只是針對你的問題。你可以試下! |

金蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +6 | WW.' 2026-03-23 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 299求調(diào)劑 +8 | 15188958825 2026-03-25 | 8/400 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +7 | 芋泥小鈴鐺 2026-03-28 | 7/350 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +10 | 幸運(yùn)的醬醬 2026-03-22 | 12/600 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|