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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】siesta 做NOSE动力学出错，请高手指教
121/2返回列表12下一页
查看: 1765  |  回复: 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的回帖

gavinliu7390

木虫 (著名写手)

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

wangsong1016

铜虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
youzhizhe(金币+2, 1ST强帖+1): 谢谢交流。 2011-05-03 10:48:12
引用回帖:
Originally posted by 轻风龙骑士 at 2011-05-02 14:59:18:
谢谢啊，后来我查了资料，明白了。现在我问题是：是不是一般认为如果跑一个恒温动力学，在温度变化范围达到设定温度的上下10%的时候，在设定温度下达到了平衡态？那么如果初始构型不是设定温度下的稳定相，是不 ...

恒温的动力学过程我没有计算过，所以对过程不太了解。至于10%是否合适我觉得也很难界定，这根据你的体系的大小不同应该是不一样的。体系越小，涨落就越大。计算50个原子时涨落应该大于500个原子，两个体系都用10%应该是不合适的。我这里也不敢给你定量的结论，很不好意思。
    还有一个是我认为你所说的相变应该是宏观的一个过程，某个温度下的稳定相应该是在固体的宏观尺度范围内进行考虑，但是看固体的微观结构，还是应该有很剧烈的涨落的。计算到百十个原子的时候，无论涨落是否很大应该都很难反映出宏观固体状态。
    我有个想法就是如果你想做一个恒温的计算，可否考虑计算时间长一点。比如你要一个100K的温度，你让体系自己运行，跑个2000或者3000个MD过程，然后再看一下他的温度变化范围，有可能就是在你100K的温度上下有较大的波动。但是你把3000个温度取个平均值试试，也许就是100K左右的温度。温度本来就是一个较长时间对能量统计平均值，如果做了这么长那个时间的模拟，温度平均值能够是你所希望得到的那个值，我觉得应该是比较有说服力了。
赞一下(2人) 回复此楼
6楼2011-05-03 10:34:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

轻风龙骑士

新虫 (初入文坛)
引用回帖:
Originally posted by wangsong1016 at 2011-05-03 10:34:49:
恒温的动力学过程我没有计算过，所以对过程不太了解。至于10%是否合适我觉得也很难界定，这根据你的体系的大小不同应该是不一样的。体系越小，涨落就越大。计算50个原子时涨落应该大于500个原子，两个体系都用 ...

谢谢指教，我后来选了实验相变温度附近的几个温度点进行恒温动力学，每个点跑1万步，现在已经出来结果了，数据分析证明得到了我需要的相变现象，呵呵！
本来我还想着我的温度涨落怎么还是不到10%以内，原来是我的体系选的没有那么大，现在我可以放心我的结果了，再次感谢点拨！
赞一下 回复此楼
自然如何呈现给我们，取决于我们用多少时间去观察他
7楼2011-05-04 11:17:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
谈不上指教，都是有了问题大家出来一起想办法，以后多交流了。
赞一下(1人) 回复此楼
8楼2011-05-05 09:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

轻风龙骑士

新虫 (初入文坛)
引用回帖:
Originally posted by wangsong1016 at 2011-05-05 09:29:59:
谈不上指教，都是有了问题大家出来一起想办法，以后多交流了。

对了，还想问一个问题，就是我之前听人说过，说SIESTA中有个什么“线性....”什么什么的一个TAG，默认值是不开的，开下来计算速度就会提高很多，你知道吗？我在手册上找到很多“LINE”有关的TAG，可是不知道是那个…
赞一下 回复此楼
自然如何呈现给我们，取决于我们用多少时间去观察他
9楼2011-05-05 23:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
好像在论坛里有人提过，应该是比较长时间了，我也不记得了。那个应该还是要付出代价的，可能精确度差一点吧。我不太清楚，你自己在查一下手册吧，到底是什么我也不太清楚。
赞一下(1人) 回复此楼
10楼2011-05-06 13:25:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 轻风龙骑士 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +14 于海海海海 2026-03-19 14/700 2026-03-19 19:40 by maocaozhuxi
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 085600材料与化工 +5 安全上岸！ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