版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

思雨G十年

金虫 (正式写手)

应助: 32 (小学生)
金币: 20
散金: 2953
红花: 17
帖子: 897
在线: 328.9小时
虫号: 807115
注册: 2009-07-11
性别: GG
专业: 理论和计算化学

[求助] 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

回复此楼

» 猜你喜欢

海南师范大学2026年博士研究生招收 (在职想提升学历人员可报考) 申请考核制已经有0人回复
锂离子电池循环寿命衰减过快，求机理分析与改进思路已经有2人回复
物理化学论文润色/翻译怎么收费? 已经有113人回复
求博导收留已经有5人回复
九江学院2026年最新高层次人才招聘公告已经有0人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士等已经有20人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有19人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有1人回复
深圳信息职业技术学院2026年教师招聘公告（最新）已经有0人回复
固态锂金属电池迁移数计算求助已经有1人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

在Gaussian09计算中考虑标量相对论和旋轨耦合已经有9人回复
全国第二届理论与计算化学软件Molpro Molcas培训班已经有19人回复

要以感性的思想和认识，理性的思考和行动对待科研。

1楼 2013-09-28 18:04:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

思雨G十年

金虫 (正式写手)

应助: 32 (小学生)
金币: 20
散金: 2953
红花: 17
帖子: 897
在线: 328.9小时
虫号: 807115
注册: 2009-07-11
性别: GG
专业: 理论和计算化学

我计算得到的正确结果的输入：
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

赞一下(1人)

回复此楼

要以感性的思想和认识，理性的思考和行动对待科研。

3楼2013-10-25 07:30:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 3 个回答

gmy1990

荣誉版主 (著名写手)

QC强帖: 1
应助: 163 (高中生)
贵宾: 0.782
金币: 13422.7
散金: 526
红花: 42
沙发: 4
帖子: 1627
在线: 1019.7小时
虫号: 745920
注册: 2009-04-11
专业: 理论和计算化学
管辖: 量子化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
思雨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的回帖

查看全部 3 个回答