计算化学公社

标题: 哈密顿矩阵lowdin transformation求助 [打印本页]

作者
Author: jszhang 时间: 2019-10-8 14:55
标题: 哈密顿矩阵lowdin transformation求助
本帖最后由 jszhang 于 2019-10-10 10:06 编辑

请问有没有人做过关于Löwdin transformation的计算呢？这里S0是重叠矩阵，H0是哈密顿矩阵，我想要得到经过Löwdin变换后的H矩阵，应该怎么计算呢？有大神做过类似的计算吗？

作者
Author: zjxitcc 时间: 2019-10-8 15:18
建议先把格式编辑好··从标题到内容都有乱码

作者
Author: sobereva 时间: 2019-10-8 16:08
发帖之后仔细检查帖子内容，不要令别人阅读困难
量化问题不要发到分子模拟版，给你移动了

作者
Author: jszhang 时间: 2019-10-8 16:43

zjxitcc 发表于 2019-10-8 15:18
建议先把格式编辑好··从标题到内容都有乱码

不好意思，已更改

作者
Author: 让你变成回忆 时间: 2019-10-8 16:47
本帖最后由让你变成回忆于 2019-10-8 16:51 编辑

感觉这个变换的形式很熟悉，估计在Multiwfn中找到相关的源代码可以实现你这个功能。
比如S0重叠矩阵，通过IOP在Gaussian的log文件中输出，然后Multiwfn中肯定有现成的可以去读这种格式（下三角矩阵）数据的代码。如果H0指的是Fock矩阵，通过F*C=S*C*epsilon，可以算出F=S*C*epsilon*C^-1，这个从Gaussian得到的fchk文件中可以提取出分子轨道系数以及能量；
最终再对H1矩阵做一个对角化。

上述仅是个人理解。

作者
Author: sobereva 时间: 2019-10-8 16:56
就是基本的矩阵运算而已
细节在Szabo的Modern quantum chemistry里写得很明白
只要会编程，知道怎么读取H、S矩阵，实现这个轻而易举

作者
Author: jszhang 时间: 2019-10-8 17:01

让你变成回忆发表于 2019-10-8 16:47
感觉这个变换的形式很熟悉，估计在Multiwfn中找到相关的源代码可以实现你这个功能。
比如S0重叠矩阵，通过 ...

我已经通过Gaussian得到了S0和H0，这里的H0是通过iop(3/33=1)计算得出的核哈密顿矩阵，想得到正交化的H，中间还需要其他的计算吗，还是Multiwfn中可以找到代码直接得出？

作者
Author: jszhang 时间: 2019-10-8 17:05

sobereva 发表于 2019-10-8 16:56
就是基本的矩阵运算而已
细节在Szabo的Modern quantum chemistry里写得很明白
只要会编程，知道怎么读取H ...

好的，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)