24小时热门版块排行榜    

查看: 1766  |  回复: 11
本帖产生 2 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

轻风龙骑士

新虫 (初入文坛)

[交流] 【求助】siesta 做NOSE动力学出错,请高手指教 已有4人参与

在做nose的温控动力学的时候,我设置的初始和目标温度分别是2700K和2600K,做降温的动力学模拟,但是我的systemlable.MDE文件中怎么显示的每一步的温度有问题啊:在第三步的时候就已经到了2620K左右,之后温度更低最低只有1400K左右,然后又开始升高,最高高达3000K以上,这是为什么?是不是我那个参数错了,请高手指教
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

自然如何呈现给我们,取决于我们用多少时间去观察他
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖
谈不上指教,都是有了问题大家出来一起想办法,以后多交流了。
8楼2011-05-05 09:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+2): 谢谢交流。 2011-04-12 23:44:36
正常 。   因为你是用热浴的方法加热。
温度本身就是一个统计平均值。  只是长时间的模拟,求平均后才是真实的温度。
一段时间后, 温度的震荡应该逐渐变小。
真理是一点点接近的!
2楼2011-04-12 23:19:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

轻风龙骑士

新虫 (初入文坛)

引用回帖:
Originally posted by gavinliu7390 at 2011-04-12 23:19:52:
正常 。   因为你是用热浴的方法加热。
温度本身就是一个统计平均值。  只是长时间的模拟,求平均后才是真实的温度。
一段时间后, 温度的震荡应该逐渐变小。

可是,我这个温度区间的中间的点的信息怎么得到呢?如您上面所述的过程得到的结果,又是什么温度的构型呢?是那个温度求平均的温度点?

再请教您一个问题,NoseMass这个参数应该怎么设置呢?跟体系应该没有什么关系吧,是控制什么的?

[ Last edited by 轻风龙骑士 on 2011-4-13 at 13:23 ]
自然如何呈现给我们,取决于我们用多少时间去观察他
3楼2011-04-13 12:52:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+2, 1ST强帖+1): 谢谢交流。 2011-05-02 19:56:53
引用回帖:
Originally posted by gavinliu7390 at 2011-04-12 23:19:52:
正常 。   因为你是用热浴的方法加热。
温度本身就是一个统计平均值。  只是长时间的模拟,求平均后才是真实的温度。
一段时间后, 温度的震荡应该逐渐变小。

我没有用过nose模块,不过我在用verletMD计算时出现过类似的温度起伏,甚至是剧烈起伏的问题。
    我想可能有如下的原因给你参考。
    首先宏观体系的温度就如二楼所说是个统计平均值,实际上是对分子热运动的一个平均描述。
    到了微观层次,按照能量均分定理,每个原子,每个自由度能量是1/2kBT,简化考虑就可以采用1/2kBT=1/2mv**2。意思是在一维无转动的情况下,这一个自由度的能量对应于粒子平动的动能。在程序中就是这么计算得到的温度。输出文件中有一个systemlable.XV里面有各个原子不同方向的速度。通过公式就可以把各个速度对应的能量转化为温度。
    你在输入文件中给定了初始的温度,程序就给各个原子不同的速度让其震动起来。但是原子一旦发生震动,势能必然变化,那么各个原子的动能(速度)也随之变化,因为总能量要保持守恒。这样就导致了体系温度的变化。这种涨落是不可避免的,不算是错误,但是可能与我们想要的所谓恒定温度情况有较大的偏差。
4楼2011-04-16 22:25:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +9 妮妮ninicgb 2026-03-15 13/650 2026-03-19 22:13 by l1006507827
[考研] 085600材料与化工 求调剂 +14 enenenhui 2026-03-13 15/750 2026-03-19 21:53 by allen-yin
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 5/250 2026-03-19 18:03 by zcl123
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 求调剂 +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
信息提示
请填处理意见