24小时热门版块排行榜    

查看: 1888  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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个分子上都加上?
求助!!
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

qizhuixia

新虫 (初入文坛)

研究的是原子间库伦decay,其中一个水分子的氧原子处于激发态,退激后把相邻水分子中的氧原子激发,计算这个氧原子的衰败宽度
2楼2019-04-08 17:22:15
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +3 Tide man 2026-08-14 4/200 2026-08-14 16:33 by Yeuchan
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +5 启萌科技 2026-08-12 16/800 2026-08-14 16:02 by 启萌科技
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] filecode +5 cratir 2026-08-14 9/450 2026-08-14 12:25 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见