应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
91/1返回列表
查看: 4089  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

江流儿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步
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
倚天照海花无数，高山流水心自知。
1楼 2013-05-21 20:45:16
已阅   回复此楼   关注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 为非周期，所以报错误
赞一下(1人) 回复此楼
专业文献求助
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
关闭
请填处理意见
关闭