24小时热门版块排行榜    

查看: 1431  |  回复: 5

ynwtg

银虫 (初入文坛)

[求助] fix npt 出错

ERROR on proc 0: Pair distance < table inner cutoff <../pair_table.cpp:112>
请教各位大神,是什么原因。不胜感激!
## in file
units        metal
atom_style    atomic
boundary    p p p

lattice        custom ... a1... a2 ... a3... &
                  basis ...
region            box block 0 10 0 10 0 4
create_box        5 box
create_atoms        5 box basis ...

mass                ...
mass                ...
...

# potential

pair_style    table linear 5000 ewald

pair_coeff    1 1 pair.table Y_Y 10.0
pair_coeff          ....                    ....
...                      ....                    ....

run             0
variable    E equal pe
variable     N equal atoms
print           "The number of atoms is $N and the system energy is $E."

dump         1 all atom  ...

neighbor    0.2 bin
neigh_modify    every 1 delay 1 check yes

min_style    cg
minimize    1.0e-10 0 10000 100000

timestep    0.00001                          #步长会不会太小?
thermo        1

fix        npt all npt temp 1000.0 1000.0 1.0 iso 0.0 0.0 100.0 drag 0.2

run        10000
print           "The final state:$N atoms, total energy is $E."

## screen output

fix npt 出错
screen output.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mz121

铜虫 (小有名气)

2楼2015-09-24 00:14:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ynwtg

银虫 (初入文坛)

Tdamp=1, Pdamp=100,都不大吧
3楼2015-09-24 08:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ynwtg

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by mz121 at 2015-09-24 00:14:49
Tdamp Pdamp怎么取那么大

Tdamp=1, Pdamp=100,都不大吧
4楼2015-09-24 08:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fenqin

铜虫 (小有名气)

我也遇到了fix npt 问题是不是
fix        npt all npt temp 1000.0 1000.0 1.0 iso 0.0 0.0 100.0 drag 0.2 这一行中 1000.0 1000.0 必须相等,可以是升温过程吗
5楼2015-09-25 15:57:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ynwtg

银虫 (初入文坛)

引用回帖:
5楼: Originally posted by fenqin at 2015-09-25 15:57:23
我也遇到了fix npt 问题是不是
fix        npt all npt temp 1000.0 1000.0 1.0 iso 0.0 0.0 100.0 drag 0.2 这一行中 1000.0 1000.0 必须相等,可以是升温过程吗

Manual中“The desired temperature at each timestep is a ramped value during the run from Tstart to Tstop”——说明Tstart Tstop可以不等
6楼2015-09-25 16:15:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ynwtg 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:24 by jCd0dEvKHShX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 13:41 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见