24小时热门版块排行榜    

查看: 2452  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小范范1989

木虫 (著名写手)

[交流] 高斯优化最低三重态的问题。 已有1人参与

之前,我对于优化单重态和三重态有自己的理解,如下:
1:优化s1的结构的时候,
#p b3lyp/6-31g opt td(singlets,root=1)
这种方法应该是对的吧,同时是不是会默认计算三个td?那和下面的计算方法有什么区别呢?
#p b3lyp/6-31g opt td(nstates=6,singlets,root=1)
这个事计算6个td?和上一个,就这一点区别?
对于这两个的区别,我还在测试。
2:就是优化三重态的问题,这个是我的主要问题:如下
第一种
#P opt b3lyp/6-31g  td(triplets,root=1)
  0  1
坐标
我的问题是:此处的自旋多重读为1,但是td中要写triplets。 同时,坐标最好为s0的基态结构坐标、这样准确。是吧?
第二种
优化t1的第二种方式:
#P b3lyp/6-31g opt
0  3
坐标
这个自旋多重读为3,
我的问题是,我测试过这两种方式,得出来的HF能量是不同的,同时,几何结构也不相同,我想知道一下,三重态优化,我的方法对不对。谢谢
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 第一种.log
  • 2014-10-26 20:12:20, 1.51 M
  • 附件 2 : 第二种方法.log
  • 2014-10-26 20:12:44, 1.84 M

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

参考别人说的,我总结的,不知道是否正确,还请指正:
激发态的优化
一:优化s1
方法就一种: #p b3lyp/6-31g*  opt  td(singlets,root=1)
0        1
              坐标  (此处坐标为优化后s0的结构,即基态)
说明:
1:对于此计算,是一种TDDFT的方法,默认计算nstates=3,
   与td(nstates=6,singlets,root=1)没有区别。
2:计算可以加入溶剂:scrf=(solvent=n-hexane)。可以加入循环次次和收敛标准:scf(maxcyc=400,conver=7)
3:计算出错:You need to solve for more vectors in order to follow this state.可以在td中加入direct来解决。










第二:优化t1
第一种方法:#p b3lyp/6-31g* opt td(triplets,root=1)
0         1
             坐标  (为优化后s0的坐标,即基态)
1:这种方法是以基态单重态作为参考态,优化第一个激发态(三重态)。利用的TDDFT的方法。
2:默认计算nstates=3。


第二种方法:#P  b3lyp/6-31g*  opt  
0         3
             坐标 (最好也是s0的坐标,但自旋多重度不同)
1:这种方法是优化三重态“基态”。利用的为DFT的方法,
2:DFT和TDDFT计算的都是同一个真实的电子态。DFT计算的电子态在所有三重态中能量最低,所以是三重态的“基态”;此外,如果把单重态也考虑进去,由于有能量更低的单重态,以其为参考态,我们要研究的该电子态又是激发态。
3:虽然描述的是同样的电子态,但两种方法得到的结果(Hartree)并不相同!因为DFT和TD-DFT都是近似的方法,而且采用了不同的一系列近似,相当于使用两种不同的方法,得到的数值自然不会相同。

4:推荐选择第二种方法(DFT),即“基态”方法。使用基态方法可以优化所有轨道能级,而TD只能以基态作为参考态优化特定的电子态。

第三种方法:#P  b3lyp/6-31g*  opt  td(triplets,root=1)
0         3
             坐标
1:这种方式为,利用T1的结构优化T2的结构。
2:自旋多重度为0 3,那么算出来的波函数是T1的,在做triplets的TDDFT,就会的到更高的三重态,即T2(root=1)

3: 当用三重态作为参考态的时候,triplet选项无效。Gaussian的U-TDDFT不能计算自旋反转,因此所有激发态也都是三重态,这时候root=1对应T2,参考态才是T1。另外,你只能得到近似的三重态,因为自旋非限制DFT有自旋污染。
4: 所以说:TDDFT是以基态为基础计算激发态的。
不论激发态单线态还是三线态,对于基态是单线态的分子,都应该把自旋多重度设置为1。
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
5楼2014-11-08 16:55:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

高斯技术支持

新虫 (初入文坛)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小范范1989: 金币+6 2014-11-07 19:56:49
gmy1990: 金币+3 2014-11-07 23:32:25
1. NStates=M
计算M个激发态,默认M=3
#p b3lyp/6-31g opt td(singlets,root=1)
计算3个激发态

#p b3lyp/6-31g opt td(nstates=6,singlets,root=1)
计算6个激发态

您说的3个TD/6个TD不严谨。高斯只做一个TD计算,只是考虑的电子态数目不同。

2. 这个问题比较有意思,这里详细讲一下。

第一种
#p b3lyp/6-31g opt td(triplets,root=1)

Title

0 1

这种方法是以基态单重态作为参考态,优化第一个激发态(三重态)

第二种
#p b3lyp/6-31g opt

Title

0 3

这种方法是优化三重态“基态”

两种方式描述的都是同一个真实的电子态。该电子态在所有三重态中能量最低,所以是三重态的“基态”;此外,如果把单重态也考虑进去,由于有能量更低的单重态,以其为参考态,我们要研究的该电子态又是激发态。

虽然描述的是同样的电子态,但两种方法得到的结果(Hartree)并不相同!因为DFT和TD-DFT都是近似的方法,而且采用了不同的一系列近似,相当于使用两种不同的方法,得到的数值自然不会相同。

推荐选择第二种方法,即“基态”方法。使用基态方法可以优化所有轨道能级,而TD只能以基态作为参考态优化特定的电子态。事实上有文献试图通过限制对称性来使用基态方法来优化其他激发态。
2楼2014-11-07 12:19:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
2楼: Originally posted by 高斯技术支持 at 2014-11-07 12:19:05
1. NStates=M
计算M个激发态,默认M=3
#p b3lyp/6-31g opt td(singlets,root=1)
计算3个激发态

#p b3lyp/6-31g opt td(nstates=6,singlets,root=1)
计算6个激发态

您说的3个TD/6个TD不严谨。高斯只做一个 ...

谢谢你,讲解的很仔细。学习了,
看来你还真是高斯的官网小木虫??
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
3楼2014-11-07 19:57:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高斯技术支持

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 小范范1989 at 2014-11-07 06:57:15
谢谢你,讲解的很仔细。学习了,
看来你还真是高斯的官网小木虫??...

4楼2014-11-08 12:29:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见