计算化学公社

标题: 请教：金属卟啉吸收光谱计算的一些问题请教 [打印本页]

作者
Author: yakir 时间: 2019-8-26 22:47
标题: 请教：金属卟啉吸收光谱计算的一些问题请教
请教大家，我计算了一个金属卟啉结构，基态是单重态，计算后的结果如下图所示，请问：
1、前四个激发态的S^2都接近2，更接近三重态，难道是激发过程中电子发生了自旋反转？
2、是因为发生了自旋反转，所以前四个激发态的振子强度都为0？那这种情况应该如何计算呢？
3、第九激发态的S^2=-0.095，这个负值怎么理解
谢谢大家的解答

作者
Author: sobereva 时间: 2019-8-26 22:55
基态是单重态时，除非你用了U或者guess=read读了U的波函数（虽然不知道你是什么目的），否则不可能算成这样。当参考态闭壳层时，TDDFT算的激发态都是自旋纯态，根本没有这个问题

作者
Author: yakir 时间: 2019-8-27 09:14

sobereva 发表于 2019-8-26 22:55
基态是单重态时，除非你用了U或者guess=read读了U的波函数（虽然不知道你是什么目的），否则不可能算成这样 ...

老师我是这样算的，先做了一个单点计算，设置的是单重态，然后用guess=read读取的chk来计算吸收光谱

作者
Author: sobereva 时间: 2019-8-27 20:55

yakir 发表于 2019-8-27 09:14
老师我是这样算的，先做了一个单点计算，设置的是单重态，然后用guess=read读取的chk来计算吸收光谱

甭read。若还这样上传整个输入输出文件

作者
Author: yakir 时间: 2019-8-27 20:58

sobereva 发表于 2019-8-27 20:55
甭read。若还这样上传整个输入输出文件

好的，谢谢老师，我把guess=read去掉计算一下

作者
Author: yakir 时间: 2019-8-28 09:22

sobereva 发表于 2019-8-27 20:55
甭read。若还这样上传整个输入输出文件

老师您好，我去掉guess=read计算吸收光谱出现了激发能为负的情况，现在正在按照您的贴子去寻找稳定的波函数。另外我计算了一个双重态的结构，计算出的S^2很多为2.7的，难道是三重态？但是基态是双重态的结构激发态怎么会出现三重态呢？请老师指点，对应的输入输出文件我放在了这个回帖的下方

作者
Author: yakir 时间: 2019-8-28 09:24
对应的输入输出文件

作者
Author: biogon 时间: 2019-8-28 10:11

yakir 发表于 2019-8-28 09:22
老师您好，我去掉guess=read计算吸收光谱出现了激发能为负的情况，现在正在按照您的贴子去寻找稳定的波函 ...

你换自旋多重度看看，卟啉金属配合物一般基态不是单重态的

作者
Author: yakir 时间: 2019-8-28 10:15

biogon 发表于 2019-8-28 10:11
你换自旋多重度看看，卟啉金属配合物一般基态不是单重态的

哦哦，谢谢您，不过如果所有电子都配对的话，就算基态不是单重态那也不可能是双重态吧

作者
Author: biogon 时间: 2019-8-28 10:59

yakir 发表于 2019-8-28 10:15
哦哦，谢谢您，不过如果所有电子都配对的话，就算基态不是单重态那也不可能是双重态吧

二重态应该不会是

作者
Author: yakir 时间: 2019-8-28 11:04

biogon 发表于 2019-8-28 10:59
二重态应该不会是

嗯嗯，我去试试

作者
Author: sobereva 时间: 2019-8-28 23:22

yakir 发表于 2019-8-28 09:22
老师您好，我去掉guess=read计算吸收光谱出现了激发能为负的情况，现在正在按照您的贴子去寻找稳定的波函 ...

2.7的没法进行确切的归属，自旋污染太大，没法说自旋多重度

作者
Author: haokkt 时间: 2019-11-28 11:00
老师你好，金属卟啉怎么在生成脚本的时候会是双重态

作者
Author: sobereva 时间: 2019-11-29 09:55

haokkt 发表于 2019-11-28 11:00
老师你好，金属卟啉怎么在生成脚本的时候会是双重态

“生成脚本”语义不明
你设多少就是多少

欢迎光临计算化学公社 (http://bbs.keinsci.com/)