24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2429  |  回复: 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的回帖

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的回帖
查看全部 4 个回答

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的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

2P态的旋轨耦合矩阵本来是3*3,你简化成1*1矩阵,肯定不行,除非三个分量之间的相互作用都为0
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2012-02-28 10:20:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602化工求调剂(331分) +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 08工科求调剂286 +4 tgs_001 2026-03-28 4/200 2026-03-30 21:03 by dophin1985
[考研] 085600 295分求调剂 +7 W55j 2026-03-30 9/450 2026-03-30 20:36 by dick_runner
[考研] 08工科,295,接受跨专业调剂 +3 lmnlzy 2026-03-30 3/150 2026-03-30 17:49 by wangjy2002
[考研] 291求调剂 +8 HanBeiNingZC 2026-03-24 8/400 2026-03-30 17:01 by 无际的草原
[考研] 298求调剂 +3 什么是胖头鱼 2026-03-30 5/250 2026-03-30 14:41 by 青海小西牛
[考研] 南京大学化学调剂 +10 景随风 2026-03-29 15/750 2026-03-30 11:21 by limeifeng
[考研] 材料与化工272求调剂 +21 阿斯蒂芬2004 2026-03-28 21/1050 2026-03-30 10:52 by 晴空210210
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历 +8 cccchenso 2026-03-29 8/400 2026-03-29 23:53 by 我是小康
[考研] 085600,材料与化工321分求调剂 +10 大馋小子 2026-03-28 10/500 2026-03-29 23:35 by 飞行日记西
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +7 AZMK 2026-03-28 11/550 2026-03-28 21:01 by xxxsssccc
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
[考研] 考研调剂 +10 呼呼?~+123456 2026-03-24 10/500 2026-03-27 11:46 by wangjy2002
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
信息提示
请填处理意见