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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于频率验证
131/2返回列表12下一页
查看: 978  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yuyaolun

木虫 (小有名气)

[交流] 【求助】关于频率验证

大家好！
我优化了甲酸分子的结构，想算一下频率，可是一直有三到四个虚频，请问是怎么回事？谢谢！

我用的INCAR:
PREC    = Accurate  
ENCUT   = 400      
LREAL   = AUTO      
ISMEAR  = 0         

EDIFF   = 1E-6
EDIFFG  = -0.01     
NSW     = 100      
                    
IBRION  = 5         
POTIM   = 0.04     
NFREE   = 2         
ISYM    = 0         
ALGO    = Fast
IWAVPR  = 1      
LORBIT  = 10     
LWAVE   = .FALSE.
LCHARG  = .FALSE.
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
allen
1楼2009-12-17 14:28:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★
TIGERYZZ(金币+2,VIP+0):谢谢专家指点 12-17 15:47
虚频是必然了
虽然平移模和旋转模理论上=0
但数值计算总是有误差的
但虚的太多就有问题了；）
赞一下(6人) 回复此楼
2楼2009-12-17 14:44:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★
TIGERYZZ(金币+1,VIP+0):谢谢交流 12-17 15:47
同意二楼，如果是过渡态，出现虚频就正常的！但是如果是稳定态，最多出现三个左右的虚频，如果太多就不正常了！！呵呵！
赞一下(6人) 回复此楼
3楼2009-12-17 15:00:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyaolun

木虫 (小有名气)
同问前面两位高手：

稳定态的频率验证不应该全是实数吗？
过渡态的话，也应该只出现一个虚频吧……？
我不理解怎么出现三个虚频还有可能是正常的呢？因为我做一些分子在金属面上吸附的情况，频率验证是可以全是实数的。当然也有一些是有虚频的，但我都认为那些构型是不会存在的。
赞一下 回复此楼
allen
4楼2009-12-17 15:18:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
★
ice_rain(金币+1,VIP+0):谢谢指点 12-17 17:10
看你虚频波数是多少，很小的话可以忽略。
赞一下(6人) 回复此楼
5楼2009-12-17 16:11:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
1 f  =   10.092273 THz    63.411620 2PiTHz  336.641971 cm-1    41.738300 meV
   2 f  =   10.038364 THz    63.072899 2PiTHz  334.843758 cm-1    41.515350 meV
   3 f  =    9.505646 THz    59.725735 2PiTHz  317.074212 cm-1    39.312206 meV
   4 f  =    8.935139 THz    56.141133 2PiTHz  298.044142 cm-1    36.952777 meV
   5 f  =    8.920263 THz    56.047668 2PiTHz  297.547948 cm-1    36.891257 meV
   6 f  =    8.880352 THz    55.796897 2PiTHz  296.216647 cm-1    36.726197 meV
   7 f  =    8.778780 THz    55.158704 2PiTHz  292.828584 cm-1    36.306130 meV
   8 f  =    8.734126 THz    54.878134 2PiTHz  291.339082 cm-1    36.121456 meV
   9 f  =    8.638383 THz    54.276559 2PiTHz  288.145419 cm-1    35.725492 meV
  10 f  =    8.594012 THz    53.997770 2PiTHz  286.665374 cm-1    35.541989 meV
  11 f  =    8.410831 THz    52.846809 2PiTHz  280.555108 cm-1    34.784413 meV
  12 f  =    8.075051 THz    50.737042 2PiTHz  269.354699 cm-1    33.395738 meV
  13 f  =    8.058464 THz    50.632820 2PiTHz  268.801404 cm-1    33.327139 meV
  14 f  =    8.039531 THz    50.513864 2PiTHz  268.169887 cm-1    33.248840 meV
  15 f  =    7.962229 THz    50.028161 2PiTHz  265.591369 cm-1    32.929145 meV
  16 f  =    7.645244 THz    48.036486 2PiTHz  255.017890 cm-1    31.618200 meV
  17 f  =    7.601196 THz    47.759725 2PiTHz  253.548612 cm-1    31.436033 meV
  18 f  =    7.570667 THz    47.567905 2PiTHz  252.530267 cm-1    31.309774 meV
  19 f  =    5.296875 THz    33.281246 2PiTHz  176.684718 cm-1    21.906121 meV
  20 f  =    5.284840 THz    33.205626 2PiTHz  176.283267 cm-1    21.856347 meV
  21 f  =    5.155353 THz    32.392041 2PiTHz  171.964076 cm-1    21.320836 meV
  22 f  =    5.149187 THz    32.353294 2PiTHz  171.758371 cm-1    21.295332 meV
  23 f  =    5.039534 THz    31.664328 2PiTHz  168.100766 cm-1    20.841846 meV
  24 f  =    5.025325 THz    31.575047 2PiTHz  167.626786 cm-1    20.783080 meV
  25 f  =    4.778604 THz    30.024857 2PiTHz  159.397081 cm-1    19.762726 meV
  26 f  =    4.778604 THz    30.024852 2PiTHz  159.397055 cm-1    19.762723 meV
  27 f  =    4.630986 THz    29.097340 2PiTHz  154.473046 cm-1    19.152224 meV
  28 f  =    4.630976 THz    29.097281 2PiTHz  154.472732 cm-1    19.152185 meV
  29 f  =    4.629826 THz    29.090054 2PiTHz  154.434361 cm-1    19.147427 meV
  30 f  =    4.564253 THz    28.678047 2PiTHz  152.247084 cm-1    18.876239 meV
  31 f  =    4.564236 THz    28.677940 2PiTHz  152.246519 cm-1    18.876169 meV
  32 f  =    4.501913 THz    28.286355 2PiTHz  150.167659 cm-1    18.618424 meV
  33 f  =    4.393300 THz    27.603918 2PiTHz  146.544710 cm-1    18.169235 meV
  34 f/i=    0.019003 THz     0.119399 2PiTHz    0.633871 cm-1     0.078590 meV
  35 f/i=    0.026485 THz     0.166411 2PiTHz    0.883450 cm-1     0.109534 meV
  36 f/i=    0.026954 THz     0.169359 2PiTHz    0.899099 cm-1     0.111474 meV
