应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】siesta 做NOSE动力学出错，请高手指教
51/1返回列表
查看: 1836  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

轻风龙骑士

新虫 (初入文坛)

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

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

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

siesta问题

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
自然如何呈现给我们，取决于我们用多少时间去观察他
1楼 2011-04-12 22:39:41
已阅   回复此楼   关注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的回帖
查看全部 12 个回答

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+2): 谢谢交流。 2011-04-12 23:44:36
正常 。   因为你是用热浴的方法加热。
温度本身就是一个统计平均值。  只是长时间的模拟，求平均后才是真实的温度。
一段时间后， 温度的震荡应该逐渐变小。
赞一下(2人) 回复此楼
真理是一点点接近的！
2楼2011-04-12 23:19:52
已阅   回复此楼   关注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里面有各个原子不同方向的速度。通过公式就可以把各个速度对应的能量转化为温度。
    你在输入文件中给定了初始的温度，程序就给各个原子不同的速度让其震动起来。但是原子一旦发生震动，势能必然变化，那么各个原子的动能（速度）也随之变化，因为总能量要保持守恒。这样就导致了体系温度的变化。这种涨落是不可避免的，不算是错误，但是可能与我们想要的所谓恒定温度情况有较大的偏差。
赞一下(2人) 回复此楼
4楼2011-04-16 22:25:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

轻风龙骑士

新虫 (初入文坛)
引用回帖:
Originally posted by wangsong1016 at 2011-04-16 22:25:16:
我没有用过nose模块，不过我在用verletMD计算时出现过类似的温度起伏，甚至是剧烈起伏的问题。
    我想可能有如下的原因给你参考。
    首先宏观体系的温度就如二楼所说是个统计平均值，实际上是对分子热运动 ...

谢谢啊，后来我查了资料，明白了。现在我问题是：是不是一般认为如果跑一个恒温动力学，在温度变化范围达到设定温度的上下10%的时候，在设定温度下达到了平衡态？那么如果初始构型不是设定温度下的稳定相，是不是平衡后相变已经发生了？
赞一下 回复此楼
自然如何呈现给我们，取决于我们用多少时间去观察他
5楼2011-05-02 14:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334求调剂 +13 Riot2025 2026-04-08 14/700 2026-04-09 00:12 by lature00
[考研] 材料334求调剂 +21 Eecho# 2026-04-03 21/1050 2026-04-08 22:55 by 猪会飞
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +5 守法公民亓纪 2026-04-08 5/250 2026-04-08 21:19 by bljnqdcc
[考研] 324求调剂 +17 想上学求调 2026-04-03 17/850 2026-04-08 20:04 by 我减肥1
[考研] 0703调剂，一志愿天津大学319分 +23 haaaabcd 2026-04-05 26/1300 2026-04-08 16:19 by luoyongfeng
[考研] 一志愿南科大生物学297分，求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 334分机械专硕求调剂 +3 蛋花紫菜汤 2026-04-03 3/150 2026-04-07 14:49 by 逍遥cocoa
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 334分控制工程求调剂 +4 姜尚真sadasd 2026-04-03 4/200 2026-04-07 09:26 by 蓝云思雨
[考研] 华南理工0703化学，总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 085405软件工程301分求调剂，专硕可跨专业，四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 材料专硕322 +11 哈哈哈吼吼吼哈 2026-04-05 11/550 2026-04-06 14:07 by lqwchd
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
信息提示
关闭
请填处理意见
关闭