关于分子占据数的DFT计算并在此基础上计算拉曼光谱

aisanahtuE

之前有幸拜读了sob老师《使用NWChem做分数占据数的DFT计算》这篇文章，发现可以应用到最近自己在做的关于拉曼的研究里。
因为一直用的都是gaussian，nwchem是完全的初学者，对于nwchem的语法什么的也不是很了解。
所以输入文件直接复制了文章里的代码，只是改动了分子坐标以及添上了计算拉曼的部分。
运行之后输出文件一直停留在下边的这个部分。
请问该怎样写输入文件才能解决这个问题，或者是说分子占据数的计算有什么无法计算振动频率的限制么？
烦请各位老师不吝赐教！




Multipole analysis of the density
     ---------------------------------

     L   x y z        total         alpha         beta         nuclear
     -   - - -        -----         -----         ----         -------
     0   0 0 0     -0.500000    -55.500000    -55.000000    110.000000

     1   1 0 0     -0.002027     -0.006197      0.004170      0.000000
     1   0 1 0      1.099773      0.238323      0.861450     -0.000000
     1   0 0 1     -0.000000     -0.000000     -0.000000      0.000000

     2   2 0 0    -78.175699   -633.510616   -626.842246   1182.177163
     2   1 1 0     -0.144283     -0.308393     -0.322930      0.487040
     2   1 0 1      0.000000      0.000000      0.000000      0.000000
     2   0 2 0    -71.337921   -902.113018   -895.195039   1725.970136
     2   0 1 1      0.000000      0.000000      0.000000      0.000000
     2   0 0 2    -69.514840    -38.015287    -37.086142      5.586589

Parallel integral file used    3446 records with       0 large values

HESSIAN: the one electron contributions are done in      24.5s

sobereva

FON算拉曼貌似不行。这个可以到nwchem的论坛上向开发者进一步确认一下