计算化学公社

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

[Gaussian/gview] 求助:如何随机抽取分子振动模式

[复制链接 Copy URL]

23

帖子

0

威望

277

eV
积分
300

Level 3 能力者

想请问一下各位老师,假如有n个原子,有3n-6个振动模式,从这3n-6个振动模式中随机抽取一个振动模式,应该怎么抽取比较合适?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

2#
发表于 Post on 2024-6-21 15:14:45 | 只看该作者 Only view this author
随机抽取一个振动模式不难,产生一个1到3n-6之间的随机数就行了,bash、python乃至excel公式都能产生。
但是很多时候需要的不是抽取一个振动模式,而是需要随机抽取各种振动模式的一个线性组合,此时需要写一个很短的脚本
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

23

帖子

0

威望

277

eV
积分
300

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-6-21 16:53:39 | 只看该作者 Only view this author
wzkchem5 发表于 2024-6-21 15:14
随机抽取一个振动模式不难,产生一个1到3n-6之间的随机数就行了,bash、python乃至excel公式都能产生。
但 ...

老师,随机抽取一个振动模式是要遵循高斯分布来抽,还是认为每个振动模式概率一样来抽?
还有一个问题,确实是像您说的需要抽取各种振动模式的线性组合,根据文献的描述,我理解的是从1到3n-6抽出一个随机数乘以正则坐标,每个正则坐标都重复上面的操作,最后叠加产生一个位移矢量,不晓得我理解的正确吗?

c1fb0a33a3decb6de8f7e8a830e5373.png (74.65 KB, 下载次数 Times of downloads: 4)

c1fb0a33a3decb6de8f7e8a830e5373.png

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

4#
发表于 Post on 2024-6-21 19:16:44 | 只看该作者 Only view this author
xdyjiay 发表于 2024-6-21 09:53
老师,随机抽取一个振动模式是要遵循高斯分布来抽,还是认为每个振动模式概率一样来抽?
还有一个问题, ...

每个振动模式都要产生一个随机数。建议看一下多维高斯分布的抽样方法
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

417

帖子

1

威望

2196

eV
积分
2633

Level 5 (御坂)

5#
发表于 Post on 2024-6-22 10:43:03 | 只看该作者 Only view this author
Daniel Lambrecht在 https://doi.org/10.1021/jp110334w 这篇文章的Appendix部分写的内容可能对你有用。

23

帖子

0

威望

277

eV
积分
300

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2024-6-24 15:56:10 | 只看该作者 Only view this author
wzkchem5 发表于 2024-6-21 19:16
每个振动模式都要产生一个随机数。建议看一下多维高斯分布的抽样方法

想请问一下如果是根据高斯分布抽取振动模式的随机数,高斯分布的具体函数应该是什么样的啊?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

7#
发表于 Post on 2024-6-24 16:31:57 | 只看该作者 Only view this author
xdyjiay 发表于 2024-6-24 08:56
想请问一下如果是根据高斯分布抽取振动模式的随机数,高斯分布的具体函数应该是什么样的啊?

仔细阅读https://en.wikipedia.org/wiki/Multivariate_normal_distribution
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 10:30 , Processed in 0.175435 second(s), 24 queries , Gzip On.

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