sobereva 发表于 2019-8-3 23:53
density=current改为density=rhoci应当能完全解决问题
但此时用的是非弛豫密度
激发态的弛豫密度和非弛 ...

再次谢谢老师，问题已经完全解决。

sobereva 发表于 2019-8-3 23:53
density=current改为density=rhoci应当能完全解决问题
但此时用的是非弛豫密度
激发态的弛豫密度和非弛 ...

明白了，谢谢老师。我先补下知识再试

一未命名一 发表于 2019-8-3 17:45
老师，我更改文件关键词试过后发现，部分文件的报错问题的确解决了，但是还有几个还会出现同样问题，不知 ...

density=current改为density=rhoci应当能完全解决问题
但此时用的是非弛豫密度
激发态的弛豫密度和非弛豫密度什么关系在此文详细说了
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434（http://bbs.keinsci.com/thread-10775-1-1.html）

sobereva 发表于 2019-7-31 22:53
如果你用的是TD形式算的，改成TDA形式（把TD关键词直接替换成TDA），应当没这个问题

十分感谢老师，我去实际操作一下

如果你用的是TD形式算的，改成TDA形式（把TD关键词直接替换成TDA），应当没这个问题

granvia 发表于 2019-7-31 16:02
好像不行吧。

另外，激发态能用基于路易斯结构的共振理论感觉有点怪啊，因为路易斯结构好像是针对基 ...

懂了，十分感谢。我再去补补相关知识

一未命名一 发表于 2019-7-31 12:37
原来如此，太感谢了
那可以请教一下不让NBO程序自动用FIXDM修复密度矩阵的方法吗？

好像不行吧。

另外，激发态能用基于路易斯结构的共振理论感觉有点怪啊，因为路易斯结构好像是针对基态分子的吧？  NRT的框架是否适用于激发态我不知道，好像没见过有人论证过。其实NRT基于密度矩阵分解的ansatz也是饱受质疑（我个人认为在单行列式波函数下它根本就是错误的）。用VB方法似乎可以处理，但计算量太大，而且从原理上说，给出共振结构的成分也并非严格可靠，因为存在两个问题：1.所有可能的共振式之间并非完全独立，所以必须要去掉线性相关的共振结构，因此在取舍上就有人为的任意性（一般采用Rumer规则，但这也是有一定任意性的）。2. 各共振式之间并非彼此正交，往往重叠积分相当可观，于是就会面对如何合理划分重叠积分的问题，如Mulliken还是各改进的Mulliken方案？这又有任意性。

granvia 发表于 2019-7-31 11:21
因为激发态波函数不是单行列式的，采用的密度矩阵是广义的，所以在NBO算法中可能会导致轨道的占据数不合 ...

原来如此，太感谢了
那可以请教一下不让NBO程序自动用FIXDM修复密度矩阵的方法吗？

一未命名一 发表于 2019-7-31 10:32
十分感谢。
我没加FIXDM关键词出现这样的问题，还存在其他可能的原因吗？我是这样写的：
$GENNBO  NAT ...

因为激发态波函数不是单行列式的，采用的密度矩阵是广义的，所以在NBO算法中可能会导致轨道的占据数不合理，如负值或大于2。因此，NBO程序会自动用FIXDM来修复密度矩阵，使得轨道占据数在0和2之间。但这种修复后的密度矩阵并不等于原来的密度矩阵，而只是尽量接近，所以有可能就使得修复后的密度矩阵给出的总电子数（等于密度矩阵的trace）与原来的有偏差。

对于单行列式波函数，由于密度矩阵的幂等性，可以证明轨道占据数是严格在0和2之间的，所以不会存在这样的问题。

granvia 发表于 2019-7-30 23:51
试试去掉FIXDM选项

十分感谢。
我没加FIXDM关键词出现这样的问题，还存在其他可能的原因吗？我是这样写的：
$GENNBO  NATOMS=43  NBAS=706  UPPER  BODM  FORMAT=PRECISE  $END
$NBO  FILE=25o nrt memory=10gb $END

试试去掉FIXDM选项