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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【讨论】频率计算收敛标准！！！！
91/1返回列表
查看: 1092  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

hzfish

金虫 (正式写手)

[交流] 【讨论】频率计算收敛标准！！！！

作频率计算收敛标准出现以下情况:
         Item               Value     Threshold  Converged?
Maximum Force            0.000022     0.000450     YES
RMS     Force            0.000005     0.000300     YES
Maximum Displacement     0.015024     0.001800     NO
RMS     Displacement     0.002646     0.001200     NO
Predicted change in Energy=-3.207076D-08

出现这种情况结果是否可以采用.采用什么方法使结果出现四个YES??

谢谢!
回复此楼

» 猜你喜欢
1楼2009-12-04 19:01:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agou

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我也遇到了这样的问题。期待
赞一下(1人) 回复此楼
小木虫三戒总理事
2楼2009-12-04 19:27:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3867826

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
可以啊  
  2个YES 可以默认收敛了！
赞一下(1人) 回复此楼
拒稿，终于体会到了！乱写接受，认真写拒稿，这就是2和4的区别吗？
3楼2009-12-04 19:53:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NUPT

新虫 (小有名气)
yjcmwgk(金币+0,VIP+0):事实上，在freq计算时，两个force值必须yes。这一点必须搞清楚，我们不能把四个yes混为一谈，要把他们分开来考虑。原因我在下面解释出来了。 12-5 16:44
引用回帖:
Originally posted by 3867826 at 2009-12-4 19:53:
可以啊  
  2个YES 可以默认收敛了！

为什么一定要4个都是yes？  2个就够了。。。
赞一下(10人) 回复此楼
4楼2009-12-04 20:14:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1,VIP+0):事实上，在freq计算时，两个force值必须yes。这一点必须搞清楚，我们不能把四个yes混为一谈，要把他们分开来考虑。原因我在下面解释出来了。 12-5 10:43
如果没有虚频，即使没有四个yes也可以认为是收敛的
赞一下(11人) 回复此楼
5楼2009-12-04 20:48:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+5,VIP+0):感谢交流 12-5 11:20
引用回帖:
Originally posted by zpppanda111 at 2009-12-4 20:48:
如果没有虚频，即使没有四个yes也可以认为是收敛的

简单解释：
楼主要搞清一点。freq的计算与Displacement无关，freq的计算仅仅与Force有关。既然Force已经收敛，你管Displacement干什么？
事实上，在freq计算时，两个force值必须yes。这一点必须搞清楚，我们不能把四个收敛标准混为一谈，要把他们分开来考虑。因为freq计算力常数的二阶导数，因此，freq计算时，两个force收敛值必须yes

稍微详细的解释（解释一下“为什么freq计算与Displacement无关”）：
首先，要搞懂什么是频率。其实，只需要中学知识我们就能理解这个问题：
中学的时候我们学过简谐振动，对应的回复力是f=-kx(k>0),对应的能量曲线，是一个开口向上的二次函数E=(1/2)k(x^2). 这样的振动，对应的x=0的点是能量极小值点（简单情况下也就是最小值点），此时二阶导数E"=k>0。这时的振动频率我们也会求：ω=(2π)[sqrt(k/m)]。显然它是一个正的频率，也就是通常意义下的振动频率。
那么，一维情况下，如果能量曲线是一个开口向下的二次曲线呢？首先，从能量上看，这是个不稳定的点，中学的物理书上称为“不稳平衡”。用现在的观点看，就是这一点一阶导数是零（受力为0），且是能量极大值。E"=a<0.如果套用上面的公式，“回复力”f=-ax（实际上已经不是回复，而是让x越来越远了），这里a<0，ω=(2π)[sqrt(a/m)]显然就是一个虚数(2π)[sqrt(-a/m)]i，即所谓的虚频。Gaussian里面给出一个负的频率，就是对应这个虚频的。

更详细的解释，自己去啃书吧，哈哈，好几章呢

[ Last edited by yjcmwgk on 2009-12-6 at 12:38 ]
赞一下(11人) 回复此楼
6楼2009-12-05 10:41:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by yjcmwgk at 2009-12-5 10:41:

简单解释：
楼主要搞清一点。freq的计算与Displacement无关，freq的计算仅仅与Force有关。既然Force已经收敛，你管Displacement干什么？
事实上，在freq计算时，两个force值必须yes。这一点必须搞清楚，我们不 ...

学习到了
膜拜
回复此楼
7楼2009-12-05 11:03:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
very good，yjcmwgk 版主的解释很清楚，很详细
赞一下(1人) 回复此楼
8楼2009-12-05 16:04:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuniuzi41

金虫 (小有名气)

消除虚频的技巧

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 12-5 20:55
yjcmwgk(金币+0,VIP+0):楼主的问题和消虚没关系 12-6 12:18
loovfnd(金币+0,VIP+0):新虫，混金币可以理解。下次不可以 12-8 13:39
在进行Gaussian计算时，经常会出现多个虚频个，每个都只有－30~－20左右，但是如果结构改动过大，则一不小心，虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试：nosymm；加大循环次数；提高收敛度；iop(1/8=1)等。通常对于势能面很平，较小的虚频很难消除。   为了解决以上的这种现象,按如下的步骤一步步调整：
  1）. 在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。
  2）. 对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。
  关于如何将内坐标改成直角坐标，可以直接看输出文件（例子如下）：
  Full point group CS NOp 2
  Largest Abelian subgroup CS NOp 2
  Largest concise Abelian subgroup CS NOp 2
   Standard orientation:
  ---------------------------------------------------------------------
  Center Atomic Atomic Coordinates (Angstroms)
  Number Number Type X Y Z
  ---------------------------------------------------------------------
   1 46 0 -0.242506 1.028743 2.204311
   2 46 0 -0.476137 2.619863 0.000000
   3 46 0 -0.242506 1.028743 -2.204311
   4 46 0 -0.242506 -1.618028 -1.382208
   5 46 0 -0.242506 -1.618028 1.382208
   6 46 0 -1.725240 0.349678 0.000000
  Coordinates这一部分就是，如果没有，看Input orientation也行。
  3）. 如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到虚频没有。比较好的做法是分子结构输入用直角坐标，将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加，可以选择权重如1/3相加等。
  4）. 实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的 东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做一般的电子结构计算，这样，应当没有问题。
  如果上面几步仍然无法解决问题，则有可能是因为计算还没有得到最稳定的结构，那么或者是分子结构有问题，或者是计算错了，再者就是游离出现在计算的范畴。
  最后如果是寻找过渡态，如果只有一个虚频，而且该虚频对应某单键旋转震动，那么产生虚频的原因就是初始结构比较接近重叠构象。这时把优化过的结构绕振动方向旋转几十度再重新优化就有可能得到正确的结果。优化中等大小的有机分子（50-100原子）的时候，应该经常碰到这种情况。
  当出现两个以上虚频时，可以先看一下振动模式。如果多余的虚频，和想要的振动模式发生在不同原子上，可以先固定过渡态模式所对应的原子，然后做opt到minimal,一般可以把多余的虚频去掉，再放开做全优化。或许会有所收获。
赞一下(31人) 回复此楼
9楼2009-12-05 19:31:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hzfish 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