24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2437  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂 +17 艾志恒 2026-04-03 18/900 2026-04-06 07:10 by jj987
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 285求调剂 +7 AZMK 2026-04-04 9/450 2026-04-06 00:06 by 永字号
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 086000生物与医药 初试274求调剂 +6 小叮当来了 2026-03-30 7/350 2026-04-05 20:30 by lys0704
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 一志愿郑州大学材料与化工085600,求调剂 +24 吃的不少 2026-04-02 24/1200 2026-04-04 23:20 by 永字号
[考研] 283求调剂 +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 0856材料与化工调剂,339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 296求调剂 +4 汪!?! 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 本科211总分289,08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
信息提示
请填处理意见