详细内容

Gaussian 讲解66-Gaussian 09 Freq(4)

振-转耦合输出。如果指定了 VibRot 选项,则在输出文件的正常热化学分析输出之后, 立即打印简谐振转耦合的分析,开头的标题是:

image.png

如果同时还指定做非谐分析(即,同时指定了VibRotAnharmonic),则在简谐振转耦合分析之后,输出非谐振转分析的结果,开头的标题是:

image.png

非谐频率计算。在正常的频率输出之后,打印 Freq=Anharmonic 计算产生的附加输出。 (如果还指定进行振转耦合计算,则数据打印在振转耦合的输出之后)。这里我们简要讨论 其中最重要的项。  

输出中显示平衡结构(即,势能曲面的最低点),接着是0 K时的非谐振动平均结构:  

image.png

注意,后一结构的键长略长。非谐零点能在其各分量项之后的输出中立即给出:

image.png

非谐频率出现在这个表之后,位于标记为E(anharm)的列:

image.png

image.png

为了方便,还列出了谐振频率。  

检查低频。完整力常数矩阵对角化(默认的 Freq=DiagFull)的输出以如下形式出现, 其中保留了转动和平移自由度:

image.png

这个输出来自甲醇的 Opt Freq 计算。注意那些基本为零的模式,最低的模式(忽略符 号)位于大约 14,19,和 25 波数。回到计算中,如果用严格的优化标准和更多的积分网格 Opt=Tight Int=UltraFine),最低的模式变为:  

image.png

现在低频模式已经相当小了,而且最低的频率也有所移动。  

对于具有小波数频率的分子体系,这个分析特别重要。例如,如果打印出的最低频率位 于大约 30,而在其上大约 25 的位置有一个低频模式,那么前一个数值十分可疑(因为分子 结构甚至可能不在极小点)。  

用不同的热化学参数运行频率计算。下面的两步计算任务先进行频率计算,之后使用 不同温度、压强和选择的同位素进行第二个热化学分析:  

image.png

还要注意,freqchk 工具可以使用存在 Gaussian 检查点文件中的频率数据,重新运行 热化学分析。

FREQ=READANHARM 的附加输入 

这个输入从单独部分读入,可以包含以下关键词:  

Fermi       还执行各向同性超精细耦合的振动平均。 

PrintGeom    打印做振动平均特性计算的几何结构。

TolFre=x     Fermi 共振的最小频率差(cm -1)(默认值为 10.0)。必须是实数。 

DaDeMin=x    Darling-Dennison 共振的最小频率差(cm -1)(默认值为 10.0)。必须是实数。 

TolCor=x     Coriolis 耦合的阈值(cm -1)(默认值为 10-3)。必须是实数。 

ScHarm=x     谐振频率做线性换算的换算因子(B3LYP/6-31+G(d)为 1.0 x 10-5)。必须是 实数。默认使用来自 Scale 关键词的值。 

FREQ= READFCHT 的附加输入 

这个输入从单独部分读入,可以包含以下关键词:  

MaxOvr=N     在计算对应于跃迁到单激发振动态的 Franck-Condon 因子时,设定达到的最 大泛频。默认值为 20。 

MaxCMB=N     设定两个态达到的最大泛频,这两个态组合成最终态。默认值为 13。 

MaxInt=N     设定在每一类跃迁中计算的积分最大个数(单位:百万)。默认值为 100。 

NoIntAn      停用 Sharp 和 Rosenstock 分析公式,它用于计算到单个泛频和到双态组合的 跃迁积分。 

NoRelI00      给定的光谱范围默认相对于 I00跃迁的能量。如果由用户给定光谱边界的绝对 能量,必须给出这个关键词。

SpecMin=x     设定最终光电子谱的下限(单位 cm-1)。必须是实数。默认值为-1000。

SpecMax=x     设定最终光电子谱的上限(单位 cm-1)。必须是实数。默认值为+8000。  

SpecRes=x     设定离散谱的两点之间间隙(单位 cm-1)。这个值会极大地影响计算时间, 非常低的值会极度减慢计算,特别是当 HWHM 设得很高时。必须是实数。默认值为 8。  

SpecHwHm=x    设定用高斯函数表示的谱带半高处的半宽度(单位 cm-1)。必须是实数。默 认值为 135。  

DeltaSP=x     设定由于收敛差而终止计算的阈值。这个值应当小于 1.0(对应于非常好的 收敛)。默认为 0.0(不结束计算)。  

AllSpectra    在 Gaussian 的输出中,除了最终光谱还打印每一套组合(类型)得到的光谱。 默认停用这个打印。  

PrtMat=N     一串数字,打印不同的矩阵用于积分计算的基:1 表示 Duschinsky 矩阵 J2 表示移动矢量 K3 表示 A4 表示 B5 表示 C6 表示 D7 表示 E,其中的 ABCDE 是 Sharp 和 Rosenstock 的矩阵。数字的顺序不重要。默认值为。  

PrtInt=x     设定在输出中打印哪些积分。阈值单位是 I00强度。必须是实数。默认值为 0.01。  

DoTemp      在光谱计算中包含温度。默认在 0 K 执行光谱计算。  

MinPop=x     设定做为跃迁起始点的振动态的最小布居。默认值为 0.1。  

InFrS0       对于初态,强制程序使用用户提供的频率。这个频率在 Freq=ReadFCHT 选项 行之后的输入中指定。  

InFrS1       对于终态,强制程序使用用户提供的频率。这个频率在 Freq=ReadFCHT 选项 行之后的输入中指定。  

JDusch, JIdent  强制程序使用标准的 Duschinsky 矩阵(默认为 JDusch),或把单位矩阵做 为 Duschinsky 矩阵(JIdent)。对于后一种情况,不考虑模式的转动。默认值为 0。

SclVec        进行换算矢量的计算,用初态频率的换算矢量和 Duschinsky 矩阵修改终态的 频率。如果给定这个关键词,会要求用户提供初态的频率,方式同 InFrS0。 

EnerInp=x      使用由用户提供的初态和终态之间的 ΔE,替换计算的值。必须是实数。默认 值为 0。

最新评论
请先登录才能进行回复登录
技术支持: CLOUD | 管理登录
seo seo