应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP频率分析
252/3返回列表上一页123下一页
查看: 6753  |  回复: 24
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lgxyz

铜虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by leila3409 at 2011-02-23 06:38:19:
这跟我第一次计算时的设置是一样的，呵呵

那是怎么回事？虚频消不掉了吗？
赞一下(1人) 回复此楼
11楼2011-02-23 14:42:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
Originally posted by cenwanglai at 2011-02-23 05:08:51:
呵呵，可能是我没有说清楚我的意思。

不是1f就代表z方向。

我是猜测你的POSCAR红，两个H原子的x,y坐标是相同的，两个原子沿着z方向排列。要不你帖出POSCAR看看？

出现这么多频率是程序自己折腾的，程序 ...

恩，版主说的没错，我确实是沿Z方向放的，POSCAR是

H2:                                    
   10     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
2
Selective dynamics
Direct
0 0 0      T   T   T
0 0 0.075 T T T

可是，我还是有点不放心那些2f,3f...，它们的出现真的可以被无视么？万一两个数据都比较接近实验值呢？我按版主前面提的矢量和得到的是4274cm-1，跟实验值更接近了呢。
所以说，我应该取哪个数据呢？1f还是矢量和呢？
赞一下 回复此楼
12楼2011-02-23 14:44:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
Originally posted by lgxyz at 2011-02-23 06:42:16:
那是怎么回事？虚频消不掉了吗？

我第一次算频率的东西，说实施虚频代表什么我都不知道。会不会像版主说的那是分子的整体振动引起的？
赞一下 回复此楼
13楼2011-02-23 14:46:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指教 2011-02-23 21:08:30
引用回帖:
Originally posted by leila3409 at 2011-02-23 06:46:15:
我第一次算频率的东西，说实施虚频代表什么我都不知道。会不会像版主说的那是分子的整体振动引起的？

实频就是说原子在某个方向在固定位置上进行振动，而虚频就是说原子在某个方向虽然是平衡态，但是是不稳定的。虚频如果很小，是可以忽略的，因为计算是有误差的

[ Last edited by bingmou on 2011-2-23 at 17:09 ]
赞一下(2人) 回复此楼
14楼2011-02-23 17:07:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+1): 谢谢交流。 2011-02-24 11:59:34
只有第一个值是有意义的吧，因为H2分子只有一个伸缩振动啊
赞一下(2人) 回复此楼
15楼2011-02-24 11:51:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢讨论 2011-02-25 08:32:24
2个原子，6个自由度，当然就是6支频率。但你计算的虚频可以忽略，非常小，可能是误差。没误差的话，最后三个频率值应该为0.

[ Last edited by stractor on 2011-2-25 at 08:05 ]
赞一下(3人) 回复此楼
16楼2011-02-25 08:00:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿黛拉

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
能否贴出INCAR文件。。。
赞一下(1人) 回复此楼
进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。
17楼2011-06-03 09:17:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿黛拉

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by leila3409 at 2011-02-23 13:00:06:
我按照你说的做了设置，计算的结果是

   1 f  =  127.818139 THz   803.105052 2PiTHz 4263.554054 cm-1   528.613521 meV
   2 f  =    3.840658 THz    24.131568 2PiTHz  128.110567 cm-1    15.883692  ...

请问vasp参数设置不同而优化得到的能量相同，但是频率相差很大？为什么呢？
赞一下(1人) 回复此楼
进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。
18楼2011-06-03 09:23:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)
★ ★
gzqdyouxia(金币+2): 谢谢提示 2011-06-03 13:07:08
后面的dx，dy，dz是振动模式和方向，按照你的计算，第一个是沿着键的方向的振动频率，第二，三个是垂直于键的振动频率，最后3个是平移的频率，理论上说应该是0
vasp的频率计算实际上是gamma点处的声子谱，即实际上是周期型结构的震动频率，所以才会有多少个自由度就会有多少频率，在这里。最后3个是声学波，为0
赞一下(1人) 回复此楼
19楼2011-06-03 13:04:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
gzqdyouxia(金币+2): 谢谢指教 2011-06-03 16:19:53
频率下面给出了本征矢量，分析本征矢量可以知道这个频率对应于什么类型的振动。第一支频率就是2个H原子的相对振动，即你所要的振动频率。因为一个分子的振动模式数是3N，N为原子个数，所以这里有6支频率。后面的5支频率是旋转和刚性平移的频率，为0。
    你计算的值和实验值差别较大，可能是如下原因：
最可能的原因出在POTIM。 POTIM=0.015太小了，这样计算的位移后构型的能量和位移前构型的能量可能没差别哦，可以设为POTIM=0.4。你可以查查资料，有限差分法计算频率就是能量对位置求二次导数，你POTIM太小，VASP计算的位移和没位移的能量没明显的差别。理想的POTIM为0.4-1.0。
    计算频率时H原子间距是取的优化的还是实验的？
    你是采用中心差分法吗？即NFREE=2吗？  最好把INCAR贴出来。

[ Last edited by stractor on 2011-6-4 at 02:05 ]
赞一下(3人) 回复此楼
20楼2011-06-03 15:45:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leila3409 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