24小时热门版块排行榜    

查看: 2880  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 cHIDTWCTtLlG 2026-08-19 5/250 2026-08-19 17:30 by oCMhZAsgNFoO
[基金申请] filecode,4个jtjc了 +10 ziyangfang 2026-08-19 11/550 2026-08-19 17:19 by wdc28
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 117/5850 2026-08-19 17:12 by 医学老男孩
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见