24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2438  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科211,293分请求调剂 +8 莲菜就是藕吧 2026-04-03 9/450 2026-04-05 19:12 by 蓝云思雨
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +4 一名26考研生 2026-04-04 4/200 2026-04-05 10:01 by barlinike
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +10 Naiko 2026-04-04 10/500 2026-04-05 09:40 by sam3303
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 11408 一志愿西电,277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 085501一志愿天工大,机械专硕求调剂,跨材料 +3 33上 2026-04-03 3/150 2026-04-03 14:08 by 1753564080
[考研] 282求调剂 +5 呼吸都是减肥 2026-03-31 5/250 2026-04-03 12:03 by 1753564080
[考研] 初试成绩337找调剂 +3 ??? ?. ? 2026-04-03 3/150 2026-04-03 11:43 by 土木硕士招生
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 一志愿山东大学,085600,344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[论文投稿] chinese chemical letters英文版投稿求助 120+4 Yishengeryi 2026-03-30 6/300 2026-04-02 17:19 by Yishengeryi
[考研] 297 地理学070500 复试求调剂 +3 小圆圈圈ooo 2026-03-30 3/150 2026-03-30 21:05 by 余震yz
信息提示
请填处理意见