计算化学公社

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

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

   关闭 [复制链接 Copy URL]

200

帖子

1

威望

1325

eV
积分
1545

Level 5 (御坂)

15106#
发表于 Post on 2022-4-4 23:28:04 | 只看该作者 Only view this author
请问下各位老师,我笔记本电脑上同时安装了G16 B01,GaussView 6和G09 D01,GaussView 5两个版本的高斯和GaussView。请问我用GaussView 5打开GaussView 6产生的gif文件,或者用GaussView 6打开GaussView 5产生的gjf文件,会出现问题(比如文本内容变化等)吗?现在有点虚......
计算小白 水平很菜
随时可喷 万望指点

1592

帖子

0

威望

5155

eV
积分
6747

Level 6 (一方通行)

15107#
发表于 Post on 2022-4-5 08:46:14 | 只看该作者 Only view this author
WaterMirror 发表于 2022-4-4 23:28
请问下各位老师,我笔记本电脑上同时安装了G16 B01,GaussView 6和G09 D01,GaussView 5两个版本的高斯和Gau ...

http://sobereva.com/439

  1. 应当注意Gaussian和GaussView的版本兼容性。根据笔者经验,Gaussian 09建议搭配GaussView 5.0.9使用,Gaussian 16建议搭配GaussView 6.0.16(目前的最新版本)使用,这样兼容性是最好的,否则会有很多问题。比如GaussView 6.0.16打开Gaussian 09的IRC任务的帧号顺序是错乱的,而GaussView 5.0.9无法正确载入Gaussian 16的振动分析输出文件等。
复制代码
又菜又爱玩

412

帖子

0

威望

2577

eV
积分
2989

Level 5 (御坂)

15108#
发表于 Post on 2022-4-5 16:02:40 | 只看该作者 Only view this author
请问一下,有哪些方法可以计算电负性

6万

帖子

99

威望

6万

eV
积分
126501

管理员

公社社长

15109#
 楼主 Author| 发表于 Post on 2022-4-5 21:46:27 | 只看该作者 Only view this author
mfdsrax2 发表于 2022-4-5 16:02
请问一下,有哪些方法可以计算电负性

看具体什么体系
如果是孤立体系,直接按标准Mulliken电负性基于IP和EA算。
北京科音自然科学研究中心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

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15110#
发表于 Post on 2022-4-6 20:44:25 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-4 18:29
STEOM-DLPNO-CCSD对基组还是比较敏感的,毕竟是post-HF方法,比DFT要敏感得多。建议可以试试把基组加大到 ...

老师,我将基组提升到def2-TZVP(-f),三重态的激发能提升了0.05eV
我还要加大基组么,但是现在提示Percentage Active Character     80.30

1万

帖子

0

威望

1万

eV
积分
22432

Level 6 (一方通行)

15111#
发表于 Post on 2022-4-6 21:20:31 | 只看该作者 Only view this author
心向暖阳 发表于 2022-4-6 13:44
老师,我将基组提升到def2-TZVP(-f),三重态的激发能提升了0.05eV
我还要加大基组么,但是现在提示Perce ...

那估计就不是基组问题而是你说的active space太小的问题了,可以调一下控制active space的阈值,多包括一些活性轨道进来
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
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15112#
发表于 Post on 2022-4-7 09:24:25 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-6 21:20
那估计就不是基组问题而是你说的active space太小的问题了,可以调一下控制active space的阈值,多包括一 ...

One can also select the active spaces manually by turning the DoCISNat to false and setting the NActIP
and NActEA (RHF STEOM calculation) or the NActIP a, NActIP b, NActEA a and NActEA b (UHF STEOM
calculation) to desired values. However, this is not recommended for general users. The following shows the
output of the active orbital selection procedure on a closed-shell molecule:

我这个该设置多大呢
DoCISNat false
NActIP N
NActEA N

还有就是下面比较详细的提示,还望老师看看
------------------------------------------
STATE AVERAGED NATURAL ORBITALS FOR ACTIVE SPACE SELECTION
------------------------------------------

Solving eigenvalue problem for the occupied space  ... done
Solving eigenvalue problem for the virtual space   ... done
After automatic selection of the active space:
No of roots active in IP calculation:    8  
No of roots active in EA calculation:   12  
Testing the validity of the selected OCC active space ...
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors

WARNING: I will add 1 orbitals to active OCC space !!

Testing the validity of the selected VIRT active space ...
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors
  --- complex eigenvalues and eigenvectors

WARNING: I will add 6 orbitals to active VIRT space !!

After auto-adaptation of the active space:
No of roots active in IP calculation:    9  
No of roots active in EA calculation:   18  
done (  362.7)

---------------------------------------------------

1万

帖子

0

威望

1万

eV
积分
22432

Level 6 (一方通行)

15113#
发表于 Post on 2022-4-7 14:29:08 | 只看该作者 Only view this author
心向暖阳 发表于 2022-4-7 02:24
One can also select the active spaces manually by turning the DoCISNat to false and setting the NA ...

