24小时热门版块排行榜

返回列表

chenjingzhi

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1001.7
散金: 97
红花: 3
帖子: 242
在线: 314.3小时
虫号: 1048640
注册: 2010-06-28
性别: GG
专业: 生物物理、生物化学与分子

[交流] 关于martyna版nose-hoover控温控压方法

在martyna的1996年的论文“Explicit reversible integrators for extended system dynamics”里面提到过，如果要使用这种刘维尔形式的方法进行控温控压的话，那么分子的构形保持算法，如SHAKE、RATTLE，就得进行更改，但是论文里面所提出的受限算法却是由上述的两种混合组成的，第一步（这里使用的是velocity-verlet积分方法）使用SHAKE进行控制，与nose-hoover chain相结合，循环直至得到的受限因子能够使nose-hoover chain操作后的粒子间位置满足键长要求，接着是第二步移动，过程与第一步相似，但是这里使用的却是RATTLE算法，到底这是怎么一回事，是我理解有问题吗。我也找过其他相关的论文，却发现他们给这种构形保持算法取名为ROLL算法，那么说似乎的确与SHAKE和RATTLE不同。求指点。

回复此楼