24小时热门版块排行榜    

查看: 3099  |  回复: 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.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 06:06 by o3EMxHUAQDtS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:56 by o3EMxHUAQDtS
[基金申请] 基金未中,这种答复是模板吗? +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
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +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
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见