计算化学公社

标题: 求助：请问如何MD模拟得到分子偶极矩自相关函数，从而求得IR光谱 [打印本页]

作者
Author: wildon 时间: 2023-5-3 22:02
标题: 求助：请问如何MD模拟得到分子偶极矩自相关函数，从而求得IR光谱
请问如何通过lammps的模拟得到分子偶极矩变化，并计算自相关函数的傅里叶变化，从而得到红外光谱？
目前了解可以通过量化程序或基于第一性原理的程序进行IR的计算，但我想通过lammps进行MD的模拟求得IR分布。

思路为:
1. 通过lammps中的 compute dipole 计算整体或某个方向上的偶极矩变化；
2. 求第一步偶极矩变化的自相关函数分布，并进行傅里叶变化，这样能得到IR分布嘛？
PS：第2步是不是要进行展宽的处理？

看了一晚上的红外光谱计算，基于MD计算红外光谱资料比较少，了解的不够深入，请有相关的老师不吝赐教，谢谢！

作者
Author: wildon 时间: 2023-5-3 22:04
目前的进展为：通过lammps中的 compute dipole 计算了电场方向的偶极矩变化，通过python代码计算了偶极矩随时间变化的自相关函数。
后续应该如何做呢。
再对自相关函数分布进行傅里叶变化吗？

作者
Author: sobereva 时间: 2023-5-4 05:35
直接就有现成的程序

(, 下载次数 Times of downloads: 40)

作者
Author: wildon 时间: 2023-5-4 09:09

sobereva 发表于 2023-5-4 05:35
直接就有现成的程序

谢谢老师

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