计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › Gaussian GTO基组，分子轨道归一化问题。

Gaussian GTO基组，分子轨道归一化问题。

查看数: 12926 | 评论数: 9 | 收藏 Add to favorites 2

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

帖子模式

I10140317

发布时间: 2017-3-20 17:34

正文摘要:

本帖最后由 I10140317 于 2017-3-20 18:10 编辑大家好：最近在写一个程序，涉及到计算分子轨道间的积分<X|Y>，写了一个查看是否归一化的小程序。流程如下， 1. 提取MO coefficient，C，按占据和非占据 ...

回复 Reply

I10140317 发表于 Post on 2017-3-21 12:20:16

非常非常感谢大家的提醒和帮助。

I10140317 发表于 Post on 2017-3-21 12:04:35

sobereva 发表于 2017-3-21 01:24
Gaussian输出的轨道都是归一化的
最好从fch里提取轨道信息，是全精度的
你再细致检查检查，并且建议让系 ...

问题已解决，但不清楚原因，我发现，之前在计算S_AO的时候，加了一个NOSYM，用这套原子轨道积分去算分子轨道积分就会出现问题。去掉之后，就是正交归一的了。

I10140317 发表于 Post on 2017-3-21 10:47:31

本帖最后由 I10140317 于 2017-3-21 10:51 编辑

niobium 发表于 2017-3-21 10:37
你只做HF，然后看看是不是归一化的呢

嗯嗯，我刚刚好像找到点问题所在，可能是原子轨道重叠积分提取的有问题，我再看看，谢谢了。

niobium 发表于 Post on 2017-3-21 10:37:14

你只做HF，然后看看是不是归一化的呢

niobium 发表于 Post on 2017-3-21 10:34:44

哦，是frozen core

niobium 发表于 Post on 2017-3-21 09:50:50

NFC是什么

I10140317 发表于 Post on 2017-3-21 08:35:53

本帖最后由 I10140317 于 2017-3-21 10:20 编辑

sobereva 发表于 2017-3-21 01:24
Gaussian输出的轨道都是归一化的
最好从fch里提取轨道信息，是全精度的
你再细致检查检查，并且建议让系 ...

嗯，谢谢sob老师，我是用rwfdump 524R从rwf文件里提取的数据，这样精度应该和fch里提取的数据精度是一致的。我再仔细检查检查。改成您说的习俗一致的格式试试看，谢谢了。有疑问希望可以再请教。数值是以指数形式读进去的，精度real（kind=8），会影响结果吗？另外，我尝试了您说的用Multiwfn检查，版本为3.4，linux,no_GUI,在wfn.f90里加了call showmatgau(matmul(matmul(transpose(CObasa),Sbas),CObasa),"test",0)，我直接编译后load .fch文件，会报Segment Default，是需要GUI吗？

sobereva 发表于 Post on 2017-3-21 01:24:18

Gaussian输出的轨道都是归一化的
最好从fch里提取轨道信息，是全精度的
你再细致检查检查，并且建议让系数矩阵的列对应各个轨道，行对应各个基函数，这样和习俗一致。
C'SC=I在Multiwfn里可以很容易地验证，在wfn.f90中的!!!--------------------- Now everything start ---------------------!!!的前头增加一句call showmatgau(matmul(matmul(transpose(CObasa),Sbas),CObasa),"test",0)，你会看到载入.fch文件后输出的矩阵是单位矩阵，非对角元的数值可忽略不计。其中CObasa是nbasis*nbasis矩阵（假定轨道数和基函数数目相同，列对应各个MO），Sbas是nbasis*nbasis尺寸的基函数间重叠矩阵。

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