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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » 求助大佬！！！关于CI计算如何输入H2O的基组
81/1返回列表
查看: 1610  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

qizhuixia

新虫 (初入文坛)

[求助] 求助大佬！！！关于CI计算如何输入H2O的基组 已有1人参与

情况是这样的，我想计算两个水分子（通过氢键相连）这个系统的衰败宽度，想用cc-pvdz这个基组，这样输入是否正确（只是data部分）：

$DATA
water ccpvdz
CN 1

O 8.0 0.61613258971905416 1.2945922650017754 0.43576112323695337
S   8
  1  11720.0000000              0.0007100
  2   1759.0000000              0.0054700
  3    400.8000000              0.0278370
  4    113.7000000              0.1048000
  5     37.0300000              0.2830620
  6     13.2700000              0.4487190
  7      5.0250000              0.2709520
  8      1.0130000              0.0154580
S   8
  1  11720.0000000             -0.0001600
  2   1759.0000000             -0.0012630
  3    400.8000000             -0.0062670
  4    113.7000000             -0.0257160
  5     37.0300000             -0.0709240
  6     13.2700000             -0.1654110
  7      5.0250000             -0.1169550
  8      1.0130000              0.5573680
S   1
  1      0.3023000              1.0000000
P   3
  1     17.7000000              0.0430180
  2      3.8540000              0.2289130
  3      1.0460000              0.5087280
P   1
  1      0.2753000              1.0000000
D   1
  1      1.1850000              1.0000000

H 1.0 1.20979641732937180 1.1399408154279724 1.2384371269171006
S   3
  1     13.0100000              0.0196850
  2      1.9620000              0.1379770
  3      0.4446000              0.4781480
S   1
  1      0.1220000              1.0000000
P   1
  1      0.7270000              1.0000000

H 1.0 0.14980528497444778 0.39660064520917754 0.20962755813354148
S   3
  1     13.0100000              0.0196850
  2      1.9620000              0.1379770
  3      0.4446000              0.4781480
S   1
  1      0.1220000              1.0000000
P   1
  1      0.7270000              1.0000000

O 8.0 -0.57771364708940753 -1.2209515091573029 -0.42399613287644144
S   8
  1  11720.0000000              0.0007100
  2   1759.0000000              0.0054700
  3    400.8000000              0.0278370
  4    113.7000000              0.1048000
  5     37.0300000              0.2830620
  6     13.2700000              0.4487190
  7      5.0250000              0.2709520
  8      1.0130000              0.0154580
S   8
  1  11720.0000000             -0.0001600
  2   1759.0000000             -0.0012630
  3    400.8000000             -0.0062670
  4    113.7000000             -0.0257160
  5     37.0300000             -0.0709240
  6     13.2700000             -0.1654110
  7      5.0250000             -0.1169550
  8      1.0130000              0.5573680
S   1
  1      0.3023000              1.0000000
P   3
  1     17.7000000              0.0430180
  2      3.8540000              0.2289130
  3      1.0460000              0.5087280
P   1
  1      0.2753000              1.0000000
D   1
  1      1.1850000              1.0000000

H 1.0 -1.5746335435163650 -1.2212884001307356 -0.39367330877139639
S   3
  1     13.0100000              0.0196850
  2      1.9620000              0.1379770
  3      0.4446000              0.4781480
S   1
  1      0.1220000              1.0000000
P   1
  1      0.7270000              1.0000000

H 1.0 -0.39966995256225191 -1.4935067198219192 -1.2426390564127234
S   3
  1     13.0100000              0.0196850
  2      1.9620000              0.1379770
  3      0.4446000              0.4781480
S   1
  1      0.1220000              1.0000000
P   1
  1      0.7270000              1.0000000

$END

