24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hzfish

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1515.8
散金: 10
红花: 3
帖子: 604
在线: 358.1小时
虫号: 753222
注册: 2009-04-20
性别: GG
专业: 理论和计算化学

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

作频率计算收敛标准出现以下情况:
      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??

谢谢!

回复此楼

» 猜你喜欢

PbS量子点如何确定浓度已经有11人回复
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）已经有22人回复
物理化学论文润色/翻译怎么收费? 已经有247人回复
PbS量子点紫外吸收已经有0人回复
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已经有2人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship 已经有0人回复
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料已经有0人回复
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建已经有0人回复
半透性容器进行高湿或者低湿研究已经有0人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生已经有11人回复
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事已经有0人回复

1楼 2009-12-04 19:01:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

niuniuzi41

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 2192.2
帖子: 112
在线: 33.3小时
虫号: 917089
注册: 2009-11-30
专业: 有机合成

消除虚频的技巧

★ ★ ★ ★
小木虫(金币+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,一般可以把多余的虚频去掉，再放开做全优化。或许会有所收获。