关于高斯读取电荷密度矩阵作为输入文件的方法的求教

req

感谢各位在百忙之中一看我遇到的疑惑。VASP可以通过参数调整来使用电荷密度矩阵作为数个输入文件之一。那么请问高斯是否也有类似的功能呢？查找IOp说明只找到了6/56这个选项，提到了"Whether to read in density matrix from input stream."。然而即使这个选项可以实现读取输入文件里电荷密度矩阵的功能，但我却并不知道电荷密度在高斯输入文件中的格式。
如果有哪位知道解答，望不吝赐教。先谢。(p.s. 有P站表情真是一大惊喜)

zjxitcc

你是指高斯里的孤立分子计算，还是周期性（晶体）计算？如果是孤立分子的计算，高斯的分子轨道保存在.chk和.fchk文件内，你换个gjf文件做输入，去guess=read它，自然就读进来了，高斯一般不使用读密度矩阵这种功能，读轨道进来它会自己构造密度矩阵。如果是周期性计算，没试过。

而且电荷密度矩阵这种说法用于第一性原理的软件，在高斯这些从头算软件里一般叫密度矩阵，不叫电荷密度矩阵。

req

非常感谢解惑。
是孤立分子计算。我对于回答的理解是原子坐标还是通过gjf提供，附加关键字guess=read和chk=某fchk来调用密度矩阵。这样就是替换掉第一步默认的初猜对吧。
关于名称的话，密度矩阵没记错的话就是电子出现概率的矩阵，和电荷密度矩阵应该是有共通点的吧（虽然不一定是一个东西）。

再次感谢你，我去尝试一下。具体情况是我想使用通过其他途径获得的电荷密度（来自cube文件）来作为高斯的输入的文件。所以我去看一看如何把cube里的mesh放入fchk文件里（如果可能的话）。

zjxitcc

req 发表于 2019-7-19 08:43
非常感谢解惑。
是孤立分子计算。我对于回答的理解是原子坐标还是通过gjf提供，附加关键字guess= ...

通过%chk=xxx.chk来实现。如果只有.fchk，需要事先转化为.chk文件。没错，这就是设置初猜。
感觉cube与SCF的密度矩阵有本质区别。。。不看好这能行···

sobereva

量化里没有“电荷密度矩阵”这个词，至少我从没听说过
fch文件根本就不记录电子密度格点数据，根本不可能拿cube里的信息当初猜

req

sobereva 发表于 2019-7-20 00:09
量化里没有“电荷密度矩阵”这个词，至少我从没听说过
fch文件根本就不记录电子密度格点数据，根本不可能 ...

谢谢sobe老师指点。

req

非常抱歉又重启话题。
请问如果只有cube，也就是原子坐标和电荷密度网格，有能够求得单点能的方法吗？
目前被堵死的路是修改fchk给高斯，目前在尝试的路是改写CHGCAR给VASP。想问有没有其他的可能方法。谢谢。

sobereva

req 发表于 2019-7-22 17:11
非常抱歉又重启话题。
请问如果只有cube，也就是原子坐标和电荷密度网格，有能够求得单点能的方法吗？
目 ...

没有。
除非已经有精确能量泛函，这样你把电子密度分布代入精确能量泛函就得到了精确的电子能量，然而众所周知，精确的能量泛函目前是没有的，谁提出来谁就能拿诺贝尔奖。
哪怕用近似的泛函，由于目前的动能泛函对于一般化学体系来说精度还很差，你也不可能得到定量准确的单点能

req

sobereva 发表于 2019-7-23 19:13
没有。
除非已经有精确能量泛函，这样你把电子密度分布代入精确能量泛函就得到了精确的电子能量，然而众 ...

非常感谢sobe老师的解惑。即使只是告知不可行也是很大的帮助。

req

sobereva 发表于 2019-7-23 19:13
没有。
除非已经有精确能量泛函，这样你把电子密度分布代入精确能量泛函就得到了精确的电子能量，然而众 ...

我又厚颜无耻的重开话题了。
即使是不精确的，近似的，得不到定量准确单点能的动能泛函，我也想尝试一下。所以敢问一下，有哪些候选的泛函可以参考，以及如果有可查阅的代码就帮助更大了。
先感谢大家。

sobereva

req 发表于 2019-8-22 15:14
我又厚颜无耻的重开话题了。
即使是不精确的，近似的，得不到定量准确单点能的动能泛函，我也 ...

我最近参与的一篇关于动能泛函对比的文章：J. Chem. Phys. 150, 204106 (2019); doi: 10.1063/1.5095072
Multiwfn能基于文中的二十多种动能泛函计算电子动能

req

sobereva 发表于 2019-8-23 22:33
我最近参与的一篇关于动能泛函对比的文章：J. Chem. Phys. 150, 204106 (2019); doi: 10.1063/1.5095072
...

感谢sobe老师。我先去拜读。