计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 哈密顿矩阵lowdin transformation求助

哈密顿矩阵lowdin transformation求助

查看数: 8054 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

jszhang

发布时间: 2019-10-8 14:55

正文摘要:

本帖最后由 jszhang 于 2019-10-10 10:06 编辑请问有没有人做过关于Löwdin transformation的计算呢？这里S0是重叠矩阵，H0是哈密顿矩阵，我想要得到经过Löwdin变换后的H矩阵，应该怎么计算呢？有大神 ...

回复 Reply

jszhang 发表于 Post on 2019-10-8 17:05:17

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

好的，谢谢老师

jszhang 发表于 Post on 2019-10-8 17:01:37

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

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

sobereva 发表于 Post on 2019-10-8 16:56:24

就是基本的矩阵运算而已
细节在Szabo的Modern quantum chemistry里写得很明白
只要会编程，知道怎么读取H、S矩阵，实现这个轻而易举

让你变成回忆 发表于 Post on 2019-10-8 16:47:28

本帖最后由让你变成回忆于 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矩阵做一个对角化。

上述仅是个人理解。

jszhang 发表于 Post on 2019-10-8 16:43:50

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

不好意思，已更改

sobereva 发表于 Post on 2019-10-8 16:08:54

发帖之后仔细检查帖子内容，不要令别人阅读困难
量化问题不要发到分子模拟版，给你移动了

zjxitcc 发表于 Post on 2019-10-8 15:18:21

建议先把格式编辑好··从标题到内容都有乱码

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register