| 5 | 1/1 | 返回列表 |
| 查看: 2582 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
lammps的minimize會(huì)丟失原子 已有3人參與
|
|||
| 今天使用了8月8日的版本運(yùn)行下面的腳本,出現(xiàn)lost atom 錯(cuò)誤。但使用以前 2013 12月的版本卻可以運(yùn)算。然后用count(all,whole)數(shù)了一下box中的原子數(shù),發(fā)現(xiàn)以前的版本中的原子確實(shí)丟了,只是由于沒檢查,所以可以運(yùn)行。最新版的要10個(gè)步長檢查一次守恒條件,所以出錯(cuò)。希望大家在計(jì)算的時(shí)候注意一下。另外有誰知道lost atoms的原因么? |


|
本帖內(nèi)容被屏蔽 |
|
# Input file for bi-layered crystal, Nickle and Aluminum, relaxed # Tikelia Wright, July 2012 # ------------------------ INITIALIZATION ---------------------- clear units metal dimension 3 boundary p p p atom_style atomic # ----------------------- ATOM DEFINITION ----------------------- lattice fcc 3.79 region upper block 0 9 0 9 0 2.999 lattice fcc 3.79 region lower block 0 9 0 9 3 7 region whole union 2 upper lower create_box 2 whole lattice fcc 3.79 create_atoms 1 region upper lattice fcc 3.79 create_atoms 2 region lower # ------------------------ FORCE FIELDS ----------------------- pair_style eam/alloy pair_coeff * * Mishin-Ni-Al-2009.eam.alloy Ni Al variable x equal count(all,whole) #---------------------------Settings---------------------------- #----------------------Run Minimization------------------------- variable length equal lx reset_timestep 0 fix 1 all box/relax iso 0.0 vmax 0.001 thermo_style custom v_x lx thermo 1 dump 1 all custom 10 test*.txt xs ys zs dump 2 all xyz 10 test*.xyz min_style cg minimize 1e-15 1e-15 5000 5000 # SIMULATION DONE print "All done" |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +5 | 腳滑的守法公民 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +3 | hanamiko 2026-03-27 | 3/150 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +10 | 幸運(yùn)的醬醬 2026-03-22 | 12/600 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|