计算化学公社

标题: 在开壳层体系中如何进行broken symmetry计算? [打印本页]

作者
Author:
司马打球    时间: 2015-10-31 10:23
标题: 在开壳层体系中如何进行broken symmetry计算?
在开壳层体系中如何进行broken symmetry计算?

作者
Author:
sobereva    时间: 2015-10-31 12:43
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
作者
Author:
司马打球    时间: 2015-11-4 19:13
sobereva 发表于 2015-10-31 12:43
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

已解决,多谢!
作者
Author:
司马打球    时间: 2016-1-7 17:11
sobereva 发表于 2015-10-31 12:43
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

您好!按照您之前写的片段组合波函数法得到的单线态双自由基计算结果仍然没有spin density,下面是我的输入文件,不知道有没问题?
%chk=singlet.chk
# opt ub3lyp/3-21g geom=connectivity guess=fragment=2

Title Card Required

0 1 0 3 0 -3
作者
Author:
sobereva    时间: 2016-1-7 22:16
司马打球 发表于 2016-1-7 17:11
您好!按照您之前写的片段组合波函数法得到的单线态双自由基计算结果仍然没有spin density,下面是我的输 ...


尝试UB3LYP/6-31G* guess(mix,always) opt
或者先用以上关键词算个单点,如果初始结构下都没有对称破缺,优化后也对称破缺不了。
作者
Author:
司马打球    时间: 2016-1-12 18:23
sobereva 发表于 2016-1-7 22:16
尝试UB3LYP/6-31G* guess(mix,always) opt
或者先用以上关键词算个单点,如果初始结构下都没有对称破 ...

直接算?不用片段组合?
作者
Author:
sobereva    时间: 2016-1-12 21:14
司马打球 发表于 2016-1-12 18:23
直接算?不用片段组合?

一般都没必要片段组合,直接用我说的关键词就能达到目的。如果发现没得到对称破缺态再考虑用略麻烦点的片段组合,或者stable=opt。
作者
Author:
司马打球    时间: 2016-1-13 08:57
sobereva 发表于 2016-1-12 21:14
一般都没必要片段组合,直接用我说的关键词就能达到目的。如果发现没得到对称破缺态再考虑用略麻烦点的片 ...

您好!按照UB3LYP/6-31G* guess(mix,always) opt关键词可以得到破缺态,但是为什么直接用guess=mix却不能得到破缺态?非常感谢您的建议!
作者
Author:
sobereva    时间: 2016-1-13 17:11
司马打球 发表于 2016-1-13 08:57
您好!按照UB3LYP/6-31G* guess(mix,always) opt关键词可以得到破缺态,但是为什么直接用guess=mix却不能 ...


你不用always的话,一开始如果是闭壳层状态,优化时候每一步初猜都会用上一步收敛的波函数,因此闭壳层状态会保持下去。always相当于每一步都产生对称破缺初猜,到了该出现对称破缺结构的时候就自然而然对称破缺了。
U也是必须用的。

作者
Author:
guoy14iccas    时间: 2017-7-11 14:28
sobereva 发表于 2016-1-7 22:16
尝试UB3LYP/6-31G* guess(mix,always) opt
或者先用以上关键词算个单点,如果初始结构下都没有对称破 ...

sob老师,请教您个问题,我要计算开壳层单线态能量,然后计算自由基特性y;在计算能量的时候 我写的关键词是 #p td UCAM-B3LYP/6-31G** guess(mix,only) pop=full,但是会报错,去掉td就可以了。这样算出来的还是单线态的能量吗?
作者
Author:
sobereva    时间: 2017-7-11 16:10
guoy14iccas 发表于 2017-7-11 14:28
sob老师,请教您个问题,我要计算开壳层单线态能量,然后计算自由基特性y;在计算能量的时候 我写的关键词 ...


干嘛要写td?你又不做电子激发计算
干嘛要写only?光产生个初猜波函数有什么意义?

作者
Author:
guoy14iccas    时间: 2017-7-11 16:36
sobereva 发表于 2017-7-11 16:10
干嘛要写td?你又不做电子激发计算
干嘛要写only?光产生个初猜波函数有什么意义?

sob老师,我这样修改之后#p UCAM-B3LYP/6-31G** guess(mix,always) pop=full 这样算出来的能量是基态能量,还是Open-shell Singlet 能量?第一次弄这个,没有经验,非常感谢您的指导。
作者
Author:
sobereva    时间: 2017-7-11 21:51
guoy14iccas 发表于 2017-7-11 16:36
sob老师,我这样修改之后#p UCAM-B3LYP/6-31G** guess(mix,always) pop=full 这样算出来的能量是基态能量 ...

计算自旋极化单重态的关键词怎么写这里说了
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
作者
Author:
d1490541507    时间: 2020-10-23 09:42
guoy14iccas 发表于 2017-7-11 14:28
sob老师,请教您个问题,我要计算开壳层单线态能量,然后计算自由基特性y;在计算能量的时候 我写的关键词 ...

前辈,我也要算这个开壳层特性,有偿求助,可以教教我吗QQ1515274108
作者
Author:
有偿求指导    时间: 2020-10-23 10:01
guoy14iccas 发表于 2017-7-11 16:36
sob老师,我这样修改之后#p UCAM-B3LYP/6-31G** guess(mix,always) pop=full 这样算出来的能量是基态能量 ...

前辈,计算完开壳层单重态后,那个自由基特性y是怎么找到的呀
作者
Author:
sobereva    时间: 2020-10-23 10:43
有偿求指导 发表于 2020-10-23 10:01
前辈,计算完开壳层单重态后,那个自由基特性y是怎么找到的呀

CASSCF计算双自由基以及双自由基特征的计算
http://sobereva.com/264http://bbs.keinsci.com/thread-322-1-1.html
作者
Author:
sobereva    时间: 2020-10-23 10:44
d1490541507 发表于 2020-10-23 09:42
前辈,我也要算这个开壳层特性,有偿求助,可以教教我吗QQ1515274108

把我给的博文里的内容都好好仔细看了,有具体问题直接发帖问,别遇到点问题就试图靠钱来解决




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3