应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 这是高斯09计算圆锥交叉点的算例，求解关键字
121/2返回列表12下一页
查看: 3163  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dkstatsra

新虫 (初入文坛)

[求助] 这是高斯09计算圆锥交叉点的算例，求解关键字 已有1人参与

算例贴出来：
%chk=test290
#p casscf(4,4)/sto-3g opt=conical test guess=alter nosymm

Gaussian Test Job 290 (Part 1):
Optimization of the conical Intersection

0 1
  X
  X  1     A00
  X  1  1.0000  2   90.00
  C  1     A01  3   90.00  2     C00  0
  C  1     A01  3   90.00  4  180.00  0
  X  2  1.0000  1   90.00  3     C01  0
  C  2     A02  6   90.00  1     C02  0
  C  2     A02  6   90.00  7  180.00  0
  H  4     A03  1     B00  2     C03  0
  H  4     A04  1     B01  2     C04  0
  H  5     A05  1     B02  2     C05  0
  H  5     A06  1     B03  2     C06  0
  H  7     A07  2     B04  1     C07  0
  H  7     A08  2     B05  1     C08  0
  H  8     A09  2     B06  1     C09  0
  H  8     A10  2     B07  1     C10  0

     A00  2.168596
     A01   .727272
     A02   .727258
     A03  1.078159
     A04  1.078218
     A05  1.070725
     A06  1.070723
     A07  1.078424
     A08  1.078441
     A09  1.070290
     A10  1.070255
     B00  119.162
     B01  119.221
     B02  120.718
     B03  120.712
     B04  119.726
     B05  119.771
     B06  120.691
     B07  120.710
     C00   68.814
     C01   -1.316
     C02   68.995
     C03  107.509
     C04 -107.700
     C05   92.605
     C06  -91.608
     C07  106.076
     C08 -106.377
     C09   93.573
     C10  -92.489

12 15

0.5       0.5


最末尾的这个“12 15”和“0.5 0.5”分别都是什么意思呢？
回复此楼

» 收录本帖的淘帖专辑推荐

经验

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-26 12:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
opt=conical时是以态平均方式做CASSCF，态平均计算中的态的权重是自己指定的，0.5 0.5就是指两个态权重相同，各占一半，这是常规做法。
赞一下(1人) 回复此楼
4楼2014-06-26 14:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dkstatsra

新虫 (初入文坛)
哦，对了，这个是test290中的第一部分
赞一下 回复此楼
2楼2014-06-26 12:07:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dkstatsra: 金币+5, ★★★很有帮助 2014-06-26 15:22:43
gmy1990: 金币+2 2014-06-27 03:15:32
12 15
意思是交换12和15号轨道，因为用了guess=alter

0.5  0.5
这个我觉得可能是opt=conical的附加输入，用来输入各个态的权重，0.5相当于态平均计算，其实应该不写也不要紧的。因自己很少用casscf来算，这个我也不是很确定，期待其他高手解答。
赞一下 回复此楼
3楼2014-06-26 14:03:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 枪下游魂 at 2014-06-26 14:03:32
12 15
意思是交换12和15号轨道，因为用了guess=alter

0.5  0.5
这个我觉得可能是opt=conical的附加输入，用来输入各个态的权重，0.5相当于态平均计算，其实应该不写也不要紧的。因自己很少用casscf来算，这个我 ...

谢谢您的指导，我还想问下，从这个例子的输入文件能否看出来这算的是哪两个态之间的圆锥交叉点？如果我想优化出S1态和T2态之间的圆锥交叉点我该怎么设置？
赞一下 回复此楼
5楼2014-06-26 15:33:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by sobereva at 2014-06-26 14:16:40
opt=conical时是以态平均方式做CASSCF，态平均计算中的态的权重是自己指定的，0.5 0.5就是指两个态权重相同，各占一半，这是常规做法。

谢谢您，怎么给不了金币了。。。。
赞一下 回复此楼
6楼2014-06-26 15:35:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)
引用回帖:
5楼: Originally posted by dkstatsra at 2014-06-26 15:33:27
谢谢您的指导，我还想问下，从这个例子的输入文件能否看出来这算的是哪两个态之间的圆锥交叉点？如果我想优化出S1态和T2态之间的圆锥交叉点我该怎么设置？...

例子里的应该是基态的单重态和三重态
激发态，结合手册和我的理解，可能是加nroot=2，nroot=3，但我没尝试过，还需要实践证实。
赞一下 回复此楼
7楼2014-06-27 08:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by 枪下游魂 at 2014-06-27 08:14:16
例子里的应该是基态的单重态和三重态
激发态，结合手册和我的理解，可能是加nroot=2，nroot=3，但我没尝试过，还需要实践证实。...

好的，谢谢您。我想知道我怎么看算完的结果是哪两个态之间的圆锥交叉点呢？
赞一下 回复此楼
8楼2014-06-27 09:17:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by sobereva at 2014-06-26 14:16:40
opt=conical时是以态平均方式做CASSCF，态平均计算中的态的权重是自己指定的，0.5 0.5就是指两个态权重相同，各占一半，这是常规做法。

请问下，我在输入文件里怎么指定我要计算的圆锥交叉点是S1态和T2态之间的呢？谢谢了
赞一下 回复此楼
9楼2014-06-27 09:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)
引用回帖:
8楼: Originally posted by dkstatsra at 2014-06-27 09:17:31
好的，谢谢您。我想知道我怎么看算完的结果是哪两个态之间的圆锥交叉点呢？...

预测圆锥交点。在执行路径使用选项Opt=Conical，把使用CASSCF 对指定态进行优化
的计算变成寻找这个态的圆锥交点或避免交叉点的计算。优化的结构将是在圆锥交点或避免
交叉点处的结构。检查CASSCF 输出文件中最后一步优化计算的本征值（在结构优化之前），
可以区分这两种可能结构：
FINAL EIGENVALUES AND EIGENVECTORS
VECTOR EIGENVALUES CORRESPONDING EIGENVECTOR
态 能量
1 -154.0503161  0.72053292 -0.48879229 ...
                          -0.16028934E-02 0.31874441E-02 ...
2 -154.0501151  0.45467877 0.77417416 ...
如果这两个本征值（在标有态编号行的第一项）几乎相同，则表示这两个态的能量相同，
为一圆锥交点。否则为避免交叉点。


以上节选自高斯手册
赞一下 回复此楼
10楼2014-06-27 09:49:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dkstatsra 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:47 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
关闭
请填处理意见
关闭