Multiwfn将fch文件转换成mkl文件失败

Oliviaw

我用的是multiwfn 3.8 dev的Windows版。我用formchk将chk转换成fch后，导入multiwfn，按照100-2-9-y依次输入，最后multiwfn闪退，没有生成mkl文件，可能是哪里出错了吗？


fch有2.36G

sobereva

Multiwfn FAQ
http://sobereva.com/452

搜 闪退

知道得能知道操作系统的报错提示才能判断
如果很类似计算条件算的小体系都能正常处理，大概率是可用内存不够。找个内存大的Linux机子转换再试

Oliviaw

sobereva 发表于 2024-3-30 16:01
Multiwfn FAQ
http://sobereva.com/452

Windows改了设置也不行，只有linux可以

Oliviaw

sobereva 发表于 2024-3-30 16:01
Multiwfn FAQ
http://sobereva.com/452

gaussian波函数作为orca初猜函数的帖子里说到非限制性波函数，所以写UKS关键词。ORCA的B3LYP和Gaussian的定义不同，所以在inp里面加/G。
请问怎么知道波函数是不是非限制性的？
请问含有过渡金属的中性体系，自选多重度是1，得到的波函数，是限制性还是是非限制性的?
ORCA的PW6B95-D3和gaussian的PW6B95D3定义相同吗？需要改变输入文件的关键词吗？
如果我orca直接用gaussian的波函数，不改关键词，影响大吗？

sobereva

Oliviaw 发表于 2024-3-30 20:21
gaussian波函数作为orca初猜函数的帖子里说到非限制性波函数，所以写UKS关键词。ORCA的B3LYP和Gaussian的 ...

你用UKS关键词（或者自旋多重度>1的情况下默认用了）自然就是非限制性的

PW6B95-D3没有什么使用价值，尤其对于ORCA用户来说。自己试一下便知

Oliviaw

sobereva 发表于 2024-3-30 20:33
你用UKS关键词（或者自旋多重度>1的情况下默认用了）自然就是非限制性的

PW6B95-D3没有什么使用价值， ...

我是不知道要不要用USK，还没有算。
您的意思是自旋多重度是1的体系就不是默认用USK，只有>1才是默认吗？
是不是不管我自旋多重度是多少，我只要用USK，轨道就会变成非限制性的呢？您是这个意思吗？

我是用的PW6B95-D3/def2-tzvp//TPSS-D3/def2-tzvp这样的组和，并且加上ATM校正和mRRHO校正和gcp。这个组合是别人benchmark推荐的用法。我对比了这个组合方法算的结合能的数值和其他文献的用M06-2X/6-311+G(d,p)//M06-2X/6-31G(d)算的结合能的数值只差0点几个kcal/mol。所以才确定用PW6B95来搭配TPSS的。

zjxitcc

Oliviaw 发表于 2024-3-30 20:21
gaussian波函数作为orca初猜函数的帖子里说到非限制性波函数，所以写UKS关键词。ORCA的B3LYP和Gaussian的 ...

不需要纠结ORCA的PW6B95-D3与gaussian的PW6B95D3定义，只要你计算过程合理正确，关键词合理，基组一致，那传轨道后ORCA自然会在5圈内收敛，你又不是非要让它1圈收敛。推荐尝试MOKIT的fch2mkl小程序。另外，D3不影响轨道，所以Gaussian计算是否带D3都无所谓，轨道系数是不变的。你想用D3，在ORCA里写D3即可。

“不改关键词”这话就令人迷惑了，关键词是自己写的，怎么叫“不改”，难道谁已经帮你写好了？如果关键词合适，自然不需要改；如果不合适，自然需要改。

sobereva

Oliviaw 发表于 2024-3-30 21:12
我是不知道要不要用USK，还没有算。
您的意思是自旋多重度是1的体系就不是默认用USK，只有>1才是默认吗 ...

那叫UKS
只要你写了UKS，哪怕是单重态、哪怕体系稳定波函数是闭壳层的，得到的波函数的形式上也是非限制性的形式，即alpha和beta两套轨道单独给出。

[size=14.6667px]M06-2X/6-311+G(d,p)算能量本来就不多准确，拿它当参照没意义。对于热力学数据的计算，wB97M-V几乎算什么类型的，平均精度都比PW6B95-D3和M06-2X强。