24小时热门版块排行榜    

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

digghost

木虫 (著名写手)

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

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

» 猜你喜欢

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

呵呵了
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
信息提示
请填处理意见