分子的振动谱(IR, Raman)可以通过计算Γ-点声子谱获得，除了SIESTA主程序外还需要另外两个工具程序：fcbuild和vibrator，前者用于构建超胞（计算力常数的输入文件）后者用于对角化力常数矩阵，获得声子谱，其中力常数矩阵的计算仍然由siesta完成，这一步很耗时，并且对结构优化精度要求较高，否则会有非物理的虚频出现。
具体细节请参见 /SIESTA/Util/Vibra/Docs中的说明文件，工具程序在/SIESTA/Util/Vibra/Src文件夹中，编译前最好将vibra.h中的参数maxa（maximum number of atoms in the unit cell）设置大一些，至少大于你所计算的体系中原子个数。
祝你成功。
回复此楼
huangll99
+1/89
