计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6087|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 采用关键词freq=intmodes,分析键长、键角、二面角对重组能贡献的疑惑

[复制链接 Copy URL]

5

帖子

0

威望

351

eV
积分
356

Level 3 能力者

如题,看到一篇分析分子的键长键角二面角对重组能的贡献的文献,文中没说处理方法。按照之前sob老师回复的帖子(如何分析键长、键角、二面角对重组能的贡献?http://bbs.keinsci.com/forum.php ... 26&fromuid=9185(出处: 计算化学公社))做了测试,有以下几个疑惑:1.每个正则模式下键长键角二面角的relative weight 总和为什么不是1?  2.具体算法是否有误:我是先用dushin计算分解了总的重组能,之后,分析每一个正则模式下,键长键角二面角的总占比(根据relative weight权重的占比计算的,不是按照1算的,总的加起来是多少就除多少),然后在找到该模式下对应的对重组能的贡献,就的到该模式下,键长键角二面角的贡献。这样就得到了所有的键长键角二面角对重组能的贡献。但是结果相差很多,尤其是键长的。
PS:我做的是THF里的。感谢各位老师同学解惑。

3

帖子

0

威望

155

eV
积分
158

Level 3 能力者

8#
发表于 Post on 2025-11-27 19:22:52 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-11-27 10:53
https://gaussian.com/overlay7

设IOp(7/75)=-1取消冗余内坐标对正则模式贡献的打印阈值。

非常感谢!

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

7#
发表于 Post on 2025-11-27 10:53:04 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2025-11-27 20:21 编辑
埃布-wang 发表于 2025-11-27 10:06
sob老师您好,您提到贡献大于特定值的项才会输出,请问如何写关键词让所有的项都输出呢?

https://gaussian.com/overlay7
IOp(7/75)

Threshold for printing redundant internal contributions to normal mode displacements.

0        Default (10%).
N        10-N.
-1        Zero (all printed).
The threshold is automatically lowered for each mode until 90% of the absolute displacements are included.

用IOp(7/75=-1)设IOp(7/75)为-1取消冗余内坐标对正则模式贡献的打印阈值。编辑:修复了括号的问题

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
wal + 1 是IOp(7/75=-1)吧

查看全部评分 View all ratings

√546=23.36664289109

3

帖子

0

威望

155

eV
积分
158

Level 3 能力者

6#
发表于 Post on 2025-11-27 10:06:52 | 只看该作者 Only view this author
sobereva 发表于 2021-9-11 03:21
1 贡献大于特定值的项才会输出。所以你会发现每个振动模式显示的项数往往是不同的
2 可能作者把各个冗余内 ...

sob老师您好,您提到贡献大于特定值的项才会输出,请问如何写关键词让所有的项都输出呢?

5

帖子

0

威望

351

eV
积分
356

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-9-15 16:09:47 | 只看该作者 Only view this author
sobereva 发表于 2021-9-15 01:31
本来默认就是在冗余内坐标下做的,intmodes给的也是冗余内坐标的贡献
考虑力常数的差异更合理

好哒,谢谢sob老师!

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

4#
发表于 Post on 2021-9-15 01:31:42 | 只看该作者 Only view this author
qinxiaodao 发表于 2021-9-14 10:18
非常感谢sob老师,很抱歉没能及时感谢您。我还有两个问题想请教您:考虑冗余内坐标的话是不是在关键词位 ...

本来默认就是在冗余内坐标下做的,intmodes给的也是冗余内坐标的贡献
考虑力常数的差异更合理
北京科音自然科学研究中心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

5

帖子

0

威望

351

eV
积分
356

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-9-14 10:18:58 | 只看该作者 Only view this author
sobereva 发表于 2021-9-11 03:21
1 贡献大于特定值的项才会输出。所以你会发现每个振动模式显示的项数往往是不同的
2 可能作者把各个冗余内 ...

非常感谢sob老师,很抱歉没能及时感谢您。我还有两个问题想请教您:考虑冗余内坐标的话是不是在关键词位置处添加opt=redundant呢?对于这种分析键长键角对重组能的贡献的研究,那种更合理一些呢?想确认一下我测试的方法是不是合理的。再次感谢sob老师!

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2021-9-11 03:21:33 | 只看该作者 Only view this author
1 贡献大于特定值的项才会输出。所以你会发现每个振动模式显示的项数往往是不同的
2 可能作者把各个冗余内坐标的力常数的差异问题考虑了进去
北京科音自然科学研究中心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

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

GMT+8, 2026-1-23 22:09 , Processed in 0.173382 second(s), 25 queries , Gzip On.

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