在NWChem源码中添加动能泛函模块请教

空山无人

设计了一个新的 meta-GGA level 动能密度泛函，已在 x2dhf 和 PROFESS 中得到了双原子分子和周期性体系下泛函的表现，想使用 NWChem 进一步计算一些有限体系大分子，现在需要在源码中添加泛函实现，并将积分格点上的动能和密度拉普拉斯量可视化。自己尝试修改代码但是一直改不对，想请教一下诸位老师同学是否有在 NWChem 中做这方面的经验，或者其他添加密度泛函形式到量化软件中的经验，多谢！

wzkchem5

hebrewsnabla 发表于 2023-7-19 11:29
你的新泛函是（1）仅修改了部分顶层参数或已有泛函的混合，还是（2）修改了已有泛函的底层参数，还是（3） ...

楼主做的是orbital-free DFT里的动能泛函，支持任意体系orbital-free DFT的程序还是比较少的。不过如果一个程序支持从libxc里调动能泛函的话，最简单的方法应该就是在libxc里仿照已有的动能泛函创建一个新的动能泛函就行了，毕竟接口都是现成的，不太容易出错

sobereva

量化程序一般使用原子中心格点，没法容易做到“动能和密度拉普拉斯量可视化”，除非你是用着色的小点来表示
这俩函数一般可视化都是等值面图、填色平面图、等值线图，Multiwfn直接就能实现而不需要额外改程序，若要在里面添加新的实空间函数来可视化也非常容易，改function userfunc函数即可

hebrewsnabla

你的新泛函是（1）仅修改了部分顶层参数或已有泛函的混合，还是（2）修改了已有泛函的底层参数，还是（3）完全从头做起的新泛函呢？

情况（1）通常在任何软件中只要能自定义输入就能实现；情况（2）（3）需要去改泛函库。

nwchem有自带的泛函库（但是我也不了解），也可以使用libxc。如果改自带泛函库进展不顺利的话，也可以尝试去改libxc。至少，libxc的开发者很活跃，应该可以给你一些帮助。

另外，有什么一定要使用nwchem的理由么，这并不是一个很好魔改的程序。魔改还是首推pyscf，psi4. 基于这两个程序来做的话，基本也是要去改libxc。