24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3673  |  回复: 8

江流儿0305

新虫 (初入文坛)

[求助] lammps实验拉伸Al纳米线的in文件的错误提示好久找不出? 已有2人参与

以下是in文件一直提示错误ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)
还有WARNING: Resetting reneighboring criteria during minimization (min.cpp:173)求高手。
# 3d EAM tensile simulation---------三维EAM拉伸仿真

#units         metal        #设定原子模型
dimension     3        #几何维度

boundary     s s p        #边界设定
atom_style     atomic

# create geometry---------------创建几何模型

lattice        fcc 0.4056 orient x 1 1 -2 orient y -1 1 0 orient z 1 1 1 origin 0.0 0.0 0.5  #晶格结构,晶格常数x,y,z坐标的晶向
region        box block 0 4 0 4 0 20    #模拟几何区域设定
create_box 1 box        #创建一个box
create_atoms 1 box         #在box中创建1一种原子
mass        1 26.981    #定义原子质量为26.981

# EAM potentials----------------势函数

pair_style    eam        #定义势函数
pair_coeff    * * Al_jnp.eam    #势所需要的参数在此文件里

neighbor    2.0 bin        #截断半径
neigh_modify    delay 1 check yes #邻近原子列表更新速度

# define groups----------------定义组

region         1 block INF INF INF INF INF 3         #定义一个叫1的区域
group        lower region 1                #此区域里的原子叫lower
region        2 block  INF INF INF INF 17 INF        #定义一个叫2的区域
group        upper region 2                #此区域里的原子叫upper
group        boundary union lower upper            #定义lower+upper = boundary
group        mobile subtract all boundary        #定义mobile=all- boundary

#initial velocities-------------设置原子初速度

velocity    lower set 0.0 0.0 0.0    #设置left区域原子初速度为0
velocity    upper set 0.0 0.0 0.3   #设置right区域z方向原子速度为0.3
velocity    mobile ramp vz 0.0 0.3 z 3 17 sum yes   

min_style sd                     #能量最小模式   
minimize 1.0e-6 1.0e-6 1000 10000        #驰豫

fix   1 lower setforce 0.0 NULL NULL     #固定下边的原子
fix   2 all npt temp 300.0 300.0 100.0 aniso 0.0 0.0 x 10.0 drag 1.0

thermo 200                 #每200步输出一次结果
timestep    0.003                       #时间步长0.003ps
run  5000                 #运行5000步

unfix2
fix   3 all nvt  
fix   4 all temp/rescale 300 1.0 1.0 5 1.0
fix   5 all deform 1000 z erate 0.002 units box  #设置沿z方向拉伸,应变率为0.002(1/ps)

dump 1 all atom 500 dump.tensile       #输出结果到dump.tensile
min_style sd                     #能量最小模式   
minimize 1.0e-6 1.0e-6 1000 10000        #驰豫
run  100000                            #运行10000步
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

倚天照海花无数,高山流水心自知。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 谢谢版主 2013-05-21 22:58:32
Illegal fix nvt/npt/nph command


错误提示很明显
2楼2013-05-21 21:06:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江流儿0305

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2013-05-21 21:06:40
Illegal fix nvt/npt/nph command


错误提示很明显

对啊,但就是找不出错误,求教!
倚天照海花无数,高山流水心自知。
3楼2013-05-22 13:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩欢欢

金虫 (初入文坛)

引用回帖:
3楼: Originally posted by 江流儿0305 at 2013-05-22 13:47:03
对啊,但就是找不出错误,求教!...

请问你的问题解决了吗?我也遇到了这个问题,不知道怎么解决。我在拉伸的时候用的是fix  1 nvt 300.0 300.0 100.0 drag 0.2 chain no 是想在300k恒温拉伸。运行也是这个错误提示,请教下怎么修改
Demi
4楼2013-07-17 09:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

【答案】应助回帖

fix   2 all npt temp 300.0 300.0 100.0 aniso 0.0 0.0 x 10.0 drag 1.0
x 为非周期,所以报错误
专业文献求助
5楼2013-07-19 17:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天泊乐

金虫 (正式写手)

【答案】应助回帖

fix npt 命令不能用于非周期性边界,如要使用,涉及压强参数要全为NULL,也就是控不了压了
音乐you我!
6楼2013-08-07 11:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

largehoo

木虫 (小有名气)

【答案】应助回帖

lattice 的单位是埃
7楼2014-07-12 15:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinaboy798

银虫 (小有名气)

【答案】应助回帖

这个in文件有几点问题,个人建议
1.非周期条件不能使用npt
2.为什么你要固定了Z轴两端的速度做能量最小化,而后在拉伸的时候是对所有原子施加npt系综做弛豫?
3.fix nvt已经有一个温度控制了,又加一个temp/rescale?
4.既然最后使用的fix deform来施行拉伸,之前施加Z轴两端速度为零的能量最小化有何意义?
8楼2014-11-11 14:45:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohu199316

新虫 (初入文坛)

引用回帖:
1楼: Originally posted by 江流儿0305 at 2013-05-21 20:45:16
以下是in文件一直提示错误ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)
还有WARNING: Resetting reneighboring criteria during minimization (min.cpp:173)求高手。
# 3d EAM tensile simulation-- ...

请问,Al拉伸in文件还有吗?

发自小木虫Android客户端
9楼2017-04-14 13:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 江流儿0305 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 297,工科调剂?河南农业大学本科 +15 河南农业大学-能 2026-04-14 15/750 2026-04-20 12:39 by Equinoxhua
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见