计算化学公社

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

[算法与编程] 基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq

[复制链接 Copy URL]

195

帖子

0

威望

4029

eV
积分
4224

Level 6 (一方通行)

16#
发表于 Post on 2021-7-25 10:23:50 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-24 21:50
normal coordinates只有方向是有物理意义的,把整个normal coordinate矢量取负号,不改变任何物理可观测 ...

明白了,谢谢

170

帖子

2

威望

1734

eV
积分
1944

Level 5 (御坂)

17#
发表于 Post on 2022-6-8 16:57:28 | 只看该作者 Only view this author
sob老师您好!
说明文档中提到:“Gaussian输出最终的振动频率前,会先从3N个频率中自动判断出整体平动和转动模式来,然后将它们对力常数矩阵的贡献投影掉使得振动频率更精确。”这里是如何判断整体平动转动模式的呢,是否也是和这个程序中一样去掉最小的六个。另外这里的投影可以给一下具体的数学表达吗?最近在做一些工作,想把一些笛卡尔坐标系下的矢量转换到内坐标下,对如何删除平动转动自由度感到很是困惑。

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

18#
发表于 Post on 2022-6-8 18:08:30 | 只看该作者 Only view this author
Kalinite 发表于 2022-6-8 09:57
sob老师您好!
说明文档中提到:“Gaussian输出最终的振动频率前,会先从3N个频率中自动判断出整体平动和 ...

我不知道sob老师的做法是什么,但是一般量化程序的做法并不是对角化以后去掉最小的6个频率,而是先推导出平动和转动模式的解析形式,再把这些模式从Hessian里投影掉,然后对角化剩下的3N-6维的Hessian。这样的好处是假如有振动频率极其小的振动模式,不至于误把这个振动扔掉。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

170

帖子

2

威望

1734

eV
积分
1944

Level 5 (御坂)

19#
发表于 Post on 2022-6-8 20:49:54 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-8 18:08
我不知道sob老师的做法是什么,但是一般量化程序的做法并不是对角化以后去掉最小的6个频率,而是先推导出 ...

请问您有关于“推导出平动和转动模式的解析形式”这方面的资料吗

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

20#
发表于 Post on 2022-6-8 20:56:16 | 只看该作者 Only view this author
Kalinite 发表于 2022-6-8 13:49
请问您有关于“推导出平动和转动模式的解析形式”这方面的资料吗

没有,这个就是经典力学的内容,找一本力学教材看看刚体平动和转动的简正模怎么算就知道了
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

21#
 楼主 Author| 发表于 Post on 2022-6-9 07:38:06 | 只看该作者 Only view this author
Kalinite 发表于 2022-6-8 20:49
请问您有关于“推导出平动和转动模式的解析形式”这方面的资料吗

Gaussian有个官方文档Vibrational Analysis in Gaussian,里面有细节
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

187

帖子

1

威望

1247

eV
积分
1454

Level 4 (黑子)

22#
发表于 Post on 2022-6-9 08:24:13 | 只看该作者 Only view this author
Gaussian官方文档Vibrational Analysis in Gaussian

Vibrational.pdf

149.04 KB, 下载次数 Times of downloads: 17

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
snljty2 + 5 谢谢

查看全部评分 View all ratings

计算小白 水平很菜
随时可喷 万望指点

170

帖子

2

威望

1734

eV
积分
1944

Level 5 (御坂)

23#
发表于 Post on 2022-6-9 13:27:26 | 只看该作者 Only view this author
感谢大家的帮助!

104

帖子

1

威望

3519

eV
积分
3643

Level 5 (御坂)

24#
发表于 Post on 2022-6-10 11:41:15 | 只看该作者 Only view this author
推荐 ajz34 的讲解和代码
  1. 频率分析 (1):从 fchk 文件得到分子频率与简正模式
  2. https://ajz34.readthedocs.io/zh_CN/latest/QC_Notes/Freq_Series/freq_1.html
复制代码

另推荐他写的 python 实现双杂化泛函梯度教程,其中的量化入门部分很友好
  1. xDH 在 Python 下实现的简易教程
  2. https://py-xdh.readthedocs.io/
复制代码

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
风起~ + 3 好物!

查看全部评分 View all ratings

十八介姑娘一蕾花呀,白白介牙齿、红红介嘴唇,得人惜

89

帖子

0

威望

1711

eV
积分
1800

Level 5 (御坂)

25#
发表于 Post on 2022-6-10 11:49:59 | 只看该作者 Only view this author
王二葛 发表于 2022-6-10 11:41
推荐 ajz34 的讲解和代码

另推荐他写的 python 实现双杂化泛函梯度教程,其中的量化入门部分很友好

谢谢

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 17:32 , Processed in 0.181705 second(s), 24 queries , Gzip On.

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