24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2436  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-03-30 10/500 2026-04-05 17:57 by jj987
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 359求调剂22408 +3 123456789qw 2026-03-31 3/150 2026-04-05 10:09 by zhq0425
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 293分求调剂,外语为俄语 +6 加一一九 2026-03-31 6/300 2026-04-04 14:57 by 聪明的大松鼠
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 266分,求材料相关专业调剂 +13 哇呼哼呼哼 2026-03-30 15/750 2026-04-03 15:24 by arrow8852
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
信息提示
请填处理意见