计算化学公社

标题: 求助：如何随机抽取分子振动模式 [打印本页]

作者
Author: xdyjiay 时间: 2024-6-21 14:19
标题: 求助：如何随机抽取分子振动模式
想请问一下各位老师，假如有n个原子，有3n-6个振动模式，从这3n-6个振动模式中随机抽取一个振动模式，应该怎么抽取比较合适？

作者
Author: wzkchem5 时间: 2024-6-21 15:14
随机抽取一个振动模式不难，产生一个1到3n-6之间的随机数就行了，bash、python乃至excel公式都能产生。
但是很多时候需要的不是抽取一个振动模式，而是需要随机抽取各种振动模式的一个线性组合，此时需要写一个很短的脚本

作者
Author: xdyjiay 时间: 2024-6-21 16:53

wzkchem5 发表于 2024-6-21 15:14
随机抽取一个振动模式不难，产生一个1到3n-6之间的随机数就行了，bash、python乃至excel公式都能产生。
但 ...

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

作者
Author: wzkchem5 时间: 2024-6-21 19:16

xdyjiay 发表于 2024-6-21 09:53
老师，随机抽取一个振动模式是要遵循高斯分布来抽，还是认为每个振动模式概率一样来抽？
还有一个问题， ...

每个振动模式都要产生一个随机数。建议看一下多维高斯分布的抽样方法

作者
Author: Daniel_Arndt 时间: 2024-6-22 10:43
Daniel Lambrecht在 https://doi.org/10.1021/jp110334w 这篇文章的Appendix部分写的内容可能对你有用。

作者
Author: xdyjiay 时间: 2024-6-24 15:56

wzkchem5 发表于 2024-6-21 19:16
每个振动模式都要产生一个随机数。建议看一下多维高斯分布的抽样方法

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

作者
Author: wzkchem5 时间: 2024-6-24 16:31

xdyjiay 发表于 2024-6-24 08:56
想请问一下如果是根据高斯分布抽取振动模式的随机数，高斯分布的具体函数应该是什么样的啊？

仔细阅读https://en.wikipedia.org/wiki/Multivariate_normal_distribution

欢迎光临计算化学公社 (http://bbs.keinsci.com/)