24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2500  |  回复: 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的回帖

高斯技术支持

新虫 (初入文坛)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+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的回帖

小范范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的回帖

高斯技术支持

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by 小范范1989 at 2014-11-08 03:55:10
参考别人说的,我总结的,不知道是否正确,还请指正:
激发态的优化
一:优化s1
方法就一种: #p b3lyp/6-31g*  opt  td(singlets,root=1)
0        1
              坐标  (此处坐标为优化后s0的结构,即基态)
说 ...

您给出的”三种方法“得到的是完全不同的电子态。 请先选择您需要计算的体系。
6楼2014-11-27 07:46:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小范范1989 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 085600,320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 278求调剂 +3 依旧! 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 309分085801求调剂 +11 MY_angel 2026-03-31 11/550 2026-04-04 19:11 by 蓝云思雨
[考研] 一志愿085404,总分291,四级已过,求调剂 +5 阿俊阿俊阿俊 2026-04-04 7/350 2026-04-04 13:23 by 莲菜就是藕吧
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
信息提示
请填处理意见