今天读到一篇论文,其中一段说relaxed force constant方法无需计算昂贵的Hessian,只需算一些单点能量即可。如果这是真的,当然就牛x了。 (, 下载次数 Times of downloads: 10)
上传 Uploaded
点击下载Click to download
然后下载了被引用的文献44,是他们自己的文章,但是对于如何避免Hessian计算只字未提,而是又引用了其他文章,都是上世纪七十年代的别人的文章(这么引用有点不要脸了)。读过之后,发现这些论文并不是讲如何避免Hessian计算,而是因为那个时代还没有发展出解析Hessian的算法,只好用梯度的数值差分来算数值Hessian,这个方法在当时已经比较先进了。比如下面的公式(4)。 (, 下载次数 Times of downloads: 13)
上传 Uploaded
点击下载Click to download
可见通过梯度(并不是能量)算Hessian是理论不成熟导致的,而不是方法的优点。有人写论文好吹牛,但也不能这么吹啊!
ps. 文中提到的relaxed force constant就是Grunenberg的compliance程序计算的relaxed force constant,用了其他程序。文献中报道的这类程序有好几个,只不过Grunenberg的用着比较方便。