24小时热门版块排行榜    

查看: 931  |  回复: 2

思雨G十年

金虫 (正式写手)

[求助] Molcas 高对称性分子 CASSCF方法下输入文件如何写

如题,最近被问到多组态效应,因而想用Molcas做一些CASSCF,但苦于无半点知识,不得不求教于大家。一下是我写的输出文件,结构是c3v对称性的,但是Molcas做一步SCF后给出的总电子数是108. 不知道怎么导致的。请高手赐教,在此先谢过。
  &GATEWAY
Title = P4Li
Symmetry
   XYZ
Basis set
LI.LITHIUM    / inline
          3.   2
* S-type functions
    11    4
            900.4600000        
            134.4330000        
             30.4365000        
              8.6263900        
              2.4833200        
              0.3031790        
              4.8689000        
              0.8569240        
              0.2432270        
              0.0635070        
              0.0243683        
      0.00228704             0.0000000              0.0000000              0.0000000        
      0.0176350              0.0000000              0.0000000              0.0000000        
      0.0873434              0.0000000              0.0000000              0.0000000        
      0.2809770              0.0000000              0.0000000              0.0000000        
      0.6587410              0.0000000              0.0000000              0.0000000        
      0.1187120              0.0000000              0.0000000              0.0000000        
      0.0000000              0.0933293              0.0000000              0.0000000        
      0.0000000              0.9430450              0.0000000              0.0000000        
      0.0000000             -0.00279827             0.0000000              0.0000000        
      0.0000000              0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000        
* P-type functions
    5    3
              4.8689000        
              0.8569240        
              0.2432270        
              0.0635070        
              0.0243683        
      0.0327661              0.0000000              0.0000000        
      0.1597920              0.0000000              0.0000000        
      0.8856670              0.0000000              0.0000000        
      0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              1.0000000        
* D-type functions
    1    1
              0.2000000        
              1.0000000        

Li1          0.000000    0.000000    0.131110
End of basis
Basis set
P.PHOSPHOROUS    / inline
          15.   2
* S-type functions
    13    6
          77492.4000000        
          11605.8000000        
           2645.9600000        
            754.9760000        
            248.7550000        
             91.1565000        
             91.1565000        
             36.2257000        
             15.2113000        
              4.7941700        
              1.8079300        
              0.3568160        
              0.1147830        
      0.0007810              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0060680              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0311600              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.1234310              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.3782090              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.5632620              0.0000000              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.1602550              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.6276470              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.2638490              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              1.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              0.0000000              0.0000000              1.0000000        
* P-type functions
    9    5
            384.8430000        
             90.5521000        
             29.1339000        
             10.8862000        
              4.3525900        
              1.7770600        
              0.6970050        
              0.2535320        
              0.0684930        
      0.0092060              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0698740              0.0000000              0.0000000              0.0000000              0.0000000        
      0.2924700              0.0000000              0.0000000              0.0000000              0.0000000        
      0.7281030              0.0000000              0.0000000              0.0000000              0.0000000        
      0.0000000              0.6283490              0.0000000              0.0000000              0.0000000        
      0.0000000              0.4280440              0.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              1.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              0.0000000              1.0000000        
* D-type functions
    1    1
              0.5500000        
              1.0000000        

P2           0.000000    0.000000    4.294807
P3           1.282336   -0.000010    2.463002
P4          -0.641159    1.110541    2.463002
P5          -0.641177   -1.110531    2.463002
End of basis
>>> Do while
  &SEWARD
>>> If ( Iter = 1 ) <<<
  &SCF
>>> EndIF <<<
  &RASSCF
    Lumorb
    spin= 2; nActEl= 13 0 0; Inacitve= 25; Ras2= 7
  &ALASKA
>>> EndDo
回复此楼
要以感性的思想和认识,理性的思考和行动对待科研。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
思雨G十年: 金币+50, 有帮助, The total number of the electron should be 63. I had maken it. Thanks a lot. 2013-10-25 07:26:31
gkf高: 金币+2 2013-10-25 09:22:04
总电子数是108?这个数字太神奇了
你说的对称性为C3V,但是定义的是Ci对称,由于不满足对称要求,我觉得得出来的电子数应该是(15*4+3)*2=126.
2楼2013-10-24 22:28:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

思雨G十年

金虫 (正式写手)

我计算得到的正确结果的输入:
1)
&gateway
Title= P4Li
coord = 4-31G.xyz
basis = 3-21G
group =c1
>>> Do while <<<
&SEWARD
&RASSCF
LUMORB
spin= 2; nactel = 13 0 0; inactive = 25; Ras2 = 13
&CASPT2; Frozen = 20
&SLAPAF
>>> EndDo <<<


4-31G.xyz
5

