24小时热门版块排行榜    

查看: 3096  |  回复: 19

牛很牛

新虫 (小有名气)

[交流] 水分子模拟 已有5人参与

我现在在做水分子的模拟,但是现在设置初始温度为150,初始速度为400,但是运行不几步之后,水分子的温度就为零了,然后整个系统中的分子都是在非常缓慢的移动,几乎都不动。不知道问题所在,求大神指点

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王精英

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
也不知道用的什么软件

发自小木虫IOS客户端
DFT计算
2楼2019-03-20 09:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arsennnic

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
几乎不动的话就说明速度为0了呀,当分子速度为0时统计出的温度也就是0了。建议你把代码贴出来看看。

发自小木虫IOS客户端
3楼2019-03-20 10:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 王精英 at 2019-03-20 09:53:57
也不知道用的什么软件

lammps
4楼2019-03-20 10:06:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
3楼: Originally posted by arsennnic at 2019-03-20 10:04:38
几乎不动的话就说明速度为0了呀,当分子速度为0时统计出的温度也就是0了。建议你把代码贴出来看看。

好,我把温度这是这块截图发出来。稍等
5楼2019-03-20 10:06:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

###########initial setting ###########
#fix    fixshake1 suan shake 0.0001 5 0 b 1 2 3 4 a 1 2 3  #这一行有问题,要问
fix 10 suan rigid molecule
velocity all  create 400 49284
timestep     1.0

###############  compute  #############
compute 1 Na  temp
compute 2 Cl  temp
compute 3 suan temp
compute 4 water temp
compute 5 all pe
compute ake all ke

compute energy all group/group suan

thermo_style  custom step temp c_1 c_2 c_3
thermo       1000
thermo_modify temp 1

fix      1  all      nvt temp 300.0 300.0 100.0
fix  rec all recenter INIT INIT INIT

dump     1   all custom   1000 ch0110.water id type x y z

run 1000000

write_data restart.water
6楼2019-03-20 10:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

看是要沉下去了,不要啊
7楼2019-03-20 15:07:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

8楼2019-03-20 22:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nater1ver

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你用的什么单位,timestep设成1.0,而且你的fix nvp的tdamp又是100.0,那频率怕不是1/100年

发自小木虫Android客户端
9楼2019-03-21 00:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
8楼: Originally posted by chuntao118 at 2019-03-20 22:00:03
自己写个程序呗!

这不就是我的程序

发自小木虫Android客户端
10楼2019-03-21 08:46:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像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 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见