sobereva 发表于 2021-6-18 01:21 老师我还有一个问题,就是我已经用这个质量加权的hessian矩阵将笛卡尔坐标下两个电子态的位移,也就是△Q的平方,我看这个计算黄昆因子的公式是s=(w/2h)*△Q2,如图所示,如果△Q是无量纲的话,那黄昆因子就有单位了是J/s2,请问一下这个单位我该如何转化呢。 |
sobereva 发表于 2021-6-18 01:21 谢谢老师 |
|
利用这里面的代码可以基于fch里的Hessian计算本征矢 基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq http://sobereva.com/328(http://bbs.keinsci.com/thread-3368-1-1.html) 没有单位 |
wzkchem5 发表于 2021-6-17 23:36 哇,谢谢啦,那我还想问一个问题就是这个本征矢有单位么 |
任嘉人amber 发表于 2021-6-17 16:34 抱歉我没看到你问的是这行字“下面”的数据是什么意思。 是mass-weighted Hessian的本征矢。不过建议从fchk里读取本征矢数据,那里读到的有效数字位数比较多 |
wzkchem5 发表于 2021-6-17 23:27 可能是我没表述清楚,我想知道这个深层含义,是不是代表Hessian矩阵的本征向量矩阵呢,因为我想将笛卡尔坐标转化为简正坐标改如何转化呢 |
|
coord 1=x, 2=y, 3=z atom就是原子编号 element是原子序数,比如6就代表碳 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-17 23:20 , Processed in 0.169746 second(s), 25 queries , Gzip On.