应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » GAMESS的inp里面$VEC这块是怎么生成的？从GAMESS的输出文件里.out里读取不到么？
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 2358  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dkstatsra

新虫 (初入文坛)

[交流] GAMESS的inp里面$VEC这块是怎么生成的？从GAMESS的输出文件里.out里读取不到么？ 已有2人参与

下面是exam19的原文：
! EXAM 19.
! Spin-orbit coupling example.
! This run duplicates the results shown in Table 3 of
! T.R.Furlani, H.F.King, J.Chem.Phys. 82, 5577-83(1985),
! GAMESS 1e-= 114.3851, 2e-= -49.4168, lit=114.38,-49.42
!
! Energies for the singlet CI are
!   State=   1 Energy =   -54.868531216 (1-delta)
!   State=   2 Energy =   -54.868531216 (1-delta)
!   State=   3 Energy =   -54.798836732 (1-sigma-plus)
! Energies for the triplet CI are
!   State=   1 Energy =   -54.938225701 (3-sigma-minus)
! Final energies of all 6 levels in the pi**2 configuration,
! after diagonalization of the spin-orbit Hamiltonian, are
!   BREIT RELATIVE E= -15296.570, -15296.432, -15296.432,
!   BREIT RELATIVE E= 0.0, 0.0, and +15296.570 wavenumbers.
! If run as OPERAT=HSO1, with ZEFF taken as true atomic Z,
! then inclusion of only the 1e- operator is 114.3851, and
!   ZEFF RELATIVE E= -15296.859, -15296.432, -15296.432,
!   ZEFF RELATIVE E= 0.0, 0.0, and +15296.859 wavenumbers.
!
! Why are there six levels? The singlet-delta is two roots,
! the singlet-sigma-plus is a third.  During the CI, the
! spatial triplet-sigma-minus is one CSF, with alpha/alpha
! spin, hence IROOTS=3,1.  The final spin-orbit Hamiltonian
! includes all three triplet spin states, namely adding the
! ab+ba and beta/beta triplets.  So, 2+1+3=6 levels.  You
! can work out for yourself these have the quantum number
! omega=0,0,1,2.  Only the omega=0 states can interact,
! raising the triplet's degeneracy and slightly affecting
! the singlet-sigma-plus state's position.
!
! Note that the lower multiplicity DRT1 is done in C1
! symmetry to generate both components of the delta state.
!
$CONTRL SCFTYP=NONE MULT=3 CITYP=GUGA RUNTYP=TRANSITN
         UNITS=BOHR $END
$SYSTEM TIMLIM=1 $END
$BASIS  GBASIS=N31 NGAUSS=6 $END
$TRANST OPERAT=HSO2 NFZC=3 NOCC=5 NUMVEC=1 NUMCI=2
         IROOTS(1)=3,1 $END
$DRT1   GROUP=C1  IEXCIT=2 NFZC=3 NDOC=1 NVAL=1 $END
$DRT2   GROUP=C4V IEXCIT=2 NFZC=3 NALP=2 $END
$DATA
Imidogen radical
Cnv 4

Nitrogen    7.0
Hydrogen    1.0     0.0   0.0   1.9748
$END

