b******h
发帖数: 2715 | 1 借这里的人气问一下有没有Guassian里计算excited states 的模板？ 我
看了很多forum，可是还是不能写好syntext，我目前的方法是ub3lyp/6-31g(d,p) cis.
我的有机分子是个自由基，希望算出产生这个自由基需要的光波波长。
如果有大牛可以帮到我请电邮联系我也行，可能需要看具体的chemdraw分子。
再谢过 | b******h
发帖数: 2715 | | y****z
发帖数: 502 | 3 还是看Gaussian的说明书比较好。最简单的：
#p td=(singlets,nstate=20) b3lyp/6-311++G(d,p) 5D 7F
cis.
【在 b******h 的大作中提到】
: 借这里的人气问一下有没有Guassian里计算excited states 的模板？ 我
: 看了很多forum，可是还是不能写好syntext，我目前的方法是ub3lyp/6-31g(d,p) cis.
: 我的有机分子是个自由基，希望算出产生这个自由基需要的光波波长。
: 如果有大牛可以帮到我请电邮联系我也行，可能需要看具体的chemdraw分子。
: 再谢过
| d****n
发帖数: 397 | 4 不知道你问的问题是指什么。
A(molecule)=>B(radical)+C(radical) 这个photodissociation过程需要的光波波长？
还是这个自由基能吸收的光波波长？
有可能你说的自由基B,是这样产生的E(radical)+F(molecule)=>B(radical)。所以，没
有上下文，“算出产生这个自由基需要的光波波长”，不知道怎么算。
cis.
【在 b******h 的大作中提到】
: 借这里的人气问一下有没有Guassian里计算excited states 的模板？ 我
: 看了很多forum，可是还是不能写好syntext，我目前的方法是ub3lyp/6-31g(d,p) cis.
: 我的有机分子是个自由基，希望算出产生这个自由基需要的光波波长。
: 如果有大牛可以帮到我请电邮联系我也行，可能需要看具体的chemdraw分子。
: 再谢过
| b******h
发帖数: 2715 | 5 dragon:
A(molecule)=>B(radical)+C(radical) 这个photodissociation过程需要的光波波长。
就是这个模式
因为是新手，现在和老板搞僵，连和老板合作的老板也不给我指导了。
你有没有template可供我参考呢？ Gaussian的CIS, TD,都没有input模板，我想这也是
那个合作老板也迟迟不摧我计算这个激发态能量的原因，他可能也还没有figure out
吧。
谢谢大家 | d****n
发帖数: 397 | 6 你要先知道基态的能量，还有那个dissociation 对应的是哪个激发态。
这样先做对基态做一个构型优化，得到基态的能量，
然后，做TD,或者CIS,就能算出vertical excitation所需要的波长。
这个其实不是所需要的最短波长，可是根据Frank-Condon原理，是吸收的cross-
section最大的波长，也就是photodissociation所需要的波长。
你查查文献，看看那个photodissociation对应的是哪个excited state,假设对应的是
第一个excited state
那CIS的输入就是：
%chk=test.chk
#p cis uhf/basis gfinput geom=check guess=read
photodissociation
0,2
还有写gaussian输入时候，最后一行一定要是空行。
因为你的是open shell (radical) system,那些singlet,triplet keywords你不用关心。
如果用TDDFT的话，应该是#p td ub3lyp/basis geom=check guess=read gfinput
其中，geom用的是你优化的基态构型，初始轨道猜测是你基态的优化好的初始轨道猜测。
当然，你可以一次多算几个激发态，可以改TD,或者CIS的subkeywords,如果有实验跟你
的计算结果比较，就更好了。
你也可以用CASSCF试试，不过这个比上面的TD,CIS都要复杂些。
【在 b******h 的大作中提到】
: dragon:
: A(molecule)=>B(radical)+C(radical) 这个photodissociation过程需要的光波波长。
: 就是这个模式
: 因为是新手，现在和老板搞僵，连和老板合作的老板也不给我指导了。
: 你有没有template可供我参考呢？ Gaussian的CIS, TD,都没有input模板，我想这也是
: 那个合作老板也迟迟不摧我计算这个激发态能量的原因，他可能也还没有figure out
: 吧。
: 谢谢大家
| f**********e
发帖数: 1994 | 7 CASSCF 的能量通常不太准，如果开错 active space, 那结果会很操蛋。
我以前这样干过：先用 CI 扫优化构型附近,看看 low-lying 的 excited
states 是不是都是一个 configuration 独大的，如果是，我就把 Hartree
Fock 的 occupied 换一换，就变成 excited state 了。那时你想用
CCSDT 高阶方法还是优化 excited state 构型随便你。
如果 excited state 是一堆 Slater det 混出来的... 那就祝好运了。
【在 d****n 的大作中提到】
: 你要先知道基态的能量，还有那个dissociation 对应的是哪个激发态。
: 这样先做对基态做一个构型优化，得到基态的能量，
: 然后，做TD,或者CIS,就能算出vertical excitation所需要的波长。
: 这个其实不是所需要的最短波长，可是根据Frank-Condon原理，是吸收的cross-
: section最大的波长，也就是photodissociation所需要的波长。
: 你查查文献，看看那个photodissociation对应的是哪个excited state,假设对应的是
: 第一个excited state
: 那CIS的输入就是：
: %chk=test.chk
: #p cis uhf/basis gfinput geom=check guess=read
| b******h
发帖数: 2715 | 8 谢谢dragon，感动的流泪中。等会儿就上机实践。
【在 d****n 的大作中提到】
: 你要先知道基态的能量，还有那个dissociation 对应的是哪个激发态。
: 这样先做对基态做一个构型优化，得到基态的能量，
: 然后，做TD,或者CIS,就能算出vertical excitation所需要的波长。
: 这个其实不是所需要的最短波长，可是根据Frank-Condon原理，是吸收的cross-
: section最大的波长，也就是photodissociation所需要的波长。
: 你查查文献，看看那个photodissociation对应的是哪个excited state,假设对应的是
: 第一个excited state
: 那CIS的输入就是：
: %chk=test.chk
: #p cis uhf/basis gfinput geom=check guess=read
| b******h
发帖数: 2715 | 9 谢谢ferdinandshe，我一定好好做，把这个问题解决了。
【在 f**********e 的大作中提到】
: CASSCF 的能量通常不太准，如果开错 active space, 那结果会很操蛋。
: 我以前这样干过：先用 CI 扫优化构型附近,看看 low-lying 的 excited
: states 是不是都是一个 configuration 独大的，如果是，我就把 Hartree
: Fock 的 occupied 换一换，就变成 excited state 了。那时你想用
: CCSDT 高阶方法还是优化 excited state 构型随便你。
: 如果 excited state 是一堆 Slater det 混出来的... 那就祝好运了。
|
|
