计算化学公社

标题: VASP算完IR该如何处理得到谱图？ [打印本页]

作者
Author: 菜鸡互啄 时间: 2019-11-18 12:25
标题: VASP算完IR该如何处理得到谱图？
本帖最后由菜鸡互啄于 2019-11-18 18:07 编辑

各位老师好，

最近使用VASP计算了材料的IR，通过脚本可以得到峰的位置和相对强度：001 3556.5 1.000像这种，但是我想画成有宽度的峰，请问该如何处理才能得到红外谱图呢？我想画成图跟实验结果对比。请各位老师指点迷津，非常感谢！

作者
Author: dummy1 时间: 2019-11-18 12:36
不是很理解你的意思，用脚本得到了位置和强度不就是红外谱图吗? 还是说你得到是一条一条的线，但想画成有一定宽度的峰？是的话你可以用高斯或者洛伦兹展宽。如果你跑的AIMD，谱本身就是展宽的了，可以直接和实验比对。

作者
Author: ORCA_in_TCC 时间: 2019-11-18 15:02
有一个后处理的shell脚本，可以得到频率对应的强度。可以先找找

作者
Author: 菜鸡互啄 时间: 2019-11-18 17:56
本帖最后由菜鸡互啄于 2019-11-18 18:15 编辑

dummy1 发表于 2019-11-18 12:36
不是很理解你的意思，用脚本得到了位置和强度不就是红外谱图吗? 还是说你得到是一条一条的线，但想画成有一 ...

我用脚本得到的是一个峰位置和一个相对强度值：001 3556.5 1.000这种，但是我想画成有一定宽度的峰。是不用脚本就可以从某个文件中读出谱图吗？我的计算和后处理方法是从这里学到的https://homepage.univie.ac.at/david.karhanek/downloads.html#Entry02

作者
Author: 菜鸡互啄 时间: 2019-11-18 17:59

ORCA_in_TCC 发表于 2019-11-18 15:02
有一个后处理的shell脚本，可以得到频率对应的强度。可以先找找

我就是用的脚本，得到频率和相对强度，但是我想画有一定宽度的那种谱图。

作者
Author: 卡开发发 时间: 2019-11-18 21:25

菜鸡互啄发表于 2019-11-18 17:59
我就是用的脚本，得到频率和相对强度，但是我想画有一定宽度的那种谱图。

展宽过程应该不是很复杂啊，假定未展宽的谱线位置在w_i处，强度为A_i，那么展宽后的谱线F(w)=∑_i Ai*M((w-wi),σ)
M是展宽函数，比如Gaussian函数的话就M((w-wi),σ)=1/a*exp(-(w-wi)^2/2σ^2)。这个过程类似于Multiwfn手册3.12.1 的部分。

作者
Author: 菜鸡互啄 时间: 2019-11-19 09:10

卡开发发发表于 2019-11-18 21:25
展宽过程应该不是很复杂啊，假定未展宽的谱线位置在w_i处，强度为A_i，那么展宽后的谱线F(w)=∑_i Ai*M(( ...

谢谢大佬，我试试看

作者
Author: ORCA_in_TCC 时间: 2019-11-21 16:49

菜鸡互啄发表于 2019-11-19 09:10
谢谢大佬，我试试看

你说的是展宽问题，这是另外一个简单的问题，再找找吧

作者
Author: 啦啦黑还黑 时间: 2019-11-25 22:28
有现成的脚本

作者
Author: sky8fall 时间: 2023-9-14 17:29

菜鸡互啄发表于 2019-11-18 17:56
我用脚本得到的是一个峰位置和一个相对强度值：001 3556.5 1.000这种，但是我想画成有一定宽度的峰。是不 ...

能请问一下如何获得峰强吗？这个链接失效了

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