计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 波函数分析与Multiwfn (Wavefunction Analysis & Multiwfn) › 关于“跃迁势”

关于“跃迁势”

查看数: 8193 | 评论数: 12 | 收藏 Add to favorites 2

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

帖子模式

肖富贵

发布时间: 2016-1-6 23:27

正文摘要:

最近着手做跃迁电荷（transition charge），Multiwfn可以计算Mulliken transition charge，但此电荷并不完美，因此想自己拟合。拟合分两步，第一步生成静电势；第二步则是用电荷拟合第一步得到的静电势。第一步 ...

回复 Reply

肖富贵 发表于 Post on 2016-1-20 20:53:28

sobereva 发表于 2016-1-20 20:06
可以

多谢多谢多谢Sob前辈，非常不好意思总是麻烦您。

sobereva 发表于 Post on 2016-1-20 20:06:53

肖富贵发表于 2016-1-20 18:56
Sob前辈，非常不好意思，再请问您一个问题：
因为wfn文件没有虚轨道的信息，所以需要用fch文件。请问我 ...

可以

肖富贵 发表于 Post on 2016-1-20 18:56:08

sobereva 发表于 2016-1-20 03:47
合理，但你得把公式里的N乘上去

Sob前辈，非常不好意思，再请问您一个问题：
因为wfn文件没有虚轨道的信息，所以需要用fch文件。请问我可不可以直接用fch文件来计算跃迁势，而不需要把fch文件转换为wfn文件？

肖富贵 发表于 Post on 2016-1-20 04:09:49

sobereva 发表于 2016-1-20 03:47
合理，但你得把公式里的N乘上去

跪谢跪谢跪谢Sob前辈.

sobereva 发表于 Post on 2016-1-20 03:47:29

肖富贵发表于 2016-1-19 19:21
Sob前辈，您好.我试着修改Multiwfn中的某一行代码，来计算跃迁势. 修改的内容如下：

合理，但你得把公式里的N乘上去

肖富贵 发表于 Post on 2016-1-18 07:08:18

sobereva 发表于 2016-1-18 06:55
正确

跪谢跪谢跪谢Sob前辈.

sobereva 发表于 Post on 2016-1-18 06:55:51

肖富贵发表于 2016-1-18 01:56
Sob前辈，您好. 我看了您的帖子后，试着重新对此公式做了调整，不知道对与不对。请问您可不可以帮忙看一 ...

正确

肖富贵 发表于 Post on 2016-1-7 00:06:13

sobereva 发表于 2016-1-7 00:03
Multiwfn当前的代码直接给不出来跃迁势，但是改改代码能得到。
function.f90里的function eleesp是计算电 ...

跪谢跪谢跪谢.

sobereva 发表于 Post on 2016-1-7 00:03:14

Multiwfn当前的代码直接给不出来跃迁势，但是改改代码能得到。
function.f90里的function eleesp是计算电子密度对静电势的贡献的。计算的时候会循环每一对儿高斯函数(GTF)，循环中会看到注释：
!Now calc "term"=<psi(iprim)|1/r_Z|psi(jprim)>
也就是说下面的term变量就是<GTF_i|1/r_Z|GTF_j>这个积分。Multiwfn利用这个积分，乘上轨道对GTF的展开系数，乘上轨道占据数再累加，得到对静电势的贡献。对于获得跃迁势，在展开系数这块动动手脚改成对应跃迁密度矩阵的情况。

肖富贵 发表于 Post on 2016-1-6 23:34:09

不好意思，帖子中有两处乱码。第一处是phi，第二处是rho。从word中直接粘贴过来就乱码了，抱歉。

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