24小时热门版块排行榜    

查看: 3102  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 面上函评意见出来了,像什么等级? 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 投票:  有多少人是今天查系统知道结果的? +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
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +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
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +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 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +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
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见