用高斯计算电荷转移及自然跃迁轨道问题求助

Jameslzy · 发表于 Post on 2017-11-16 15:29:26

各位老师好，最近想重现一个文献中的一些计算结果，掌握一下得到这些数据的方法。在咱们论坛上找了一些相关帖子，学习到了很多知识，可是有些问题还是没能解决，想在这里请教下各位老师

这个是我想重现的文献中的计算内容，目前有以下几个问题：
1）（a）图中的结构优化我已经完成了，但是不知道这个电荷转移量是如何计算出来的。之前看了一个帖子，里面sob老师说是把整个体系分成两个片段，将转移前后片段上的电荷加和，看电荷变化量。文章中计算方法里写到是用NBO计算的电荷转移量。我想的是对S0构型下的结构算一个NBO，再对S0构型同时算TDDFT和NBO比较二者差别，但最后得出的结果是二者的NBO没有任何差别，不知道是不是我的想法有误，下面我上传的是我的输入和输出文件，还请各位老师指点下问题在哪里。对于（a）图我还有一个问题就是计算方法中没有写明是用b3lyp-D还是wb97x做的NBO那么一般情况下认为是用的那个方法做的NBO呢？

a图的计算及输出文件.rar (531.11 KB, 下载次数 Times of downloads: 67)

2）经过最近在论坛上的学习，大概掌握了些计算自然跃迁轨道的流程，能够做出图（d）所示的轨道图，但是对图（d）中写的 hole-particle 对激发的贡献那个百分比数值是如何得到的还没有掌握，希望各位老师能指点一下。还有就是如何能看出各个状态对应的是LE，CT 还是两者混合的状态呢，希望各位老师帮忙解答一下。下面我上传的是我计算NTO的输入和输出文件。

NTO的输入输出文件.rar (87.17 KB, 下载次数 Times of downloads: 71)

3）第三个问题是个小问题，最近正在学习sob老师的Multiwfn软件，不知道这个软件有没有复制粘贴功能，因为每次打开文件要输入一个文件路径，如果能复制粘贴文件路径的话感觉会比较省时

第一次发帖，不是很熟练，有什么不妥的地方还希望各位老师能够指出来，我下次一定注意和改正。多谢各位老师

Jameslzy · 发表于 Post on 2017-11-21 18:23:51

让你变成回忆发表于 2017-11-20 20:57
额，我不是老师。也是一个刚接触量化不久的学生~~~~

喔喔喔，好的，谢谢同学

以后还要向同学你多学习呀

让你变成回忆 · 发表于 Post on 2017-11-20 20:57:53

Jameslzy 发表于 2017-11-20 19:31
嗯嗯好的，谢谢老师的解答，这样我对这方面又有了更深的了解

额，我不是老师。也是一个刚接触量化不久的学生~~~~

Jameslzy · 发表于 Post on 2017-11-20 19:31:55

让你变成回忆发表于 2017-11-20 08:27
我去核实了一下，楼下有人说的是对的。
仅仅写density关键词就行：
对于后HF，只有写了density关键词传 ...

嗯嗯好的，谢谢老师的解答，这样我对这方面又有了更深的了解

让你变成回忆 · 发表于 Post on 2017-11-20 08:27:00

Jameslzy 发表于 2017-11-16 19:58
十分感谢您的回复
1）我也觉得是自己的输入文件有问题，之前一直做的基态NBO，可能激发态的 ...

我去核实了一下，楼下有人说的是对的。
仅仅写density关键词就行：
对于后HF，只有写了density关键词传递给NBO模块的才是当前级别的密度矩阵，否则，传递的是HF级别的，分析结果也是HF级别的。
对于CIS、TDDFT，只有写了density关键词，传递给NBO模块的才是用root指定的激发态的密度矩阵，否则仅仅是把参考态的密度矩阵传递给了NBO，这样分析结果是基态的。

Jameslzy · 发表于 Post on 2017-11-18 22:22:37

zjxitcc 发表于 2017-11-18 22:00
sorry， TDDFT研究CT我也不是很清楚···没做过···你还是反复阅读一下sob的文章吧

哦哦好的，谢谢老师，我再找找sob老师的文章好好看一下

zjxitcc · 发表于 Post on 2017-11-18 22:00:16

Jameslzy 发表于 2017-11-18 09:03
恩恩，十分感谢老师您的回答。确实我用的方法和文献中不太一样，等一下用wb97x再试一下老师 ...

sorry， TDDFT研究CT我也不是很清楚···没做过···你还是反复阅读一下sob的文章吧

Jameslzy · 发表于 Post on 2017-11-18 09:03:53

zjxitcc 发表于 2017-11-17 23:34
只写density。你这写的是对的，与文献不同的话，可能是（1）我看文献是wb97x，但你用的是wb97xd；（2）你写 ...

恩恩，十分感谢老师您的回答。确实我用的方法和文献中不太一样，等一下用wb97x再试一下

老师，文献中算出来的是0.05e的转移，而用我这个方法算出来的是0.94e的转移，感觉这样就和文献中说的只是弱的电荷转移应该是对不上了吧？0.94e是不是算比较强的电荷转移了呢？还有个问题想问老师就是，在处理这种情况下是不是wb97xd比wb97x更好一点呀，因为之前看sob老师的文章中好像说加上色散校正更好，因为想看看以后自己做体系的时候用什么方法更可靠，所以想问您一下

谢谢老师。

zjxitcc · 发表于 Post on 2017-11-17 23:34:29

只写density。你这写的是对的，与文献不同的话，可能是（1）我看文献是wb97x，但你用的是wb97xd；（2）你写的态数15不一定与文献作者计算时一致，但这一般不会导致定性错误，只是结果可能定量上有点不同。

Jameslzy · 发表于 Post on 2017-11-17 21:27:45

zjxitcc 发表于 2017-11-17 14:36
只写density试试看

老师，我加上density关键词后确实得到了不同的数值，但是和文献报道上的差的好像很多。我把输入文件的命令行写出来，您看看有没有问题

基态的输入文件
# wb97xd/6-31g** scrf(solvent=acetonitrile) pop(reg,nboread)

激发态的输入文件
# wb97xd/6-31g** td(nstates=15) density
scrf(solvent=acetonitrile) pop(reg,nboread)

老师您看写的有问题吗，谢谢老师

Jameslzy · 发表于 Post on 2017-11-17 16:02:23

zjxitcc 发表于 2017-11-17 14:36
只写density试试看

density=1是吗，我试一下

zjxitcc · 发表于 Post on 2017-11-17 14:36:34

只写density试试看

Jameslzy · 发表于 Post on 2017-11-17 11:18:08

zjxitcc 发表于 2017-11-17 09:54
有没加density关键词？

老师您好，我加上了density=transition=1关键词来算S1态，结果还是和基态的NBO结果相同，您看这可能是什么情况呢？

Jameslzy · 发表于 Post on 2017-11-17 10:31:27

zjxitcc 发表于 2017-11-17 09:54
有没加density关键词？

喔喔喔，谢谢您的提醒，我算nto的时候加上了，刚才看了一下，nbo里面没有加，等下加上试试。十分感谢

zjxitcc · 发表于 Post on 2017-11-17 09:54:28

有没加density关键词？

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

[新手求助] 用高斯计算电荷转移及自然跃迁轨道问题求助

评分 Rate

浏览过的版块