如果想用cc-pvdz基组加上1s1p1d基组让计算更精确，我是应该只在两个氧分子加，还是在6个分子上都加上？
求助！！
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2019-04-08 16:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)
研究的是原子间库伦decay，其中一个水分子的氧原子处于激发态，退激后把相邻水分子中的氧原子激发，计算这个氧原子的衰败宽度
赞一下 回复此楼
2楼2019-04-08 17:22:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+3 2019-04-09 09:03:22
qizhuixia: 金币+40, ★★★★★最佳答案 2019-04-09 16:14:12
这样写是可以的。如果直接加一些基函数进去 也就是相当于给基组加弥散。这个确实算氢键这种弱相互作用是必要的。  要加就给所有原子都加。  不过与其这样不如直接用 aug-cc-pvdz 基组  也就是 对H 加一个 s 一个p 对而O加 一个s 一个p 一个d。
赞一下 回复此楼
3楼2019-04-08 22:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by paramecium86 at 2019-04-08 22:28:22
这样写是可以的。如果直接加一些基函数进去 也就是相当于给基组加弥散。这个确实算氢键这种弱相互作用是必要的。  要加就给所有原子都加。  不过与其这样不如直接用 aug-cc-pvdz 基组  也就是 对H 加一个 s 一个p 对 ...

如果都加的话到7s7p7d，计算就慢的不行。。。如果只给氧原子加合理么？
赞一下 回复此楼
4楼2019-04-08 23:14:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by paramecium86 at 2019-04-08 22:28:22
这样写是可以的。如果直接加一些基函数进去 也就是相当于给基组加弥散。这个确实算氢键这种弱相互作用是必要的。  要加就给所有原子都加。  不过与其这样不如直接用 aug-cc-pvdz 基组  也就是 对H 加一个 s 一个p 对 ...

还有一个问题，如果想要结果更准确，用cc-pvtz基组，可不可以只给氧原子换成tz，其他氢原子还是dz，这样行么……
赞一下 回复此楼
5楼2019-04-08 23:54:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
★ ★ ★ ★ ★
zhou2009: 金币+5 2019-04-09 09:03:53
引用回帖:
4楼: Originally posted by qizhuixia at 2019-04-08 10:14:24
如果都加的话到7s7p7d，计算就慢的不行。。。如果只给氧原子加合理么？...

倒是也不是不行 类似于 jul-cc-pVnZ 系列基组，  就是不对H加弥散 但是保持对其它原子加弥散的一系列基组。 这样做说得通。不过 我不确定 单纯加弥散加到7s 7p 7d 是不是合理   这也太不好收敛了。而且 单纯的加spd 而不升高基组最高的角动量是不是对结果有很大帮助就不好说了。 差不多大小的基组里，这算出来 感觉结果不一定有
aug-cc-pVQZ 这种含有高角动量的基组的结果好。
赞一下 回复此楼
6楼2019-04-08 23:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
★ ★ ★ ★ ★
zhou2009: 金币+5 2019-04-09 09:04:23
引用回帖:
5楼: Originally posted by qizhuixia at 2019-04-08 10:54:14
还有一个问题，如果想要结果更准确，用cc-pvtz基组，可不可以只给氧原子换成tz，其他氢原子还是dz，这样行么……...

分子不大的话 真的很难说服别人 (比如以后发文章的审稿人 )这样做的原因。。。。而且这分子里 H也并不是无关紧要的部分。并且前提是之前已经用了超大的弥散的dz基组算的。  哪怕是把弥散函数减下来 O H 都用 ang-cc-pvtz 也感觉上比 俩氧原子把弥散加到7s 7p 7d,7f的tz基组 而一边6个H用小小的dz基组更加合理。
赞一下 回复此楼
7楼2019-04-09 00:05:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qizhuixia

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by paramecium86 at 2019-04-09 00:05:22
分子不大的话 真的很难说服别人 (比如以后发文章的审稿人 )这样做的原因。。。。而且这分子里 H也并不是无关紧要的部分。并且前提是之前已经用了超大的弥散的dz基组算的。  哪怕是把弥散函数减下来 O H 都用 ang-c ...

谢谢老哥，我再想想
赞一下 回复此楼
8楼2019-04-09 16:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qizhuixia 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 289材料与化工（085600）B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
信息提示
关闭
请填处理意见
关闭