24小时热门版块排行榜    

CyRhmU.jpeg
查看: 491  |  回复: 0

tristan9362

银虫 (初入文坛)

[交流] 【求助】Nose hoover chain 的code问题(希望写过NHC的可以路过指教一下)

本人在按照molecular physics 1996 87 1117写NHC的时候看到其中用了call shake,这个是必须要用吗,主要是Daan Frenkel 的Molecular Simulation书里case 13的附注没有提这件事
现在我的结果仍是一个环,能量也不守恒

我写的思路是这样的(和文献一样):
定义初始的位移,动量,热浴的位移,动量
体系是一个simple harmonic oscillator

第一个dt/2 循环,
然后是velocity verlet 跟新位移和动量
第二个dt/2循环

如果需要的话,我可以贴出code,或者大家讲一下注意点也可以。。
回复此楼

» 猜你喜欢

化学好像不是我的唯一出路....
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 tristan9362 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见