简单Multiwfn使用问题和波函数分析问题答疑专贴

sobereva · 发表于 Post on 2018-7-11 05:03:40

Adair 发表于 2018-7-10 17:18
各位坛友，老师:我在用Multiwfn绘制跃迁密度矩阵平面图（TDM）的时候，按照操作步骤，中间有一步提示“If s ...

看最新版手册3.21.2节的说明

Adair · 发表于 Post on 2018-7-11 10:12:57

sobereva 发表于 2018-7-11 05:03
看最新版手册3.21.2节的说明

谢谢老师

mly12345 · 发表于 Post on 2018-7-14 16:15:50

老师好，我在计算电子激发过程中片段简单电荷转移百分比的时候发现整个分子的各个原子一起计算总电荷量就是0，但是我分成各个片段后，能构成整个分子的所有的片段的电荷量加在一起总和并不是0，这是因为什么呢？

sobereva · 发表于 Post on 2018-7-15 05:52:30

mly12345 发表于 2018-7-14 16:15
老师好，我在计算电子激发过程中片段简单电荷转移百分比的时候发现整个分子的各个原子一起计算总电荷量就是 ...

我不清楚你具体怎么做的，最好能给个简单例子

mly12345 · 发表于 Post on 2018-7-15 12:33:36

sobereva 发表于 2018-7-15 05:52
我不清楚你具体怎么做的，最好能给个简单例子

具体情况是这样的，我计算的是一个分子“D-π-A-π-A-π-D”的电荷转移，计算整个分子的基态时候的总电荷量为0.00000286.当我分别计算了D，π，A部分的结构（比如D是一个甲氧基）之后，我将各个计算的电荷量加在一起并不是0.00000286了，这是我计算错误吗？

mly12345 · 发表于 Post on 2018-7-15 12:35:50

还有一个问题，就是我在计算分子中的某一部分比如苯环，激发态的电荷量减去基态的，结果是是正的，那是不是就可以说明我这个苯环是在激发后是获得电荷的呢？谢谢老师

平辉正 · 发表于 Post on 2018-7-16 09:46:06

Sob老师，我在使用multiwfn 加载VASP的CHGCAR转换成的cube文件后，发现有的文件能打开（2.cube），有的会闪退(1.cube)。我看过您关于cube文件格式的帖子，感觉格式没问题，不明白为什么会加载不了，特来向您求助。

闪退信息：
forrt1: severe(24):end-of-file during read, unit 10

附件：cube文件地址 https://pan.baidu.com/s/1hfH_QgttGXGe8-PeD5aHIw

sobereva · 发表于 Post on 2018-7-16 10:27:10

平辉正发表于 2018-7-16 09:46
Sob老师，我在使用multiwfn 加载VASP的CHGCAR转换成的cube文件后，发现有的文件能打开（2.cube），有的会闪 ...

上传cube文件前应当压缩一下，体积能小N倍

sobereva · 发表于 Post on 2018-7-16 10:27:46

mly12345 发表于 2018-7-15 12:35
还有一个问题，就是我在计算分子中的某一部分比如苯环，激发态的电荷量减去基态的，结果是是正的，那是不是 ...

是的（即丢失了电子）

sobereva · 发表于 Post on 2018-7-16 10:28:15

mly12345 发表于 2018-7-15 12:33
具体情况是这样的，我计算的是一个分子“D-π-A-π-A-π-D”的电荷转移，计算整个分子的基态时候的总电荷 ...

如果只是零点零零零几后面有差异，这属于数值积分误差，不用管

ABetaCarw · 发表于 Post on 2018-7-16 11:29:29

老师好，我用#p b3lyp/genecp（6-31G* Lanl2dz） scrf=smd TD opt  计算含Fe配合物激发态，主要是研究荧光猝灭。研究发现S1 振子强度为0，NTO变换后发现该跃迁主要是Fe的d-d跃迁。应该是传说的宇称禁阻。
但是因为b3lyp可能存在ghost态问题，所以想用Multiwfn验证一下，但是出现上述问题，请问该怎么办。
我看老师给的例子里面那个没显示的数字是 1st/2nd terms: 后面的数字之差。我的计算“ 1st/2nd terms:”后面的数字是不是大到不正常了？
还有一个量化的问题不知道能否在这个波函数专贴中提问。
S1激发涉及的几个轨道里面没有HOMO225-LUMO226 是否正常？是不是方法太差基组太烂定性错误了
Excited State 1:    Singlet-A    1.5878 eV  780.85 nm  f=0.0000  <S**2>=0.000
   221 -> 236       0.16662
   221 -> 238       0.13921
   223 -> 234       0.11646
   223 -> 236       0.49486
   223 -> 238       0.41335
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -2722.64080136
Copying the excited state density for this state as the 1-particle RhoCI density.

平辉正 · 发表于 Post on 2018-7-16 12:51:05

sobereva 发表于 2018-7-16 10:27
上传cube文件前应当压缩一下，体积能小N倍

已重新上传压缩包，https://pan.baidu.com/s/1hfH_QgttGXGe8-PeD5aHIw

qingqing1223 · 发表于 Post on 2018-7-16 18:18:46

请问社长，我在用Multiwfn计算HF与NH4阳离子的dimer体系进行原子盆积分的时候，在选择2 The same as 1, but with exact refinement of basin boundary然后进行电子密度的盆积分计算时，屏幕上出现好多Warning: Exceeded the step limit of steepest ascent process!，请问在这种情况下输出的盆积分结果还可靠吗？如果不可靠的话应该怎样调整参数设置呢？

mly12345 · 发表于 Post on 2018-7-17 08:33:33

sobereva 发表于 2018-7-16 10:28
如果只是零点零零零几后面有差异，这属于数值积分误差，不用管

老师，是这样的。我的D1/2-π1-A-π3-A-π2-D3/4结构，总的电荷量算完为=0.00000286
但是各的部分
D1=0.1744%
D2=0.0880%
D3=0.1744%
D4=0.0880%
π1=0.3292%
π2=0.3292%
π3=-0.8164%
A1=-1.1057%
A2=-1.1063%
这几个加在一起的电荷量=-1.8451%，为什么不是0呢？是不是我有什么没注意到的地方？
求老师指点一下~

sobereva · 发表于 Post on 2018-7-17 08:40:15

mly12345 发表于 2018-7-17 08:33
老师，是这样的。我的D1/2-π1-A-π3-A-π2-D3/4结构，总的电荷量算完为=0.00000286
但是各的部分
D1=0 ...

你最好把你的完整操作步骤贴一下

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

[Multiwfn使用咨询] 简单Multiwfn使用问题和波函数分析问题答疑专贴