--- ROHF ORBITALS --- GENERATED AT 12:04:18 29 MAR 90 ( 88)
IMIDOGEN RADICAL
E(ROHF)= -54.9382257007, E(NUC)=    3.5446627507,    8 ITERS
$VEC1
1  1 9.97073281E-01 1.91214515E-02 0.00000000E+00 0.00000000E+00 2.36016032E-03
1  2-9.66111886E-03 0.00000000E+00 0.00000000E+00-2.13067963E-03-5.62009080E-04
1  3 3.44278483E-03
2  1-1.79300889E-01 4.10265484E-01 0.00000000E+00 0.00000000E+00 2.27571384E-01
2  2 3.83219392E-01 0.00000000E+00 0.00000000E+00 1.18321655E-01 2.61608099E-01
2  3 7.12373427E-02
3  1-1.26216963E-01 2.78613459E-01 0.00000000E+00 0.00000000E+00-4.52995544E-01
3  2 4.82312187E-01 0.00000000E+00 0.00000000E+00-2.72710767E-01-2.31321603E-01
3  3-2.10846345E-01
4  1 0.00000000E+00 0.00000000E+00-6.71321829E-01 0.00000000E+00 0.00000000E+00
4  2 0.00000000E+00-4.70445403E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00
4  3 0.00000000E+00
5  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.71321829E-01 0.00000000E+00
5  2 0.00000000E+00 0.00000000E+00 4.70445403E-01 0.00000000E+00 0.00000000E+00
5  3 0.00000000E+00
6  1 7.36851433E-02-1.43290334E-01 0.00000000E+00 0.00000000E+00-3.26478078E-01
6  2-1.02184944E+00 0.00000000E+00 0.00000000E+00-8.07276620E-01 7.82959641E-02
6  3 1.76226457E+00
7  1 0.00000000E+00 0.00000000E+00-9.58419043E-01 0.00000000E+00 0.00000000E+00
7  2 0.00000000E+00 1.07141084E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
7  3 0.00000000E+00
8  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 9.58419043E-01 0.00000000E+00
8  2 0.00000000E+00 0.00000000E+00-1.07141084E+00 0.00000000E+00 0.00000000E+00
8  3 0.00000000E+00
9  1-6.00738465E-02 5.99272598E-01 0.00000000E+00 0.00000000E+00-7.39008372E-01
9  2-3.36089417E-01 0.00000000E+00 0.00000000E+00 1.18000819E+00-5.28801067E-01
9  3-3.52885757E-02
10  1 1.39256694E-04 6.71797254E-01 0.00000000E+00 0.00000000E+00-3.93535986E-01
10  2-9.13508814E-01 0.00000000E+00 0.00000000E+00-3.04066289E-01 1.21292686E+00
10  3-4.07553827E-01
11  1 5.08094255E-02-1.41498587E+00 0.00000000E+00 0.00000000E+00-6.27815314E-01
11  2 2.22349514E+00 0.00000000E+00 0.00000000E+00 1.15113704E+00 3.26442418E-01
11  3-1.58015531E+00
$END
$VEC这部分是怎么generate出来的？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-20 11:05:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bittersmile

新虫 (初入文坛)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+1 2014-06-23 05:59:11
$vec在那个PUNCH或者.dat文件里。
赞一下 回复此楼
2楼2014-06-20 22:27:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★
gmy1990: 金币+1 2014-06-23 05:59:18
$VEC之前的这段信息已经说了生成方法：用ROHF做的。
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2014-06-21 01:09:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by beefly at 2014-06-21 01:09:12
$VEC之前的这段信息已经说了生成方法：用ROHF做的。

好吧，我可能没说明白，我想知道在哪里找到这段信息。$VEC的
赞一下 回复此楼
4楼2014-06-24 17:08:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by bittersmile at 2014-06-20 22:27:11
$vec在那个PUNCH或者.dat文件里。

$SCF NPUNCH=2 ？还是说怎么设置？这个是在.out的尾部输出MO吗？
赞一下 回复此楼
5楼2014-06-24 17:17:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dkstatsra

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by bittersmile at 2014-06-20 22:27:11
$vec在那个PUNCH或者.dat文件里。

还有，这个 .dat文件在哪里出来的？是和 .out在一个文件夹里？我怎么没有这个
赞一下 回复此楼
6楼2014-06-24 17:19:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dkstatsra 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学0703-一志愿211-338分求调剂 +8 vants 2026-04-05 8/400 2026-04-06 06:17 by houyaoxu
[考研] 0703调剂，一志愿天津大学319分 +8 haaaabcd 2026-04-05 8/400 2026-04-06 06:16 by houyaoxu
[考研] 081700，311，求调剂 +14 冬十三 2026-04-04 15/750 2026-04-06 01:57 by BruceLiu320
[考研] 材料与化工371求调剂 +13 陪琳看海 2026-04-04 14/700 2026-04-06 00:58 by fmesaito
[考研] 332求调剂 +16 小小孟... 2026-04-05 17/850 2026-04-06 00:52 by fmesaito
[考研] 生物学308分求调剂（一志愿华东师大） +8 相信必会光芒万� 2026-04-05 10/500 2026-04-05 12:19 by Hdyxbekcb
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-03-31 5/250 2026-04-04 19:45 by 1753564080
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 生物学308分求调剂（一志愿华东师大） +7 相信必会光芒万� 2026-04-02 7/350 2026-04-03 16:48 by rzh123456
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[考研] 085602化工求调剂（331分） +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 302求调剂一志愿北航070300，本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 285求调剂 +7 AZMK 2026-03-30 13/650 2026-04-01 17:00 by 七度不信任
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
信息提示
关闭
请填处理意见
关闭