这个没有一定规则,比程序自动选的轨道多就行
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
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

1646

eV
积分
1714

Level 5 (御坂)

15114#
发表于 Post on 2022-4-7 15:06:42 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-7 14:29
这个没有一定规则,比程序自动选的轨道多就行

好的 谢谢老师

412

帖子

0

威望

2577

eV
积分
2989

Level 5 (御坂)

15115#
发表于 Post on 2022-4-8 11:12:29 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-31 17:02
可能有,可能没有
对断裂的那根键做一个柔性扫描,其中最长的键长至少等于两个原子范德华半径之和。如果 ...

C-Cl范德华半径之和大概是3.5左右,不同计算方法稍有差异,我从1.4扫到了3.4,还需要继续扫吗?

202204081111164461..png (21.04 KB, 下载次数 Times of downloads: 21)

202204081111164461..png

1万

帖子

0

威望

1万

eV
积分
22432

Level 6 (一方通行)

15116#
发表于 Post on 2022-4-8 14:18:54 | 只看该作者 Only view this author
mfdsrax2 发表于 2022-4-8 04:12
C-Cl范德华半径之和大概是3.5左右,不同计算方法稍有差异,我从1.4扫到了3.4,还需要继续扫吗?[/backcol ...

扫到范德华半径之和就不需要扫了,虽然这样不能严格保证没有过渡态,但是如果过渡态里面的键长超过了范德华半径之和,就可以认为这是一个协同异步反应,其中过渡态的时候键还完全没有形成,那么此时就可以合理地说起码键的形成/断裂这个过程没有过渡态,过渡态是其他非共价作用导致的
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
本团队长期招收研究生,有意者可私信联系

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

15117#
发表于 Post on 2022-4-9 17:51:52 | 只看该作者 Only view this author
请教老师,我对m062x泛函优化好的结构采用CBS-QB3计算能量,为防止结构变化,我固定了分子坐标,out文件中ZPE=0,原本用0k-ZPE计算能量,那么现在我该如何得到能量呢?谢谢老师

202204091751062848..png (16.99 KB, 下载次数 Times of downloads: 20)

202204091751062848..png

6万

帖子

99

威望

6万

eV
积分
126501

管理员

公社社长

15118#
 楼主 Author| 发表于 Post on 2022-4-11 00:31:29 | 只看该作者 Only view this author
糖糖糖豆9988 发表于 2022-4-9 17:51
请教老师,我对m062x泛函优化好的结构采用CBS-QB3计算能量,为防止结构变化,我固定了分子坐标,out文件中Z ...

先说清楚你要得到什么
你如果要得到电子能量,就不用管ZPE这些热力学校正量
如果要以混合方式得到热力学量,把M06-2X做opt freq得到的热力学校正量加到CBS-QB3的电子能量上
如果要要以一般方式得到CBS-QB3的热力学量,事先根本不需要M06-2X优化,直接让CBS-QB3一条龙做完
北京科音自然科学研究中心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

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

15119#
发表于 Post on 2022-4-11 14:43:16 | 只看该作者 Only view this author
本帖最后由 糖糖糖豆9988 于 2022-4-11 15:36 编辑
sobereva 发表于 2022-4-11 00:31
先说清楚你要得到什么
你如果要得到电子能量,就不用管ZPE这些热力学校正量
如果要以混合方式得到热力 ...

谢谢sob老师的回复,再追问您一些简单的问题,希望您不要嫌弃。①CBS-QB3方法的CBS-QB3(0K)=SPE + ZPE,所以之前用CBS-QB3(0K)- E(ZPE)得到电子能量,如果我只想要得到电子能量,那么现在直接采用CBS-QB3(0K)的值是不是就可以啦?②如果我就是想要以第二种方法得到热力学量,我之前已经用m062x 6-311G(d,p)做过opt+freq计算,就从相应输出文件把“Zero-point correction=                           0.212262 (Hartree/Particle)”和CBS-QB3(0K)的值相加相加即可吗?③您说的热力学量是指什么呢?和那个U(0)=Eele+Ezpe是一样的吗?麻烦老师您解释一下,谢谢。

6万

帖子

99

威望

6万

eV
积分
126501

管理员

公社社长

15120#
 楼主 Author| 发表于 Post on 2022-4-11 18:48:31 | 只看该作者 Only view this author
糖糖糖豆9988 发表于 2022-4-11 14:43
谢谢sob老师的回复,再追问您一些简单的问题,希望您不要嫌弃。①CBS-QB3方法的CBS-QB3(0K)=SPE + ZPE, ...

1 显然得把CBS-QB3(0K)减去ZPE
2 显然是错的。CBS-QB3(0K)里本来就包含ZPE
3 把下文,以及Shermo手册附录看了,了解热力学计算基本常识
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552http://bbs.keinsci.com/thread-17494-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

本版积分规则 Credits rule

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

GMT+8, 2026-4-27 06:25 , Processed in 0.206651 second(s), 23 queries , Gzip On.

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