| 查看: 369 | 回复: 0 | ||
relonfbq金虫 (正式写手)
|
[求助]
我这样的程序到底对不对
|
|
我的程序是做分子动力学的，我编程的思想是这样的，用的p-r恒压模型，积分算法是gear算法。 我现在通过一个算例来验证我的程序的正确性（首先我查不出我程序的错误了） 我引入了5*5*5的249个原子的体系（完美晶体应该是250个原子），设置初始温度0K，外压保持0pa，然后运行，运行结果是，前面28000步（步长是1fs）体系的温度在50K以下震荡（或说涨落），因而动能也在很小的数值下涨落，同样，在前面的28000步长里势能，和总能也在约百分之一的范围里涨落。（注意我的数据是200步一平均）。 但是在约30000步左右，体系的温度开始逐步向上升，因而引起动能和总能量逐渐上升（势能和总能的绝对值下降），我运行了40000步，没有看见停止上升的趋势。 当然体系的体积却一直在百分之10的范围内涨落！ 我这样的结果有什么问题呢?而最关键的我的问题（我的程序）出在哪里呢？？？？ [ Last edited by relonfbq on 2011-4-26 at 11:13 ] |
回复此楼» 猜你喜欢
博士去二本高校当辅导员，如何调整心态？
已经有8人回复
-
27细胞生物学申博
已经有5人回复
-
面上有专家说实验设备不是我们单位的
已经有5人回复
-
江苏省自然基金 什么时候出结果
已经有6人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
闲聊
已经有3人回复
-
如何查询2026年获批项目名称？
已经有3人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
