求助：用合金吸附模型计算拉曼谱POLAR报错

wuwei1

各位老师，我用CP2K来模拟合金吸附模型的拉曼谱，我先把vasp里面优化好的cif拉入Multiwfn，输出CP2K的输入文件。没有切换OT，开启smearing，固定属于合金的原子，输入和输出文件如下。但算到POLAR这一步出现报错了，也对比了Sob老师在“http://bbs.keinsci.com/thread-26800-1-1.html”中回复的输入文件，我属于菜鸟看不出啥毛病。。。请求各位大佬帮忙看一下
LINRES| Properties to be calculated:
                                                        Electric Polarizability
LINRES| Optimization algorithm                             Conjugate Gradients
LINRES| Preconditioner                                                FULL_ALL
LINRES| EPS                                                            1.0E-06
LINRES| MAX_ITER                                                            50


  Total energy ground state:                              -11249.49509322690210


POLAR| Starting polarizability calculation
POLAR| Initialization of the polar environment
POLAR| Self consistent optimization of the response wavefunctions




POLAR| Response to the perturbation operator Berry phase_x


  Iteration    Method   Restart      Stepsize      Convergence         Time
  ------------------------------------------------------------------------------


*******************************************************************************
*   ___                                                                       *
*  /   \                                                                      *
* [ABORT]                                                                     *
*  \___/                             CPASSERT failed                          *
*    |                                                                        *
*  O/|                                                                        *
* /| |                                                                        *
* / \                                                 qs_linres_methods.F:730 *
*******************************************************************************

wuwei1

我现在打算用CP2K重新优化一下，还不行的话就没辙了

cpeng11

https://groups.google.com/g/cp2k/c/qOCvYy3TRTE/m/wWtq291eBAAJ
来自Jürg Hutter的回答。smearing开了没法算。请问你有找到解决办法吗？

this error is triggered when two arrays have different dimensions.
Without more information, I can only guess that you are trying to
get dipole moments for a periodic system while using smearing of
electronic states. This combination of options is not possible.

zjxitcc

cpeng11 发表于 2024-8-24 14:04
https://groups.google.com/g/cp2k/c/qOCvYy3TRTE/m/wWtq291eBAAJ
来自Jürg Hutter的回答。smearing开了 ...

你可以尝试逐渐降低电子温度，看是否可以降到0 K。比如300K的计算读取600K下收敛的轨道。