24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

维特的日记

木虫 (小有名气)

应助: 7 (幼儿园)
金币: 3809.9
红花: 2
帖子: 168
在线: 55.5小时
虫号: 3159711
注册: 2014-04-24
专业: 工程热物理相关交叉领域

[求助] 二氧化硅包裹烷烃的核壳结构模拟温度和能量不守恒问题已有1人参与

本人最近一段时间在模拟二氧化硅包裹烷烃的核壳结构，模型是在烷烃的xy方向周围包裹一层二氧化硅，即xy方向是非周期的、z方向是周期的，类似于长方体形状的管道。势函数二氧化硅采用tersoff，烷烃使用Trappe，他们之间的相互作用用lj势。目前在270K下弛豫100万步了，再转到nve弛豫100万步，时间步长为0.1fs。然后出现问题了，模型在270K的nvt弛豫温度还算正常在270上下小幅度振荡，但是体系的总能量不断减小，烷烃的总能量不断增大，二氧化硅的总能量不断减小。nve系综的问题更大，温度随着运行时间缓慢上升，总能量是守恒的，二氧化硅的总能量不断下降，烷烃的总能量不断升高。一开始我以为是时间步长的问题，所以才改成了0.1fs，起初是1fs的，但是改成0.1fs还是不能解决问题。希望虫友们能够帮帮忙，解决这个问题。如下是nvt系综的几个图，nve没有运行万就被我关了，因为已经出现问题了。

1.png

2.png

3.png

4.png

回复此楼

» 猜你喜欢

1楼 2015-12-26 19:25:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

专家经验: +218
模拟EPI: 4
应助: 422 (硕士)
金币: 2294.8
散金: 1720
红花: 88
帖子: 1423
在线: 670.2小时
虫号: 3048965
注册: 2014-03-14
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 分子模拟

【答案】应助回帖

我不觉得有什么问题。nvt过程中总能降低是很自然的。nve过程中总能守恒也是正常的。至于一个势能升高，一个势能降低，那是因为还未达到平衡。你的初始结构可能远离平衡。你能做的是尽量准备接近目标状态的初始构型以及尝试更长时间的模拟。另外，也可以尝试在周期方向用npt算法。

[ 发自手机版 http://muchong.com/3g ]