计算化学公社

标题: 求助resp电荷的模拟 [打印本页]

作者
Author: kjxwl3 时间: 2017-9-21 14:18
标题: 求助resp电荷的模拟
各位大神，
   我在进行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
不知道哪一种做法是对的？我感觉两种做出来的静电势很不一样。

作者
Author: diaok 时间: 2017-9-21 15:24
IOp(6/50)

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

作者
Author: sobereva 时间: 2017-9-21 18:52
靠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）

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