24小时热门版块排行榜    

查看: 2622  |  回复: 3

chj787

金虫 (小有名气)

[求助] molpro计算自旋轨道耦合输入求助

刚学习molpro, 试做一个F原子的自旋轨道耦合计算,输入文件如下:
geometry={angstrom;
  F
  }
   {hf
   occ,2,1,1,,1
   wf,9,5,1
   open,1.5}
  
    {casscf
    occ,2,1,1,,1
    wf,9,2,1
    wf,9,3,1
    wf,9,5,1
    }
  
  {ci;core;occ,2,1,1,,1;wf,9,5,1;save,6000.2}
   
   text,casscf,occ,9,5,1
   {ci;hlsmat,ecp,6000.2}
该输入文件能够顺利的执行完,输出结果如下:

*** casscf,occ,9,5,1


Reading current basis input from record  610.1
1PROGRAM * CI (Multireference internally contracted CI)     Authors: H.-J. Werner, P.J. Knowles, 1987


   ******************************
   *** Spin-orbit calculation ***
   ******************************


   Spin-orbit matrix elements
   ==========================

Wavefunction restored from record  6000.2  Symmetry=5  S=0.5  NSTATE=1

No symmetry adaption


Spin-orbit eigenstates   (energies)
======================

   Nr         E             E-E0         E-E0           E-E(1)      E-E(1)      E-E(1)
            (au)            (au)        (cm-1)           (au)       (cm-1)        (eV)
   1   -99.52523070     0.00000000        0.00      0.00000000        0.00      0.0000
   2   -99.52523070     0.00000000        0.00      0.00000000        0.00      0.0000


Property matrices transformed in SO basis (not sym. adapted)
============================================================
   No matrix element larger than 1E-10
   No matrix element larger than 1E-10
   No matrix element larger than 1E-10
   No matrix element larger than 1E-10
   No matrix element larger than 1E-10
   No matrix element larger than 1E-10


**********************************************************************************************************************************
DATASETS  * FILE   NREC   LENGTH (MB)   RECORD NAMES
              1      22        2.21       500      610      700      900      950      970     1000      129      960     1100   
                                          VAR    BASINP    GEOM    SYMINP    ZMAT    AOBASIS   BASIS     P2S    ABASIS      S
                                         1400     1410     1200     1210     1080     1600     1650     1300     1700      960(1)
                                           T        V       H0       H01     AOSYM     SMH    MOLCAS    ERIS     OPER    ABASIS   
                                         1380     1700(1)
                                         JKOP     OPER   

              2       5        0.30       700     1000     2100     2140     6000   
                                         GEOM     BASIS     RHF     MCSCF    MRCI   

PROGRAMS   *        TOTAL        CI        CI    CASSCF        HF       INT
CPU TIMES  *         0.50      0.06      0.07      0.01      0.01      0.01
REAL TIME  *         0.54 SEC
DISK USED  *        11.06 MB      
**********************************************************************************************************************************

        CI              CASSCF          HF-SCF  
    -99.52523070    -99.37108022    -99.37186194
**********************************************************************************************************************************
Variable memory released

请问高手为什么输出中有 No symmetry adaption,而且能量一样,好像没有进行SO运算,求高手指点,怎么解决,输入文件怎么写?刚开始研究molpro, 急需帮助,谢谢指点啊!

[ Last edited by chj787 on 2012-2-27 at 20:56 ]
回复此楼

» 猜你喜欢

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

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

小红豆(金币+5, 专家考核): 谢谢 2012-02-27 21:19:50
因为ci只计算了一个wf,9,5,1分量
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-02-27 21:02:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chj787

金虫 (小有名气)

引用回帖:
: Originally posted by beefly at 2012-02-27 21:02:00:
因为ci只计算了一个wf,9,5,1分量

谢谢这位老师,弱弱的问下,molpro不能对单一组态进行自旋耦合计算吗?计算LS耦合,是不是只能进行多参考组态CI计算?
3楼2012-02-28 09:36:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

2P态的旋轨耦合矩阵本来是3*3,你简化成1*1矩阵,肯定不行,除非三个分量之间的相互作用都为0
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2012-02-28 10:20:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chj787 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见