求助：请问如何MD模拟得到分子偶极矩自相关函数，从而求得IR光谱

wildon · 发表于 Post on 2023-5-3 22:02:09

请问如何通过lammps的模拟得到分子偶极矩变化，并计算自相关函数的傅里叶变化，从而得到红外光谱？
目前了解可以通过量化程序或基于第一性原理的程序进行IR的计算，但我想通过lammps进行MD的模拟求得IR分布。

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

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

wildon · 发表于 Post on 2023-5-3 22:04:41

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

sobereva · 发表于 Post on 2023-5-4 05:35:15

直接就有现成的程序

wildon · 发表于 Post on 2023-5-4 09:09:40

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

谢谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Lammps] 求助：请问如何MD模拟得到分子偶极矩自相关函数，从而求得IR光谱

评分 Rate

浏览过的版块