求助resp电荷的模拟

kjxwl3

各位大神，
     我在进行resp电荷的模拟的时候，在gaussian中有两种模拟方式：
--Link1--
%nproc=1
%chk=molecule
%mem=1024MB
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) iop(6/50=1) opt

remark line goes here

1   1
    C  -38.9208000000       15.3014000000      -24.0302000000   
    N  -39.4587000000       15.9888000000      -25.0480000000     
    C  -38.8204000000       17.2099000000      -25.1653000000     

antechamber-ini.esp

antechamber.esp
以及下面这种：
--Link1--
%nproc=1
%chk=molecule
%mem=1024MB
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt

remark line goes here

1   1
    C  -38.9208000000       15.3014000000      -24.0302000000   
    N  -39.4587000000       15.9888000000      -25.0480000000     
    C  -38.8204000000       17.2099000000      -25.1653000000   
不知道哪一种做法是对的？我感觉两种做出来的静电势很不一样。  

diaok

IOp(6/50)

Whether to write Antechamber file during ESP charge fitting.
这个关键词只是把esp的文件单独写出来
然后antechamber的时候用的输入的文件不一样
antechamber -L 或者 直接antechamber 有一堆各种信息

sobereva

靠iop(6/50=1) 产生esp文件是多余的，对于最新版本情况，antechamber直接就能从Gaussian输出文件中读取计算RESP电荷要用的静电势信息

2019-May-15补充：用ambertools计算RESP电荷的做法已经完全过时，目前算RESP电荷的最理想做法是用Multiwfn程序，极度快速、灵活和方便，见
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）