计算化学公社

标题: 求助：用合金吸附模型计算拉曼谱POLAR报错 [打印本页]

作者
Author: wuwei1 时间: 2024-6-21 09:30
标题: 求助：用合金吸附模型计算拉曼谱POLAR报错
各位老师，我用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 *
*******************************************************************************

作者
Author: wuwei1 时间: 2024-6-21 16:41
我现在打算用CP2K重新优化一下，还不行的话就没辙了

作者
Author: cpeng11 时间: 2024-8-24 14:04
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.

作者
Author: zjxitcc 时间: 2024-8-25 17:21
本帖最后由 zjxitcc 于 2024-8-25 21:11 编辑

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

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

作者
Author: wuwei1 时间: 2024-12-29 10:02

zjxitcc 发表于 2024-8-25 17:21
你可以尝试逐渐降低电子温度，看是否可以降到0 K。比如300K的计算读取600K下收敛的轨道。

好的，我尝试一下

作者
Author: wuwei1 时间: 2025-2-21 16:31
更新一下，这个问题我目前就只能不开smear，加上降低SCF收敛标准，强行算

欢迎光临计算化学公社 (http://bbs.keinsci.com/)