Amber进阶-黑塞矩阵在哪里

  1. 找手册

    • 1-7章提过了
    • 8章:经验共价键键模拟 计算反应
      基函数的概述 单个基函数由一个或多个原始高斯函数构成。例如,s-型基函数 φµ(r)是一个高斯分布型函数称为s-型基函数,一个壳层是一套有共同指数的基函数φµ。Gaussian支持任意角动量的壳层:s,p,d,f,g,h五层。一个s-壳层包含一个单独的s-型基函数。一个p-壳层包含三个基函数pX,pY和pZ。sp-壳层包含四个具有共同高斯指数的基函数:一个s-型函数和三个p-函数pX,pY和pZ。
    • AMBER18里搜hessian发现有一个关于second derivative的,但似乎用到mme2,不懂,也没搜到
      有一个CartHess2FC.py的模块用了output Hessian 以为是它,实际上利用的是Gaussian输出的fchk文件
      小木虫上看到很成熟的利用Gaussian解决hessian乃至找本征频率的做法,http://muchong.com/html/201504/8775391_2.html,15年就有了
    • 改变搜索思路,找fchk,找到了量子动力学模块
      PIMD可以做模拟,一种是PRIMPIMD和NMPIMD(normal mode path-integral),后者是一种更好的方式,解耦谐波项(在谐波项看到了正则化中的频率ω)

      流泪了好兄弟

      实践证明,以hessian为关键词搜一遍,发现,但凡涉及到的都是偷Gaussian的现成文件来用,也太垃圾了amber,呜呜呜,浪费我一个月啊!!!!!!!!!

      hessian的本质是3n-6个振动自由度的频率值,3n个自由度,但去掉3个振动自由度,3个转动自由度,他们的改变不影响体系相对位置,也就不改变体系能量。