24小时热门版块排行榜    

查看: 2474  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见