计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 22202 | 评论数: 21 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-11-16 15:29

正文摘要:

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

回复 Reply

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关键词?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 04:17 , Processed in 0.195023 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list