multi mrci计算问题，态平均？ - 量子化学 - Molpro/Molcas - 第2页小木虫论坛-学术科研互动平台

11楼2017-12-11 16:47:32

kent1022

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 566
散金: 20
红花: 1
帖子: 1278
在线: 397.8小时
虫号: 680880
注册: 2008-12-24
专业: 原子和分子物理

引用回帖:

11楼: Originally posted by 18766315007 at 2017-12-11 16:47:32
像occ和wf这样的设置也让我很费解，不太清楚？不知道哪里有讲的？

说明书上讲到了的

12楼2017-12-11 16:55:51

18766315007

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 6018.5
散金: 100
红花: 6
帖子: 1248
在线: 262.1小时
虫号: 2979660
注册: 2014-02-20
专业: 原子和分子物理

引用回帖:

12楼: Originally posted by kent1022 at 2017-12-11 16:55:51
说明书上讲到了的...

好的，谢谢！还有一个问题。就是OCC卡后面的轨道是怎么确定的？像同核双原子分子D2h点群，ag,b3u,b2u,b1g,b1u,b2g,b3g,au对应分子的什么轨道。也就是说，给定一个双原子分子的电子态，怎么确定occ卡后面的数值？

13楼2017-12-11 18:43:23

kent1022

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 566
散金: 20
红花: 1
帖子: 1278
在线: 397.8小时
虫号: 680880
注册: 2008-12-24
专业: 原子和分子物理

引用回帖:

13楼: Originally posted by 18766315007 at 2017-12-11 18:43:23
好的，谢谢！还有一个问题。就是OCC卡后面的轨道是怎么确定的？像同核双原子分子D2h点群，ag,b3u,b2u,b1g,b1u,b2g,b3g,au对应分子的什么轨道。也就是说，给定一个双原子分子的电子态，怎么确定occ卡后面的数值？...

让程序计算hf，输出里面就有默认的轨道

赞一下(1人)

» 本帖已获得的红花（最新10朵）

18766315007

14楼2017-12-11 21:59:53

18766315007

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 6018.5
散金: 100
红花: 6
帖子: 1248
在线: 262.1小时
虫号: 2979660
注册: 2014-02-20
专业: 原子和分子物理

送红花一朵

引用回帖:

14楼: Originally posted by kent1022 at 2017-12-11 21:59:53
让程序计算hf，输出里面就有默认的轨道...

谢谢！我试着用这种默认的轨道算两个态的跃迁偶极矩时，却是出现了错误：说是建议我使用，option, nstate=2；
然后我在某些轨道又加入活性空间，没有出现这种错误？这是怎么回事？这样算出来的结果合理吗？

还有，像你这个主题里提的问题：1Delta对应C2v中的A1+A2，是用态平均处理是可以的嘛？然后ci的计算是分开两个写，还是两个wf卡写在同一个ci下？

15楼2017-12-12 20:23:42

kent1022

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 566
散金: 20
红花: 1
帖子: 1278
在线: 397.8小时
虫号: 680880
注册: 2008-12-24
专业: 原子和分子物理

引用回帖:

15楼: Originally posted by 18766315007 at 2017-12-12 20:23:42
谢谢！我试着用这种默认的轨道算两个态的跃迁偶极矩时，却是出现了错误：说是建议我使用，option, nstate=2；
然后我在某些轨道又加入活性空间，没有出现这种错误？这是怎么回事？这样算出来的结果合理吗？

还 ...

是的，刚开始是活性空间不够
A1、A2对称性分开算，一个ci下对一个wf

» 本帖已获得的红花（最新10朵）

18766315007

16楼2017-12-12 22:07:54

18766315007

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 6018.5
散金: 100
红花: 6
帖子: 1248
在线: 262.1小时
虫号: 2979660
注册: 2014-02-20
专业: 原子和分子物理

送红花一朵

引用回帖:

16楼: Originally posted by kent1022 at 2017-12-12 22:07:54
是的，刚开始是活性空间不够
A1、A2对称性分开算，一个ci下对一个wf...

非常感谢您的耐心解答！还有一个问题，如果一个ci对应一个wf，而1Delta对应C2v中的A1+A2，两ci运算得到的结果如何处理才能表示1Delta态？

17楼2017-12-13 11:21:07

kent1022

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 566
散金: 20
红花: 1
帖子: 1278
在线: 397.8小时
虫号: 680880
注册: 2008-12-24
专业: 原子和分子物理

引用回帖:

17楼: Originally posted by 18766315007 at 2017-12-13 11:21:07
非常感谢您的耐心解答！还有一个问题，如果一个ci对应一个wf，而1Delta对应C2v中的A1+A2，两ci运算得到的结果如何处理才能表示1Delta态？...

A1A2都要算，得到的能量应该是简并的，就可以确定是Delta了

18楼2017-12-13 12:18:54

18766315007

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 6018.5
散金: 100
红花: 6
帖子: 1248
在线: 262.1小时
虫号: 2979660
注册: 2014-02-20
专业: 原子和分子物理

送红花一朵

你好！还想问您一个问题，我想算跃迁偶极矩，采用如下命令：

{ci;trans,7000.2,7001.2;}
Trdm1(i)=trdmx;
Trdm2(i)=trdmy;
Trdm3(i)=trdmz;

table,rn2,eci1,eci2,edi1,edi2,Trdm1,Trdm2,Trdm3;
但是输出的偶极矩和计算的ci；trans<1.6|DMX|1.7>的结果不一样，该如何修改才能输出正确的跃迁偶极矩？

19楼2017-12-23 10:14:48