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

金蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)
金蟲(chóng) (著名寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專(zhuān)碩,293,求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 0703化學(xué) +10 | 妮妮ninicgb 2026-03-27 | 10/500 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +8 | 辰沐5211314 2026-03-26 | 8/400 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專(zhuān)業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|