24小时热门版块排行榜    

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

digghost

木虫 (著名写手)

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

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

» 猜你喜欢

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

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

gulubaozi

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
附上相关命令看看
4楼2014-08-12 11:52:24
已阅   回复此楼   关注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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕297已过四六级求调剂推荐 +4 adaie 2026-03-19 4/200 2026-03-20 09:31 by yuncha
[考研] 085600材料与化工 求调剂 +17 enenenhui 2026-03-13 18/900 2026-03-20 08:55 by xingguangj
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 0703化学调剂 ,六级已过,有科研经历 +12 曦熙兮 2026-03-15 12/600 2026-03-19 19:42 by maocaozhuxi
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
信息提示
请填处理意见