24小时热门版块排行榜    

查看: 2554  |  回复: 6

digghost

木虫 (著名写手)

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

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

» 猜你喜欢

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

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

huachao90

禁虫 (正式写手)


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

2楼2014-08-12 00:32:48
已阅   回复此楼   关注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的回帖

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的回帖

gulubaozi

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by digghost at 2014-08-12 13:27:12
# Input file for bi-layered crystal, Nickle and Aluminum, relaxed
# Tikelia Wright, July 2012

# ------------------------ INITIALIZATION ----------------------
clear
units                 metal
dimension        3
...

可能是minimize 1e-15 1e-15 5000 5000里边的min极限太小所致
6楼2014-08-13 08:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanna_zhu

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问楼主问题解决了吗?我也遇到同样的问题了,跪求帮助啊
7楼2018-10-30 19:21:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 digghost 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +4 小贾同学123 2026-03-15 6/300 2026-03-17 07:06 by laoshidan
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 5/250 2026-03-17 06:56 by Muziwl
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
请填处理意见