| 查看: 2230 | 回复: 6 | |||
维特的日记木虫 (小有名气)
|
[求助]
二氧化硅包裹烷烃的核壳结构模拟温度和能量不守恒问题 已有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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有13人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
fuhuidong
铁杆木虫 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 8755.3
- 红花: 12
- 帖子: 6933
- 在线: 483.1小时
- 虫号: 2433013
- 注册: 2013-04-23
- 性别: GG
- 专业: 胶体与界面化学
2楼2015-12-26 21:55:34
维特的日记
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 3814.9
- 红花: 2
- 帖子: 169
- 在线: 55.5小时
- 虫号: 3159711
- 注册: 2014-04-24
- 专业: 工程热物理相关交叉领域
3楼2015-12-26 22:11:45
fuhuidong
铁杆木虫 (知名作家)
- 应助: 1 (幼儿园)
- 金币: 8755.3
- 红花: 12
- 帖子: 6933
- 在线: 483.1小时
- 虫号: 2433013
- 注册: 2013-04-23
- 性别: GG
- 专业: 胶体与界面化学
4楼2015-12-27 09:51:54
Liupeng-MD
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 358.5
- 帖子: 81
- 在线: 41.6小时
- 虫号: 3926764
- 注册: 2015-06-16
- 性别: GG
- 专业: 机械摩擦学与表面技术
5楼2018-01-18 23:54:44
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
6楼2018-01-19 00:23:47
7楼2019-03-25 21:06:11
50