计算化学公社

标题: 为什么要用gjf文件输入Multiwfn得到RESP电荷？ [打印本页]

作者
Author: 黄舒伟 时间: 2021-3-28 20:16
标题: 为什么要用gjf文件输入Multiwfn得到RESP电荷？
新人想请教一下，将gjf文件输入guassian不就出来了.gesp文件可用[size=1.2em]antechamebr
生成resp电荷么？为什么还要用Multiwfn来输入.chk文件去得到RESP呢？是计算更快还是更准确呢？

还有想请教一下我这个gjf文件哪错了？guassian16运行不了：
%nprocshared=2
%mem=1GB
%chk=HP2.chk
# opt=tight hf/6-31g(d) nosymm scf=qc

HP2

0 1
原子略

HP2_ini.gesp

HP2.gesp

（两行空格）

提前多谢大佬回复

作者
Author: sobereva 时间: 2021-3-28 22:43
仔细看此文了解Multiwfn算RESP电荷的便利和强大
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

用antechamber你还得装Ambertools，编译就够费时间的而且可能中途出错，而且对于antechamber处理不了的杂七杂八的体系你根本没法通过antechamber获得RESP电荷（antechamber一般只能搞成键方式普通的有机体系）。相反，Multiwfn对任何体系都是普适的，有机/无机/分子/团簇全都能用。Multiwfn算RESP电荷还有极强的灵活性和可控性，片段电荷约束、等价性约束、多构象权重都可以非常方便地定义。Multiwfn还用起来特方便，照着屏幕提示敲几下键盘即可，都不用记忆任何命令、查看帮助。而且antechamber只能基于Gaussian算的静电势得到RESP电荷，而几乎所有主流量子化学程序的用户都可以靠Multiwfn算RESP电荷。

Multiwfn还可以一行命令直接得到RESP电荷，做法真是没法更简单：
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）

而且Multiwfn还能一行命令算出来比RESP明显更好的RESP2(0.5)电荷：
RESP2原子电荷的思想以及在Multiwfn中的计算
http://sobereva.com/531（http://bbs.keinsci.com/thread-16190-1-1.html）

antechamber直接给你的关键词太垃圾了，scf=qc哪能随便用，仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61
opt带着tight完全多余。用HF/6-31G*优化在如今来看根本无法接受，在HF/6-31G*下获得RESP电荷远不如用像样的级别算RESP2(0.5)电荷。写nosymm根本没意义。

记清楚Gaussian的拼写，老拼错是丢人的事

作者
Author: 黄舒伟 时间: 2021-3-28 23:22
好的，谢谢社长指教。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)