24小时热门版块排行榜    

查看: 2563  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

digghost

木虫 (著名写手)

[交流] lammps的minimize会丢失原子 已有3人参与

今天使用了8月8日的版本运行下面的脚本,出现lost atom 错误。但使用以前 2013 12月的版本却可以运算。然后用count(all,whole)数了一下box中的原子数,发现以前的版本中的原子确实丢了,只是由于没检查,所以可以运行。最新版的要10个步长检查一次守恒条件,所以出错。希望大家在计算的时候注意一下。另外有谁知道lost atoms的原因么?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

呵呵了
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digghost

木虫 (著名写手)

引用回帖:
2楼: Originally posted by huachao90 at 2014-08-12 00:32:48
要么是初始构型有问题或者是势能参数有问题,导致作用力过大,位移过大lammps会是别为丢失

但是,我使用的是周期性边界条件呀,而且用的是minimize命令。原子出了box,自然就有其他原子进来。另外刚发现只有32位的系统有这个问题。64位的没有
呵呵了
3楼2014-08-12 06:43:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

huachao90

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

2楼2014-08-12 00:32:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gulubaozi

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
附上相关命令看看
4楼2014-08-12 11:52:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digghost

木虫 (著名写手)

引用回帖:
4楼: Originally posted by gulubaozi at 2014-08-12 11:52:24
附上相关命令看看

# 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"
呵呵了
5楼2014-08-12 13:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 296求调剂 +3 www_q 2026-03-18 6/300 2026-03-19 22:28 by zhq0425
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
信息提示
请填处理意见