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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP频率分析
52/1返回列表
查看: 6797  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

leila3409

捐助贵宾 (正式写手)

[交流] 【求助】VASP频率分析 已有11人参与

请问高手们，我想得到氢分子的振动频率vibrational frequency（cm-1），设置了NSW=1, IBRION=5, POTIM=0.015，OUTCAR中得到下面的结果，请问哪个才是我所需要的氢分子的频率？
1 f  =  127.990789 THz   804.189843 2PiTHz 4269.313033 cm-1   529.327544 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.000036   -0.000036    0.707238  
      0.000000  0.000000  0.750000     0.000283   -0.000011   -0.706976  

   2 f  =    4.690886 THz    29.473704 2PiTHz  156.471101 cm-1    19.399951 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.748316   -0.000024    0.000172  
      0.000000  0.000000  0.750000    -0.663342    0.000014   -0.000131  

   3 f  =    4.396660 THz    27.625032 2PiTHz  146.656802 cm-1    18.183133 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.000022   -0.759554   -0.000017  
      0.000000  0.000000  0.750000     0.000015    0.650444    0.000012  

   4 f/i=    0.023611 THz     0.148351 2PiTHz    0.787570 cm-1     0.097646 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.004764   -0.000534   -0.706959  
      0.000000  0.000000  0.750000     0.005330   -0.000629   -0.707218  

   5 f/i=    0.283258 THz     1.779760 2PiTHz    9.448456 cm-1     1.171460 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.663323    0.001788    0.004920  
      0.000000  0.000000  0.750000     0.748294    0.002094    0.005187  

   6 f/i=    0.340255 THz     2.137885 2PiTHz   11.349685 cm-1     1.407182 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.001819    0.650441   -0.000575  
      0.000000  0.000000  0.750000    -0.002059    0.759551   -0.000620  

Eigenvectors after division by SQRT(mass)

Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   1 f  =  127.990789 THz   804.189843 2PiTHz 4269.313033 cm-1   529.327544 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.000036   -0.000036    0.707238  
      0.000000  0.000000  0.750000     0.000283   -0.000011   -0.706976  

   2 f  =    4.690886 THz    29.473704 2PiTHz  156.471101 cm-1    19.399951 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.748316   -0.000024    0.000172  
      0.000000  0.000000  0.750000    -0.663342    0.000014   -0.000131  

   3 f  =    4.396660 THz    27.625032 2PiTHz  146.656802 cm-1    18.183133 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.000022   -0.759554   -0.000017  
      0.000000  0.000000  0.750000     0.000015    0.650444    0.000012  

   4 f/i=    0.023611 THz     0.148351 2PiTHz    0.787570 cm-1     0.097646 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.004764   -0.000534   -0.706959  
      0.000000  0.000000  0.750000     0.005330   -0.000629   -0.707218  

   5 f/i=    0.283258 THz     1.779760 2PiTHz    9.448456 cm-1     1.171460 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.663323    0.001788    0.004920  
      0.000000  0.000000  0.750000     0.748294    0.002094    0.005187  

   6 f/i=    0.340255 THz     2.137885 2PiTHz   11.349685 cm-1     1.407182 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000    -0.001819    0.650441   -0.000575  
      0.000000  0.000000  0.750000    -0.002059    0.759551   -0.000620  

Finite differences POTIM=  1.500000000000000E-002
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2011-02-23 11:04:25
已阅   回复此楼   关注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的回帖
查看全部 25 个回答

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+1): 谢谢提示 2011-02-23 11:17:28
应该是第一个吧?实验H2振动频率是多少呢?

从数据来看, 你的H2是沿着Z轴放置的,沿键长方向振动是主要的.

虚频看来是分子整体振动造成的.

但是为什么会出现虚频呢?
赞一下(2人) 回复此楼
2楼2011-02-23 11:16:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢提示 2011-02-23 12:06:00
引用回帖:
Originally posted by cenwanglai at 2011-02-23 03:16:34:
应该是第一个吧?实验H2振动频率是多少呢?

从数据来看, 你的H2是沿着Z轴放置的,沿键长方向振动是主要的.

虚频看来是分子整体振动造成的.

但是为什么会出现虚频呢?

把EDIFFG设的更高，如1e-7，同时POTIM设的更小，如0.002，虚频可能就没有了。可以试试看。
赞一下(2人) 回复此楼
3楼2011-02-23 11:55:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ellsaking(金币+2): 谢谢指导 2011-02-23 15:11:58
引用回帖:
Originally posted by lgxyz at 2011-02-23 03:55:34:
把EDIFFG设的更高，如1e-7，同时POTIM设的更小，如0.002，虚频可能就没有了。可以试试看。

做频率计算的时候，H2可以沿着轴放置。表面吸附分子方向一般不是这么”规则“的，是不是要求矢量和呢？

—————————————————————————————————————

自己算了一个例子，实际上无论分子如何放置，频率计算是根据分子本身的对称性来做振动分析的，跟相对坐标无关。

[ Last edited by cenwanglai on 2011-2-23 at 14:27 ]
赞一下(2人) 回复此楼
4楼2011-02-23 12:21:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