前言：

对于比较大的分子体系，我们有可能会为了减小计算量，在优化结构的时候，固定某些原子。这样会带来导致计算频率的时候，产生一些虚频；有时候即使优化的很好，也会得到一些波数很小、振幅也很小的频率，其中有的波数很小的频率其实是平动，而虚频则是负数，我们计算Franck-Condon Factors（FCF）的时候，并不希望把这些频率考虑进去。

解决方法：

ADF计算FCF的时候，会自动去掉波数（绝对值）很小的振动模式，但并不会自动去掉虚频。那么我们要去掉虚频以及很小的频率的时候，就可以通过MODES关键字来指定需要计算的频率范围。以下为一个例子：

我们以甲醛为例，分别计算了两种状态下的甲醛分子的振动。分别得到12个振动模式，

FCF的*.run文件如下：

$ADFBIN/fcf << eor
STATES D:\ADF_DATA\w1.t21 D:\ADF_DATA\w2.t21
QUANTA 0 0
MODES 8 10
eor

注意此处QUANTA 0 0表示计算涉及振动模式都只考虑振动基态，实际计算的时候，不能设置为0。

那么我们看FCF计算的*.out文件：

那么事实上，对于两个状态的甲醛分子，都只纳入了第2、3、4个振动模式。那么我们现在详细说明MODES的后面的数字的含义：

MODES M N表示计算两个状态的分子分别从第M-6到第N-6个振动模式。

如果分子中的原子个数为2或者3，则表示计算两个状态的分子分别从第M-3到第N-3个振动模式。原子个数为1，实际上就没有振动谱，也谈不上FCF。

6和3，实际上表示这个分子具有的6个或者3个平移运动，波数实际上为0。

ADF软件提供免费试用，试用申请方式参见费米科技维基百科：AMS免费试用