Uus/pMeC6H4-/キ 发表于 2025-1-19 18:25
硬要再说一个CP2K的不足的话，不用k点计算时记录波函数信息的二进制文件用了.wfn后缀应该算一个。不仅不方 ...

CP2K可以把各k点的轨道展开系数在输出文件里打印出来，可以由此在分析时利用，但是输出格式不够紧凑、浪费空间

硬要再说一个CP2K的不足的话，不用k点计算时记录波函数信息的二进制文件用了.wfn后缀应该算一个。不仅不方便提取分析（虽然参考博文http://sobereva.com/651用molden能弥补一下），还很容易和正宗的plaintext可读的.wfn文件搞混。

Google讨论区的这个帖子里，有人分享了叫cp2k_wfntools的python小工具，先用scipy.io.FortranFile(wfn_file, "r")读取文件再分析波函数信息。不过这个2018年的工具还能不能用于新版本的CP2K输出文件都难说，而且是否兼容考虑k点计算产生的.kp文件也不清楚。

学习了，真是深入细致。

sobereva 发表于 2024-12-16 02:21
还有CPMD、JDFTx、OpenMX、DFTB+、GPAW、FLEUR、Dacapo、ACE-Molecule、BigDFT
ABCUS应该是ABACUS

点赞，已经更改。卢老师厉害、牛！牛！牛！

gog 发表于 2024-12-15 22:51
免费开源的第一性原理理论计算包：QE、ELK、ABCUS、Abinit、Siesta。供参考

还有CPMD、JDFTx、OpenMX、DFTB+、GPAW、FLEUR、Dacapo、ACE-Molecule、BigDFT
ABCUS应该是ABACUS

免费开源的第一性原理理论计算包：QE、ELK、ABACUS、Abinit、Siesta。供参考

Stars 发表于 2024-12-9 09:21
VASP也没有解析Hessian吧？VASP官网上的内容是
“In general, the DFPT routines in VASP are somewhat ru ...

支持DFPT算q=gamma点的解析Hessian，其中仅部分利用有限差分

VASP也没有解析Hessian吧？VASP官网上的内容是
“In general, the DFPT routines in VASP are somewhat rudimentary and only support displacements commensurate with the supercell, i.e., so-called q=0 phonons. Therefore, this approach offers few advantages over computing phonons from finite differences. In particular, the DFPT routines are limited to LDA and GGA functionals, and they do not determine the elastic tensors, since the perturbation with respect to the strain tensor is not implemented. The only advantage of the linear response routines is that they eliminate the need to choose the magnitude of the finite displacement POTIM. Therefore, it might be helpful to first calculate phonon frequencies within DFPT and then switch to the finite differences approach in order to determine the largest displacement that will produce results compatible with the linear response routines.”
并且下文
“A few technical comments: VASP solves the linear Sternheimer equation to determine the linear response of the orbitals. Hence, unoccupied orbitals are not required. Internally, the VASP routines for linear response rely on finite differences in two places:

The first place is the determination of the second derivative of the exchange-correlation functional: Since most functionals do not support an algebraic determination of second derivatives, VASP always resorts to finite differences to determine the second-order change of the exchange correlation-potential and the PAW one-center terms for each atomic displacement.
Second, after VASP has determined the first-order change of the orbitals, it computes all second derivatives using finite displacements.”
表明这依然部分使用有限差分，更别提IBRION=5或6的情况了。

ma455173220 发表于 2024-11-30 07:00
CP2K还接有deepMD-kit Machine learning的接口，可惜没有wannier90的接口

其实cp2k应该存在所谓的“wannier90接口”，只是这个功能一直是实验性的，还没有公开的使用方法。

CP2K还接有deepMD-kit Machine learning的接口，可惜没有wannier90的接口

Uus/pMeC6H4-/キ 发表于 2024-11-29 19:05
CP2K虽然擅长AIMD，但又仿佛在努力控制情绪不要时而热情时而冷淡，另一边又害羞得在手册里没有多言。结果就 ...

我是很希望官方加上类似GROMACS那样退火设置的。在官方google讨论区提议的话有可能有开发者愿意加上，举手之劳的事，原理上自己魔改也不困难。若谁有时间魔改下在论坛里提供个patch能方便许多人。

牧生 发表于 2024-11-29 15:47
社长提到的这个问题，

需要额外指定，官网上有SIRIUS的用法说明

CP2K虽然擅长AIMD，但又仿佛在努力控制情绪不要时而热情时而冷淡，另一边又害羞得在手册里没有多言。结果就是分子动力学只有专一的目标平均温度和不断波动的瞬时温度，没有对稳定速率线性升温/降温的支持，还很容易引来“怎么像其他软件那样设置退火/变温”的日经问题...{ /) (\ }

http://bbs.keinsci.com/thread-39214-1-1.html

也是一年前的问题，关于用CP2K做自旋极化FPMD的，最后还是放弃了。