CP2K2025.1版本涉及钠离子的能量似乎存在bug

ljh123

同一个轨迹，完全相同的计算级别和参数设置。CP2K2025用retraj计算得到的能量和CP2K2024、2023差距巨大（差不多0.6 a.u.)。看了下输出文件，似乎只有对于钠离子的初猜有不同。
我还尝试过包含C,H,O,Cl,S的体系，和CP2K之前版本得到的值基本一致。
目前来看是钠离子有问题？各位有没有注意过这个情况呢？感觉和能量相关的计算都要注意一下了（比如AIMD，

################################################################################
v2025.1:
Guess for atomic kind: Na

Electronic structure
    Total number of core electrons                                          2.00
    Total number of valence electrons                                       9.00
    Total number of electrons                                              11.00
    Multiplicity                                                   not specified
    S   [  2.00] 2.00  1.00
    P      6.00


*******************************************************************************
                  Iteration          Convergence                     Energy [au]
*******************************************************************************
                          1        0.540290                     -47.490094127456
                          2        0.896650E-01                 -47.533798426511
                          3        0.890137E-04                 -47.535723644572
                          4        0.144851E-04                 -47.535723653965
                          5        0.858547E-06                 -47.535723654019

#########################################################################
v2024.1以及2023.1:
Guess for atomic kind: Na


Electronic structure
    Total number of core electrons                                          2.00
    Total number of valence electrons                                       9.00
    Total number of electrons                                              11.00
    Multiplicity                                                   not specified
    S   [  2.00] 2.00  1.00
    P      6.00




*******************************************************************************
                  Iteration          Convergence                     Energy [au]
*******************************************************************************
                          1        0.708123                     -46.859771414279
                          2        0.131650                     -46.927516281371
                          3        0.147329E-03                 -46.931632581654
                          4        0.318297E-04                 -46.931632628122
                          5        0.657977E-06                 -46.931632628352

Uus/pMeC6H4-/キ

一年前的Commit 3204c26中，平时默认q9的钠的赝势如GTH-PBE有所更改。所以如果要比较2024.2及之后版本与2024.1及之前版本涉及钠的计算结果，应在新版本中明确采用如GTH-PBE-q9_old的旧款赝势。

ljh123

Uus/pMeC6H4-/キ 发表于 2025-4-2 09:44
一年前的Commit 3204c26中，平时默认q9的钠的赝势如GTH-PBE有所更改。所以如果要比较2024.2及之后版本与202 ...

原来如此，感谢回复。