chj787

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 933
帖子: 104
在線: 42.2小時(shí)
蟲號: 798564
注冊: 2009-06-24
專業(yè): 原子和分子物理

[求助] molpro計(jì)算自旋軌道耦合輸入求助

剛學(xué)習(xí)molpro, 試做一個(gè)F原子的自旋軌道耦合計(jì)算，輸入文件如下：
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}
該輸入文件能夠順利的執(zhí)行完，輸出結(jié)果如下：

*** 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，而且能量一樣，好像沒有進(jìn)行SO運(yùn)算，求高手指點(diǎn)，怎么解決，輸入文件怎么寫？剛開始研究molpro, 急需幫助，謝謝指點(diǎn)�。�

[ Last edited by chj787 on 2012-2-27 at 20:56 ]

回復(fù)此樓

» 猜你喜歡

理論與計(jì)算化學(xué)方向調(diào)劑——初試388分（數(shù)學(xué)144）已經(jīng)有0人回復(fù)
上海工程技術(shù)大學(xué)生態(tài)與功能紡織品開發(fā)研究團(tuán)隊(duì)招收學(xué)碩專碩研究生已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有95人回復(fù)
青島大學(xué)化學(xué)化工學(xué)院2026年招收調(diào)劑碩士研究生已經(jīng)有0人回復(fù)
新能源材料方向考研調(diào)劑招生（長江師范學(xué)院）已經(jīng)有0人回復(fù)
0854電子信息272求調(diào)劑！已經(jīng)有0人回復(fù)
化學(xué)-化工-材料-食品-生物-環(huán)境2026年研究生調(diào)劑（可以不考數(shù)學(xué)）已經(jīng)有0人回復(fù)
化學(xué)-化工-材料-食品-生物-環(huán)境-農(nóng)業(yè) 2026年研究生調(diào)劑（可以不考數(shù)學(xué)）已經(jīng)有0人回復(fù)
070300一志愿211，312分求調(diào)劑院校已經(jīng)有4人回復(fù)
26年儲能、電池方向博士申請已經(jīng)有0人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

全國第二屆理論與計(jì)算化學(xué)軟件Molpro Molcas培訓(xùn)班已經(jīng)有19人回復(fù)
求助 molpro運(yùn)行出錯(cuò)！急已經(jīng)有5人回復(fù)
關(guān)于Molpro自旋軌道耦合輸出的問題已經(jīng)有9人回復(fù)
【討論】molpro, molcas和gaussian的優(yōu)缺點(diǎn) 已經(jīng)有3人回復(fù)
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧已經(jīng)有105人回復(fù)
【求助】請教molpro計(jì)算中出現(xiàn)的一個(gè)問題已經(jīng)有7人回復(fù)
【求助】請問Gaussian計(jì)算帶電荷有機(jī)物時(shí)電荷和自旋多重度如何選擇已經(jīng)有18人回復(fù)
【求助】請教使用orca計(jì)算出的這個(gè)自旋密度圖是怎么畫的？已經(jīng)有5人回復(fù)
【求助】VASP考慮自旋軌道耦合的話，如何編譯已經(jīng)有17人回復(fù)
【求助】一個(gè)二價(jià)鈷配合物的自旋多重度--附不同自旋多重度的計(jì)算結(jié)果已經(jīng)有11人回復(fù)
【求助】Molpro如何執(zhí)行旋軌耦合計(jì)算如何執(zhí)行態(tài)平均計(jì)算? 已經(jīng)有10人回復(fù)
【求助】怎樣用高斯計(jì)算自旋密度已經(jīng)有15人回復(fù)
【求助】分子篩量化計(jì)算求助-關(guān)于自旋污染問題已經(jīng)有8人回復(fù)
【求助】td計(jì)算三重態(tài)問題已經(jīng)有8人回復(fù)

1樓 2012-02-27 20:13:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

beefly

專家顧問 (職業(yè)作家)

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
QC強(qiáng)帖: 12
應(yīng)助: 88 (初中生)
金幣: 13022.2
散金: 27307
紅花: 300
帖子: 4673
在線: 927.1小時(shí)
蟲號: 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 計(jì)算模擬

小紅豆(金幣+5, 專家考核): 謝謝 2012-02-27 21:19:50

因?yàn)閏i只計(jì)算了一個(gè)wf,9,5,1分量

贊一下(1人)

回復(fù)此樓

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

2樓2012-02-27 21:02:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

chj787

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 933
帖子: 104
在線: 42.2小時(shí)
蟲號: 798564
注冊: 2009-06-24
專業(yè): 原子和分子物理

引用回帖:

樓: Originally posted by beefly at 2012-02-27 21:02:00:
因?yàn)閏i只計(jì)算了一個(gè)wf,9,5,1分量

謝謝這位老師，弱弱的問下，molpro不能對單一組態(tài)進(jìn)行自旋耦合計(jì)算嗎？計(jì)算LS耦合，是不是只能進(jìn)行多參考組態(tài)CI計(jì)算？

贊一下

回復(fù)此樓

3樓2012-02-28 09:36:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

beefly

專家顧問 (職業(yè)作家)

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
QC強(qiáng)帖: 12
應(yīng)助: 88 (初中生)
金幣: 13022.2
散金: 27307
紅花: 300
帖子: 4673
在線: 927.1小時(shí)
蟲號: 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 計(jì)算模擬

2P態(tài)的旋軌耦合矩陣本來是3*3，你簡化成1*1矩陣，肯定不行，除非三個(gè)分量之間的相互作用都為0

贊一下(1人)

回復(fù)此樓

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

4樓2012-02-28 10:20:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑 +6	青春裁為三截 2026-03-29	6/300	2026-03-30 12:48 by 159357hjz
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 297求調(diào)劑 +13	田洪有 2026-03-26	13/650	2026-03-30 11:33 by 替代品000
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	7/350	2026-03-30 11:02 by 1172367218
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6	來財(cái)aa 2026-03-25	6/300	2026-03-30 09:48 by 青海小西牛
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 330分求調(diào)劑 +5	qzenlc 2026-03-29	5/250	2026-03-29 07:37 by 無際的草原
[考研] 學(xué)碩274求調(diào)劑 +9	Li李魚 2026-03-26	9/450	2026-03-28 21:42 by bymhappy
[考研] 085701求調(diào)劑初試286分 +4	secret0328 2026-03-28	4/200	2026-03-28 21:09 by 15366876211
[考研] 322求調(diào)劑 +5	舊吢 2026-03-24	5/250	2026-03-28 13:26 by Iveryant
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3	熠然好運(yùn)氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 343求調(diào)劑 +4	贈我一本書 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin