VASP算完IR该如何处理得到谱图？

菜鸡互啄 · 发表于 Post on 2019-11-18 12:25:31

本帖最后由菜鸡互啄于 2019-11-18 18:07 编辑

各位老师好，

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

dummy1 · 发表于 Post on 2019-11-18 12:36:45

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

ORCA_in_TCC · 发表于 Post on 2019-11-18 15:02:59

有一个后处理的shell脚本，可以得到频率对应的强度。可以先找找

菜鸡互啄 · 发表于 Post on 2019-11-18 17:56:04

本帖最后由菜鸡互啄于 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

菜鸡互啄 · 发表于 Post on 2019-11-18 17:59:47

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

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

卡开发发 · 发表于 Post on 2019-11-18 21:25:10

菜鸡互啄发表于 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 的部分。

菜鸡互啄 · 发表于 Post on 2019-11-19 09:10:29

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

谢谢大佬，我试试看

ORCA_in_TCC · 发表于 Post on 2019-11-21 16:49:44

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

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

啦啦黑还黑 · 发表于 Post on 2019-11-25 22:28:47

有现成的脚本

sky8fall · 发表于 Post on 2023-9-14 17:29:27

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

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

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

[VASP] VASP算完IR该如何处理得到谱图？

评分 Rate

浏览过的版块