计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

31

帖子

0

威望

187

eV
积分
218

Level 3 能力者

16231#
发表于 Post on 2022-8-23 11:32:32 | 只看该作者 Only view this author

请教GaussView画图

老师们好,阅读文献时看到下面这张图,请问此图通过GaussView可以画出来吗?还是需要别的软件进行辅助?

1016

帖子

0

威望

2126

eV
积分
3142

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

16232#
发表于 Post on 2022-8-23 12:03:23 | 只看该作者 Only view this author
可以。GaussView可以在原子上标元素名,但是那个R(H-Ring)和中间的虚线要用图像处理软件弄一下。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16233#
发表于 Post on 2022-8-23 16:26:14 | 只看该作者 Only view this author
chands 发表于 2022-8-23 10:26
假如你用结合自由能衡量配位能力,是要算(Fe3+) + L = (Fe3+)L的结合自由能(我这里只算气相一个配体的情 ...

好的,谢谢老师,我准备先用低基组优化再用高基组算单点能,也就是说我算不同元素的时候基组不统一得出来的配位能之间是可以比较的吗?

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16234#
发表于 Post on 2022-8-23 16:28:55 | 只看该作者 Only view this author
sobereva 发表于 2022-8-23 10:44
不同元素不需要基组统一
def、def2对Fe都不是赝势基组。直接用def2全电子基组就完了,如果想节约时间考 ...

好的,谢谢老师,老师也就是说不同元素基组不统一算出来的配位结合能是可以比较的是吗?

1016

帖子

0

威望

2126

eV
积分
3142

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

16235#
发表于 Post on 2022-8-23 16:45:19 | 只看该作者 Only view this author
KCM@SHU 发表于 2022-8-23 16:26
好的,谢谢老师,我准备先用低基组优化再用高基组算单点能,也就是说我算不同元素的时候基组不统一得出来 ...

当然可以,否则混合基组就不能用了。不过注意算某个反应所有反应物和产物的级别要一致。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

16236#
发表于 Post on 2022-8-23 17:34:20 | 只看该作者 Only view this author
KCM@SHU 发表于 2022-8-23 09:26
好的,谢谢老师,我准备先用低基组优化再用高基组算单点能,也就是说我算不同元素的时候基组不统一得出来 ...

可以,但大小最好不要差太多,举个极端的例子,比如一个元素用6-31G*,另一个用aug-cc-pVQZ,那估计不行。如果不确定差多少是多,那就统一把不同元素的基组都取得比较大,因为基组越大,结果对基组大小越不敏感。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16237#
发表于 Post on 2022-8-23 17:54:30 | 只看该作者 Only view this author
chands 发表于 2022-8-23 16:45
当然可以,否则混合基组就不能用了。不过注意算某个反应所有反应物和产物的级别要一致。

好的,谢谢老师

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16238#
发表于 Post on 2022-8-23 17:54:52 | 只看该作者 Only view this author
wzkchem5 发表于 2022-8-23 17:34
可以,但大小最好不要差太多,举个极端的例子,比如一个元素用6-31G*,另一个用aug-cc-pVQZ,那估计不行 ...

好的,谢谢老师

145

帖子

1

威望

1766

eV
积分
1931

Level 5 (御坂)

16239#
发表于 Post on 2022-8-23 18:49:13 | 只看该作者 Only view this author
请问,我用6-311g**没有找到自旋极化单重态的稳定波函数,用def2-tzvpp能够找到,将这个找到的波函数读取用6-311g**去算程序报错,要怎么解决呀?

Y(W~VMXLL174@EAW~Q9NF@R.png (107.17 KB, 下载次数 Times of downloads: 9)

Y(W~VMXLL174@EAW~Q9NF@R.png

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16240#
发表于 Post on 2022-8-23 21:14:00 | 只看该作者 Only view this author
各位老师我想问一下做计算的时候计算机的那些配置参数会影响计算时间,和对应计算机配置条件下,各个泛函和基组的耗时大概是多少呢?感觉每次选择基组都没有概念不能选择最适合自己条件的计算方式?还想问一下大家在正式发文章的时候选用的泛函和基组进行结构优化和单点能的时候一般都用了多久呢?

1016

帖子

0

威望

2126

eV
积分
3142

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

16241#
发表于 Post on 2022-8-23 21:28:57 | 只看该作者 Only view this author
KCM@SHU 发表于 2022-8-23 21:14
各位老师我想问一下做计算的时候计算机的那些配置参数会影响计算时间,和对应计算机配置条件下,各个泛函和 ...

首先得说多大的体系,有什么原子。这个入门贴最后提到的博文建议看一看。
谈谈学量子化学如何正确地入门
http://sobereva.com/355

硬盘速度与内存容量对量子化学计算速度影响的测试
http://sobereva.com/397
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

73

帖子

0

威望

551

eV
积分
624

Level 4 (黑子)

16242#
发表于 Post on 2022-8-23 22:47:52 | 只看该作者 Only view this author
请问DM21泛函或者pySCF可以用gpu来计算吗?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

16243#
 楼主 Author| 发表于 Post on 2022-8-24 08:55:43 | 只看该作者 Only view this author
KCM@SHU 发表于 2022-8-23 21:14
各位老师我想问一下做计算的时候计算机的那些配置参数会影响计算时间,和对应计算机配置条件下,各个泛函和 ...

不存在“一般都用了多久”的说法
不同体系不同问题的计算情况千差万别
顶多说一般情况下,一个任务可接受的耗时撑死了也就两三天,还算不完的话应当考虑以其它方式计算或者用更好的计算条件

谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

16244#
 楼主 Author| 发表于 Post on 2022-8-24 08:56:08 | 只看该作者 Only view this author
WOOOOWOOOO 发表于 2022-8-23 18:49
请问,我用6-311g**没有找到自旋极化单重态的稳定波函数,用def2-tzvpp能够找到,将这个找到的波函数读取用 ...

截图里没体现报错
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

201

帖子

0

威望

892

eV
积分
1093

Level 4 (黑子)

16245#
发表于 Post on 2022-8-24 12:21:42 | 只看该作者 Only view this author
各位老师我想问一下,我在算K离子和水配位的结构优化的时候,我用显式水配位六个个后再加隐式的水算结构优化,最后报错了,我想问一下我在用gaussian view画构图的时候k离子与水上氧配位成键的键长选择会不会对计算有影响,如果不知道大概键长一般怎么设置呢?还有我在计算的时候我把报错的文件的out文件发现k离子与周围的氧原子的连接键只有五个了,但是别人的文章里面说k离子与水配位式6配位,我用gaussian view 构建的时候也是在K离子周围连了6个水的氧,这是我构型构建问题吗?下面是我的输入文件和报错,我还想问一下报错里面那个ERROR IS A HARDY PLANT... IT FLOURISHETH IN EVERY SOIL.
            MARTIN FARQUHAR TUPPER是什么意思呢?麻烦老师们有空指导一下
%chk=C:\Users\sxs\Desktop\MCK\k\k_h20.chk
# opt freq pbe1pbe/gen scrf=(solvent=water) geom=connectivity

Title Card Required

1 1
K                  0.90717301    0.18987341    0.00000000
O                  1.54008441    1.45569618    0.00000000
H                  2.50008441    1.45569618    0.00000000
H                  1.21962982    2.36063201    0.00000000
O                  0.16129686   -1.10803051    0.23271563
H                 -0.44728397   -0.41412671    0.49677412
H                 -0.31991579   -1.79658513   -0.23196374
O                  2.22700282   -0.54577486   -0.25001721
H                  3.00311366    0.01202218   -0.15993707
H                  2.45627425   -1.33658536   -0.74362978
O                 -0.28600726    0.77469466    0.40833885
H                 -0.35616633    1.69465912    0.14311612
H                 -0.98415661    0.56737050    1.03380202
O                  1.45968445    0.05188979    1.44748623
H                  0.98168249   -0.09921212    2.26619396
H                  2.39415321   -0.11913913    1.58575132
O                  0.44075773    0.35464907   -1.46054240
H                  0.91007220   -0.23311854   -2.05709345
H                 -0.03387615    1.01579723   -1.96966498

1 2 1.0 11 1.0 17 1.0 8 1.0 5 1.0 14 1.0
2 3 1.0 4 1.0
3
4
5 6 1.0 7 1.0
6
7
8 9 1.0 10 1.0
9
10
11 12 1.0 13 1.0
12
13
14 15 1.0 16 1.0
15
16
17 18 1.0 19 1.0
18
19

K    0
S  6  1.00
  31478.746764      .39838653994E-02
  4726.8876066      .30501759762E-01
  1075.4345353      .15073752622
  303.39811023      .51912939801
  98.327112831      1.0366957005
  33.636222177      .76398963199
S  3  1.00
  65.639209962     -.28242617106
  7.3162592218      1.6914935860
  2.8902580135      1.2965331953
S  3  1.00
  4.5459748965     -.76343555273E-02
  .70404124062      .25635718960E-01
  .28266888959      .16606859208E-01
S  1  1.00
  .29058164020E-01  1.0000000000
S  1  1.00
  .12111638157E-01  1.0000000000
P  5  1.00
  361.22492154      .20906479823E-01
  84.670222166      .15043641740
  26.469088236      .55440061077
  9.2658077615      1.0409009991
  3.3423388293      .67825341194
P  3  1.00
  1.5100876104      .75248191146
  .56568375163      1.3708585031
  .20817008495      .66047633079
P  1  1.00
  .041737   1.000000
D  1  1.00
  0.353            1.0
D  1  1.00
  0.098            1.0
****
H O 0
6-31G*
****

202208241220409166..png (9.32 KB, 下载次数 Times of downloads: 11)

202208241220409166..png

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 09:06 , Processed in 0.281438 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list