赞一下 回复此楼
6楼2009-12-17 16:34:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★ ★ ★ ★ ★
aylayl08(金币+3,VIP+0):感谢热心交流 12-17 18:55
yuyaolun(金币+2,VIP+0): 12-17 20:04
上面这个是俺计算某个体系(稳定的)的结果，最后三个就是虚的！！！
出现这样的结果，正如二楼所说的，是数值计算的误差导致的！
但是这三个频率很小，直接可以忽略的！！

俺也用过VASP+ATAT结合的方法，得到过的gamma点的振动频率，总有三个很小(与其他的相比较）或正或负，所以直接忽略！

还有一个问题就是，有的结构能量上不是稳定的，但力学上是稳定的（比如某些亚稳相），那么得到的结果有可能没有虚频！！

至于过渡态理论上是不是只有一个虚频，我也不是十分肯定。我算过得扩散中的一些过渡态有四个很小的频率，或正或负。从扩散的角度出发，总有一个能量最低路径MEP，就对应着这个虚频！

以上拙见，不对的，请各位指正啊！！呵呵！！
赞一下(10人) 回复此楼
7楼2009-12-17 16:54:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyaolun

木虫 (小有名气)
这是我算的HCOOH结果

   1 f  =  108.202797 THz   679.858226 2PiTHz 3609.256704 cm-1   447.490960 meV   
   2 f  =   89.396732 THz   561.696235 2PiTHz 2981.953922 cm-1   369.715299 meV   
   3 f  =   52.973751 THz   332.843894 2PiTHz 1767.014079 cm-1   219.081903 meV   
   4 f  =   40.281139 THz   253.093862 2PiTHz 1343.634136 cm-1   166.589461 meV   
   5 f  =   37.465077 THz   235.400024 2PiTHz 1249.700427 cm-1   154.943161 meV   
   6 f  =   32.130510 THz   201.881946 2PiTHz 1071.758402 cm-1   132.881154 meV   
   7 f  =   29.949268 THz   188.176801 2PiTHz  999.000021 cm-1   123.860261 meV   
   8 f  =   20.437823 THz   128.414630 2PiTHz  681.732377 cm-1    84.524073 meV   
   9 f  =   17.842455 THz   112.107450 2PiTHz  595.160210 cm-1    73.790488 meV   
  10 f  =    1.768010 THz    11.108733 2PiTHz   58.974455 cm-1     7.311903 meV   
  11 f  =    0.541482 THz     3.402232 2PiTHz   18.061895 cm-1     2.239390 meV   
  12 f  =    0.032257 THz     0.202676 2PiTHz    1.075971 cm-1     0.133403 meV   
  13 f/i=    0.006576 THz     0.041318 2PiTHz    0.219351 cm-1     0.027196 meV   
  14 f/i=    0.153763 THz     0.966121 2PiTHz    5.128979 cm-1     0.635913 meV   
  15 f/i=    1.141286 THz     7.170909 2PiTHz   38.069191 cm-1     4.719980 meV

还有，你说虚频小就可以忽略，有什么判断标准吗？一般要小于多少才可以这样认为？

[ Last edited by yuyaolun on 2009-12-17 at 17:16 ]
赞一下 回复此楼
allen
8楼2009-12-17 17:12:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★ ★ ★ ★ ★
aylayl08(金币+3,VIP+0):谢谢专家解答 12-17 18:55
yuyaolun(金币+2,VIP+0): 12-17 20:03
yuyaolun(金币+1,VIP+0): 12-17 20:04
引用回帖:
Originally posted by yuyaolun at 2009-12-17 15:18:
同问前面两位高手：

稳定态的频率验证不应该全是实数吗？
过渡态的话，也应该只出现一个虚频吧……？
我不理解怎么出现三个虚频还有可能是正常的呢？因为我做一些分子在金属面上吸附的情况，频率验证是可以全 ...

理论上是0频的情况
数值计算可能会给出很小虚频
因此忽略他们

0频的数目不一定是3个
取决与你计算的系统
如果是晶体，自然是3个
如果是分子，应该是5个；）
赞一下(12人) 回复此楼
9楼2009-12-17 18:01:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★
aylayl08(金币+1,VIP+0):谢谢交流 12-17 18:56
引用回帖:
Originally posted by yuyaolun at 2009-12-17 17:12:
这是我算的HCOOH结果

   1 f  =  108.202797 THz   679.858226 2PiTHz 3609.256704 cm-1   447.490960 meV   
   2 f  =   89.396732 THz   561.696235 2PiTHz 2981.953922 cm-1   369.715299 meV   
   3 f ...

判断标准没有硬性规定
要自己分析
具体可以参照这里
http://blog.sina.com.cn/s/blog_5f15ead20100cyaq.html
赞一下(10人) 回复此楼
10楼2009-12-17 18:03:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuyaolun 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