| 5 | 1/1 | 返回列表 |
| 查看: 3270 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
風(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) (小有名氣)
金蟲(chóng) (小有名氣)
鐵桿木蟲(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) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚(yú)ly 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫(xiě)郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|