Li                 0.00000000    0.00000000   -2.75270300 Angstrom
P                  0.00000000    0.00000000    1.49928400 Angstrom
P                  0.00000000    1.28010200   -0.31624800 Angstrom
P                 -1.10860100   -0.64005100   -0.31624800 Angstrom
P                  1.10860100   -0.64005100   -0.31624800 Angstrom
2)
> export MOLCASDISK = 20480 MB
> export MOLCASMEM = 15000 MB
&GATEWAY
  Title
  P4LI

Symmetry
    X
Basis set
LI.LITHIUM    / inline
          3.   2
* S-type functions
    11    5
           1469.0000000        
            220.5000000        
             50.2600000        
             14.2400000        
              4.5810000        
              1.5800000        
              0.5640000        
              0.0734500        
              0.0280500        
              0.9060000        
              0.0086400        
      0.0007660             -0.0001200              0.0000000              0.0000000              0.0000000        
      0.0058920             -0.0009230              0.0000000              0.0000000              0.0000000        
      0.0296710             -0.0046890              0.0000000              0.0000000              0.0000000        
      0.1091800             -0.0176820              0.0000000              0.0000000              0.0000000        
      0.2827890             -0.0489020              0.0000000              0.0000000              0.0000000        
      0.4531230             -0.0960090              0.0000000              0.0000000              0.0000000        
      0.2747740             -0.1363800              0.0000000              0.0000000              0.0000000        
      0.0097510              0.5751020              0.0000000              0.0000000              0.0000000        
     -0.0031800              0.5176610              1.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              0.0000000              1.0000000        
* P-type functions
    6    4
              1.5340000        
              0.2749000        
              0.0736200        
              0.0240300        
              3.2420000        
              0.0057900        
      0.0227840              0.0000000              0.0000000              0.0000000        
      0.1391070              0.0000000              0.0000000              0.0000000        
      0.5003750              0.0000000              0.0000000              0.0000000        
      0.5084740              1.0000000              0.0000000              0.0000000        
      0.0000000              0.0000000              1.0000000              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000        
* D-type functions
    2    2
              0.1239000        
              0.0725000        
      1.0000000              0.0000000        
      0.0000000              1.0000000   
Li1          0.00000000    0.00000000    2.65000000 Angstrom
End of basis
Basis set
P.PHOSPHOROUS    / inline
          15.   1
* S-type functions
    16    4
          19413.3000000        
           2909.4200000        
            661.3640000        
            185.7590000        
             59.1943000        
             20.0310000        
            339.4780000        
             81.0101000        
             25.8780000        
              9.4522100        
              3.6656600        
              1.4674600        
              2.1562300        
              0.7489970        
              0.2831450        
              0.0998317        
      0.0018516              0.0000000              0.0000000              0.0000000        
      0.0142062              0.0000000              0.0000000              0.0000000        
      0.0699995              0.0000000              0.0000000              0.0000000        
      0.2400790              0.0000000              0.0000000              0.0000000        
      0.4847620              0.0000000              0.0000000              0.0000000        
      0.3352000              0.0000000              0.0000000              0.0000000        
      0.0000000             -0.00278217             0.0000000              0.0000000        
      0.0000000             -0.0360499              0.0000000              0.0000000        
      0.0000000             -0.1166310              0.0000000              0.0000000        
      0.0000000              0.0968328              0.0000000              0.0000000        
      0.0000000              0.6144180              0.0000000              0.0000000        
      0.0000000              0.4037980              0.0000000              0.0000000        
      0.0000000              0.0000000             -0.2529230              0.0000000        
      0.0000000              0.0000000              0.0328517              0.0000000        
      0.0000000              0.0000000              1.0812500              0.0000000        
      0.0000000              0.0000000              0.0000000              1.0000000        
* P-type functions
    10    3
            339.4780000        
             81.0101000        
             25.8780000        
              9.4522100        
              3.6656600        
              1.4674600        
              2.1562300        
              0.7489970        
              0.2831450        
              0.0998317        
      0.00456462             0.0000000              0.0000000        
      0.03369360             0.0000000              0.0000000        
      0.13975500             0.0000000              0.0000000        
      0.33936200             0.0000000              0.0000000        
      0.45092100             0.0000000              0.0000000        
      0.23858600             0.0000000              0.0000000        
      0.0000000             -0.01776530             0.0000000        
      0.0000000              0.27405800             0.0000000        
      0.0000000              0.78542100             0.0000000        
      0.0000000              0.0000000              1.00000000      

P2           0.00000000    0.00000000   -1.50670417 Angstrom
P3           0.00000000    1.28233653    0.32510185 Angstrom
P4          -1.11053601   -0.64116827    0.32510185 Angstrom
End of basis
>>> Do While <<<
&SEWARD
&RASSCF
  Symmetry =1
  spin= 2; Nactel= 13 0 0; Inactive= 17 8; Ras2= 8 5
  LUMORB
&SLAPAF
>>> EndDo <<<
All; AScii
要以感性的思想和认识,理性的思考和行动对待科研。
3楼2013-10-25 07:30:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 思雨G十年 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见