请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 129|回复: 7

[算法与编程] 哈密顿矩阵lowdin transformation求助

[复制链接]

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

发表于 6 天前 | 显示全部楼层 |阅读模式
本帖最后由 jszhang 于 2019-10-10 10:06 编辑

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

2WYTG7}28YS]C$Z3H}4IN_7.png

473

帖子

0

威望

2434

eV
积分
2907

Level 5 (御坂)

发表于 6 天前 | 显示全部楼层
建议先把格式编辑好··从标题到内容都有乱码

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 6 天前 | 显示全部楼层
发帖之后仔细检查帖子内容,不要令别人阅读困难
量化问题不要发到分子模拟版,给你移动了
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
zjxitcc 发表于 2019-10-8 15:18
建议先把格式编辑好··从标题到内容都有乱码

不好意思,已更改

360

帖子

2

威望

1916

eV
积分
2316

Level 5 (御坂)

发表于 6 天前 | 显示全部楼层
本帖最后由 让你变成回忆 于 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矩阵做一个对角化。

上述仅是个人理解。

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 6 天前 | 显示全部楼层
就是基本的矩阵运算而已
细节在Szabo的Modern quantum chemistry里写得很明白
只要会编程,知道怎么读取H、S矩阵,实现这个轻而易举
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

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

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

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
sobereva 发表于 2019-10-8 16:56
就是基本的矩阵运算而已
细节在Szabo的Modern quantum chemistry里写得很明白
只要会编程,知道怎么读取H ...

好的,谢谢老师
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-10-14 22:27 , Processed in 0.167342 second(s), 27 queries .

快速回复 返回顶部 返回列表