# 自动计算ESP和RESP电荷(AMBER and G09)

## 高斯计算电荷和antechamber拟合电荷

### 关键词:

`#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2, 6/42=6, 6/50=1) opt`

• Pop是population analysis,用于输出分子轨道,电荷等信息.这里,MK(Merz-Kollman, also known as ESP) atomic charges will be printed out (atomic charges are fitted to the quantum-chemically generated electrostatic potential)
• `iop(6/33=2)`: coordinates of points along with the values of quantum-chemically calculated electrostatic potential in each of them will be printed out. 用于后续计算resp等.
• `iop(6/41=6)` : Number of layers in esp charge fit. 缺省4.也有设6的.现在的antechamber转换都不用这项, 旧教程有用6.
• `iop(6/42=6)` : Density of points per unit area in esp fit.缺省1. 一般使用6. 越大点越多.
• `iop(6/50=1)`: Whether to write Antechamber file during ESP charge fitting. 缺省0. 再算RESP的时候,正常需要使用该项, 这时在输入文件末尾写入两个文件名,可以分别储存开始结构和优化结构的ESP电荷文件.在G09B01中失效..
##### 输入文件例子:
``````--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

``````

``````# Get Gaussian Input file for esp/resp calculation
antechamber -i ligand.mol2 -fi mol2 -o ligand.gjf -fo gcrt -pf y -gn "%nproc=8" -gm "%mem=1000MB" -ch "ligand" -gk "#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) iop(6/50=1) opt" -ge ligand.gesp -gv 1

# Run Gaussian, can also use > ligand.out to assign output file
g09 ligand.gjf

# Use Gaussian output to get resp, also "-c esp" to get esp charge
antechamber -i ligand.log -fi gout -o ligand_resp.mol2 -fo mol2 -pf y -c resp
# The following use gesp file to get resp charge
# antechamber -i ligand.gesp -fi gesp -o ligand_resp.mol2 -fo mol2 -pf y -c resp

# To get amber parameter file
antechamber -i ligand_resp.mol2 -fi mol2 -o ligand_resp.prep -fo prepi
parmchk -i ligand_resp.prep -f prepi -o ligand_resp.frcmod
``````

• `-fo gcrt`指明输出高斯输入类型文件, `-fi`是输入文件类型, `-i`, `-o`分别是输入和输出文件, `-pf y`是不产生中间文件.
• `-gn`指明cpu数量, `-gm`是内存指明,`-ch`是chk文件名(默认”molecule”),均需要双引号!
• `-gk`是关键词指明(不输入就会自动产生我上面那句话); `-gv 1`是写出resp文件的选项,使用后会加入”#iop(6/50=1)”, g09开始支持, 也可以使用-gk自己写入; `-ge`是输出的gesp文件名(默认g09.gesp,写在高斯输入文件名后面),用于后续拟合,必须在`-gv 1`下才有效.
• `-c`指明拟合原子电荷的方法.

Gaussian 09 fix
In Gaussian09 rev B.01, the facility to write out the electrostatic potential on a grid of points was inadvertently deleted. This means that antechamber and resp jobs won’t work as they should. Fernando Clemente of Gaussian has kindly provided a script to work around the problem. Download the fixreadinesp.sh file, and follow the instructions there. (Note: you will have to make the script executable by typing `chmod +x fixreadinesp.sh`.)

``````#!/bin/bash
######################################################################
#
# "fixreadinesp.awk" by Fernando R. Clemente (Gaussian, Inc.)
#
#  Workaround for G09 rev. B.01 not printing the list of
#   electrostatic potential values at the grid points in
#   "Pop=ESP" jobs with "IOp(6/33=2)":
#
#  1) Run the "Pop=ESP" jobs with "IOp(6/33=2)" as you would
#     normally do.
#  2) Run this script on the output of job 1) and redirect to
#     a new Gaussian input file. This will extract the list of
#     grid points and setup a new job in which this list is given
#     explicitly in the input file. This second Gaussian job
#     uses the results saved in the checkpoint file of job 1),
#     so the QM calculation does not need to be redone. Example:
#
#
#  3) Run the new Gaussian job with the generated input file in
#     step 2) in order to obtain a G09 output file that contains
#     the values of the electrostatic potential at each grid point.
#
#  4) Re-run this script on the G09 output file generated by step 3)
#     to fix the format so "espgen" is able to read this type of
#     Gaussian output file. Example:
#
#
#  This final Gaussian output file should work with "espgen".
#
######################################################################
#
# Extracts list of grid points from a G09 output file
#  and sets up an Gaussian input file to re-run the job
#  with the list of points given explicitly.
#
if [ "\$(grep "Charges from ESP fit" \${1})" != "" ] ; then
printf "%s\n" "\$(grep -i "%chk=" \${1})"
printf "#p geom=allcheck chkbas guess=(read,only) density=check\n"
grep "ESP Fit Center" \${1} | cut -c32- -
printf "\n"
else
# Fixes the output of a G09 job with
#  read-in values to make "espgen" be
#  able to recognize the list of centers
#  and the electrostatic potential values.
awk '
{
if (\$0 ~ /Electric Field/) {
while (\$0 !~ /Atom/) { print \$0 ; getline }
while (\$0 ~ /Atom/) { print \$0 ; getline }
while (\$0 !~ /^ ------/) {
CENTERID=\$1
OLDSTR=sprintf("%i    ",CENTERID)
NEWSTR=sprintf("%i Fit",CENTERID)
sub(OLDSTR,NEWSTR,\$0)
print \$0
getline
}
}
print \$0
}' \${1}
fi
######################################################################
``````

``````chmod +x fixreadinesp.sh
g09 tmpesp.gjf > tmpesp.log
antechamber -i ligand.out -fi gout -o ligand_resp.mol2 -fo mol2 -pf y -c resp
``````

`````` %chk=ligand

1.208000  0.697100  2.100000
................
``````

Another Problem is the calculation result is very strange if you just use the standard method and the fixreadinesp.sh tool. Because of the use of “nosymm” in the second part extraction of fit charge from chk file, the coordinates may be different to the origin calculation. So, the final charge on atom summaried from the grid may be wrong. To solve this problem, just need to add nosymm in the first calculation, or remove the nosymm in the second gaussian input file generated by fixreadinesp.sh.

``````        OLDSTR=sprintf("%i    ",CENTERID)
NEWSTR=sprintf("%i Fit",CENTERID)
``````

``````        OLDSTR=sprintf("%s    ",CENTERID)
NEWSTR=sprintf("%s Fit",CENTERID)
``````

Update: 给作者写邮件指出问题后, 作者很快就回复了并给出了新的脚本, 效率很快啊~

``````#!/bin/bash
######################################################################
#
# "fixreadinesp.sh" by Fernando R. Clemente (Gaussian, Inc.)
#                      ( September 2015 )
#
#  Workaround for G09 rev. B.01 not printing the list of
#   electrostatic potential values at the grid points in
#   "Pop=ESP" jobs with "IOp(6/33=2)":
#
#  1) Run the "Pop=ESP" jobs with "IOp(6/33=2)" as you would
#     normally do.
#  2) Run this script on the output of job 1) and redirect to
#     a new Gaussian input file. This will extract the list of
#     grid points and setup a new job in which this list is given
#     explicitly in the input file. This second Gaussian job
#     uses the results saved in the checkpoint file of job 1),
#     so the QM calculation does not need to be redone. Example:
#
#
#  3) Run the new Gaussian job with the generated input file in
#     step 2) in order to obtain a G09 output file that contains
#     the values of the electrostatic potential at each grid point.
#
#  4) Re-run this script on the G09 output file generated by step 3)
#     to fix the format so "espgen" is able to read this type of
#     Gaussian output file. Example:
#
#
#  This final Gaussian output file should work with "espgen".
#
######################################################################
#
# Extracts list of grid points from a G09 output file
#  and sets up an Gaussian input file to re-run the job
#  with the list of points given explicitly.
#
if [ "\$(grep "Charges from ESP fit" \${1})" != "" ] ; then
printf "%s\n" "\$(grep -i "%chk=" \${1} | cut -c2-)"
printf "ESP fit with read-in centers\n\n"
awk '/Charge =... Multiplicity =../ { printf("%s %s\n",substr(\$0,10,3),substr(\$0,28,2)) }' \${1}
awk '
ATNO=0
while (\$1 !~ /^1/) { getline }
while (\$0 !~ /^ Generate VDW surfaces/) { ATNO++ ; AT[ATNO]=\$2 ; getline }
while (\$0 !~ /Atomic Center/) { getline }
ATNO=0
while (\$0 ~ /Atomic Center/) { ATNO++ ; XYZ[ATNO]=substr(\$0,32,31) ; getline }
NATOMS=ATNO
PT=0
while (\$0 ~ /ESP Fit Center/) { PT++ ; FITXYZ[PT]=substr(\$0,32,31) ; getline }
NPT=PT
}
END {
for (I=1;I<=NATOMS;I++) { printf("%3i  %s\n",AT[I],XYZ[I]) }
printf("\n")
for (I=1;I<=NPT;I++) { printf(" %s\n",FITXYZ[I]) }
printf("\n")
}' \${1}
else
# Fixes the output of a G09 job with
#  read-in values to make "espgen" be
#  able to recognize the list of centers
#  and the electrostatic potential values.
awk '
{
if (\$0 ~ /Z-Matrix orientation:/) { sub(/Z-Matrix orientation:/,"   Input orientation:",\$0) }
if (\$0 ~ /Electric Field/) {
while (\$0 !~ /Atom/) { print \$0 ; getline }
while (\$0 ~ /Atom/) { print \$0 ; getline }
while (\$0 !~ /^ ------/) {
FIRST=substr(\$0,1,5)
LAST=substr(\$0,11)
printf("%s Fit %s\n",FIRST,LAST)
getline
}
}
print \$0
}' \${1}
fi
######################################################################
``````
##### 自动进行ESP/RESP电荷计算脚本

Gaussian B.01使用. nosymm处理, 已更改fixreadinesp.sh中`%i``%s`.

``````#! /bin/bash
# Usage: Calculate Molecule ESP/RESP charge via Amber/G09.
# Input: molecule file in pdb/mol2 format
# Author: Zhixiong Zhao 2015-9-17

if [ -z \$1 ];then
echo "Please assign the input file!"
exit
fi

# If you use G09 Rev B.01, uncomment the following line to fix bug
# G09B01="TRUE"

#Setup the amber/Gaussian environment. You should modify by your own environment
if [ -z \$AMBERHOME ];then
source /mnt/home/zhaozx/AmberTools/amber.sh
fi

if [ -z \$g09root ];then
source /mnt/home/zhaozx/Software/g09/gau_setup.sh
fi

basename=\${1%.*}
exdname=\${1##*.}

if [ \$exdname = "mol2" ];then
antechamber -fi mol2 -fo gcrt -pf y -i \$1  -o \${basename}.gjf -gn "%nproc=8" -gm "%mem=1000MB" -gk "#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) iop(6/50=1) opt nosymm" -ch "\${basename}" -ge \${basename}.gesp -gv 1
elif [ \$exdname = "pdb" ];then
antechamber -fi pdb -fo gcrt -pf y -i \$1 -o \${basename}.gjf -gn "%nproc=8" -gm "%mem=1000MB" -gk "#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) iop(6/50=1) opt nosymm" -ch "\${basename}" -ge \${basename}.gesp -gv 1
else
echo "Only support for pdb or mol2 files!"
exit
fi

g09 \${basename}.gjf
antechamber -fi gout -fo mol2 -pf y -i \${basename}.log -o \${basename}_esp.mol2 -c esp

# Fix bug in G09 Rev B.01
if [ ! -z \$G09B01 ];then
# Generate bugfix tool fixreadinesp.sh developed by Fernando Clemente of Gaussian

if [ ! -f fixreadinesp.sh ];then
echo "if [ \"\\$(grep \"Charges from ESP fit\" \\${1})\" != \"\" ] ; then" >> fixreadinesp.sh
echo "  printf \"%s\n\" \"\\$(grep -i \"%chk=\" \\${1})\"" >> fixreadinesp.sh
echo "  grep \"ESP Fit Center\" \\${1} | cut -c32- -" >> fixreadinesp.sh
echo "  printf \"\n\"" >> fixreadinesp.sh
echo "  awk '" >> fixreadinesp.sh
echo "    if (\\$0 ~ /Electric Field/) { " >> fixreadinesp.sh
echo "      while (\\$0 !~ /Atom/) { print \\$0 ; getline } " >> fixreadinesp.sh
echo "      while (\\$0 ~ /Atom/) { print \\$0 ; getline }" >> fixreadinesp.sh
echo "      while (\\$0 !~ /^ ------/) {" >> fixreadinesp.sh
echo "        OLDSTR=sprintf(\"%s    \",CENTERID)" >> fixreadinesp.sh
echo "        NEWSTR=sprintf(\"%s Fit\",CENTERID)" >> fixreadinesp.sh
echo "        print \\$0" >> fixreadinesp.sh
echo "    print \\$0" >> fixreadinesp.sh
echo "  }' \\${1}" >> fixreadinesp.sh
fi

g09 tmpesp.gjf
antechamber -fi gout -fo mol2 -pf y -i \${basename}.out -o \${basename}_resp.mol2 -c resp
rm punch qout QOUT esout tmpesp.gjf tmpesp.log

else
# G09 Rev C01/D01
antechamber -fi gesp -fo mol2 -pf y -i \${basename}.gesp -o \${basename}_resp.mol2 -c resp
rm punch qout QOUT esout

fi # End different Version for resp
echo "Done for \$basename!"
``````

## 测试

1. D.01版相同分子测试,带nosymm
2. D.01版相同分子测试,不带nosymm
3. B.01版相同分子测试,带nosymm
4. B.01版相同分子测试,不带nosymm

### ESP, 跟log输出文件一致

``````      1 C1          1.1550    0.6590   -0.0420 ca        1 MOL     -0.192863
2 C2          2.4130    1.4800    0.0700 c3        1 MOL      0.371372
3 H1          2.1600    2.5350    0.0440 h1        1 MOL      0.001563
4 O1          3.2910    1.1610   -0.9900 oh        1 MOL     -0.673059
5 H2          4.0950    1.6510   -0.8880 ho        1 MOL      0.412009
6 H3          2.8990    1.2760    1.0210 h1        1 MOL     -0.026595
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.117862
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.115217
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.072899
10 S1         -2.8330    1.3820   -0.3880 sh        1 MOL     -0.300690
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.161915
12 C5         -1.1910   -0.8320   -0.1410 ca        1 MOL     -0.030450
13 H6         -2.0930   -1.4150   -0.1850 ha        1 MOL      0.108941
14 C6          0.0300   -1.4710   -0.0390 ca        1 MOL     -0.442495
15 H7          0.0680   -2.5460   -0.0010 ha        1 MOL      0.208991
16 C7          1.2150   -0.7420    0.0100 ca        1 MOL      0.498549
17 N1          2.4400   -1.3810    0.1590 nh        1 MOL     -0.958993
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.401655
19 H9          2.4130   -2.3620   -0.0210 hn        1 MOL      0.389897

1 C1          0.8450   -0.4030   -0.1980 ca        1 MOL     -0.204288
2 C2          1.9750   -1.3720   -0.4250 c3        1 MOL      0.357895
3 H1          1.5770   -2.3780   -0.5140 h1        1 MOL      0.006281
4 O1          2.8910   -1.3020    0.6480 oh        1 MOL     -0.670388
5 H2          3.6190   -1.8840    0.4770 ho        1 MOL      0.412521
6 H3          2.4830   -1.1300   -1.3560 h1        1 MOL     -0.021984
7 C3         -0.4610   -0.8520   -0.1310 ca        1 MOL     -0.099205
8 H4         -0.6550   -1.9080   -0.1960 ha        1 MOL      0.110730
9 C4         -1.5340    0.0210    0.0110 ca        1 MOL      0.062863
10 S1         -3.2040   -0.6010    0.1340 sh        1 MOL     -0.299507
11 H5         -3.5300   -0.5740   -1.1540 hs        1 MOL      0.162118
12 C5         -1.2700    1.3780    0.1110 ca        1 MOL     -0.026887
13 H6         -2.0810    2.0710    0.2370 ha        1 MOL      0.108836
14 C6          0.0290    1.8480    0.0610 ca        1 MOL     -0.445057
15 H7          0.2170    2.9040    0.1440 ha        1 MOL      0.208296
16 C7          1.1000    0.9730   -0.0910 ca        1 MOL      0.506061
17 N1          2.4020    1.4480   -0.1890 nh        1 MOL     -0.961354
18 H8          3.0790    0.8510    0.2360 hn        1 MOL      0.402173
19 H9          2.5130    2.3950    0.1020 hn        1 MOL      0.390895

1 C1          1.1550    0.6580   -0.0420 ca        1 MOL     -0.192720
2 C2          2.4130    1.4790    0.0710 c3        1 MOL      0.370896
3 H1          2.1600    2.5340    0.0440 h1        1 MOL      0.001774
4 O1          3.2910    1.1600   -0.9890 oh        1 MOL     -0.673062
5 H2          4.0950    1.6510   -0.8870 ho        1 MOL      0.412057
6 H3          2.8990    1.2760    1.0220 h1        1 MOL     -0.026466
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.117905
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.115115
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.073215
10 S1         -2.8330    1.3830   -0.3870 sh        1 MOL     -0.300670
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.161789
12 C5         -1.1920   -0.8320   -0.1410 ca        1 MOL     -0.030626
13 H6         -2.0930   -1.4150   -0.1860 ha        1 MOL      0.108998
14 C6          0.0300   -1.4710   -0.0400 ca        1 MOL     -0.442596
15 H7          0.0680   -2.5460   -0.0020 ha        1 MOL      0.209051
16 C7          1.2150   -0.7420    0.0090 ca        1 MOL      0.498516
17 N1          2.4390   -1.3820    0.1590 nh        1 MOL     -0.958883
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.401638
19 H9          2.4130   -2.3620   -0.0220 hn        1 MOL      0.389881

1 C1         -0.8450   -0.4030    0.1980 ca        1 MOL     -0.198304
2 C2         -1.9750   -1.3730    0.4250 c3        1 MOL      0.355612
3 H1         -1.5770   -2.3780    0.5140 h1        1 MOL      0.006917
4 O1         -2.8910   -1.3020   -0.6480 oh        1 MOL     -0.671186
5 H2         -3.6190   -1.8850   -0.4760 ho        1 MOL      0.413065
6 H3         -2.4830   -1.1310    1.3560 h1        1 MOL     -0.022323
7 C3          0.4610   -0.8520    0.1310 ca        1 MOL     -0.105486
8 H4          0.6550   -1.9080    0.1950 ha        1 MOL      0.113032
9 C4          1.5340    0.0210   -0.0110 ca        1 MOL      0.062586
10 S1          3.2040   -0.6010   -0.1340 sh        1 MOL     -0.299580
11 H5          3.5290   -0.5730    1.1550 hs        1 MOL      0.162243
12 C5          1.2700    1.3780   -0.1110 ca        1 MOL     -0.026940
13 H6          2.0810    2.0710   -0.2370 ha        1 MOL      0.108969
14 C6         -0.0290    1.8480   -0.0610 ca        1 MOL     -0.442838
15 H7         -0.2170    2.9040   -0.1440 ha        1 MOL      0.207678
16 C7         -1.1000    0.9730    0.0910 ca        1 MOL      0.505230
17 N1         -2.4020    1.4480    0.1890 nh        1 MOL     -0.962495
18 H8         -3.0790    0.8510   -0.2360 hn        1 MOL      0.402874
19 H9         -2.5130    2.3950   -0.1020 hn        1 MOL      0.390946
``````

### RESP,opt structure(from gesp)

``````      1 C1          1.1550    0.6590   -0.0420 ca        1 MOL     -0.156662
2 C2          2.4130    1.4800    0.0700 c3        1 MOL      0.325359
3 H1          2.1600    2.5350    0.0440 h1        1 MOL     -0.001825
4 O1          3.2910    1.1610   -0.9900 oh        1 MOL     -0.658991
5 H2          4.0950    1.6510   -0.8880 ho        1 MOL      0.410987
6 H3          2.8990    1.2760    1.0210 h1        1 MOL     -0.001825
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.125750
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.115495
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.073429
10 S1         -2.8330    1.3820   -0.3880 sh        1 MOL     -0.299423
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.161763
12 C5         -1.1910   -0.8320   -0.1410 ca        1 MOL     -0.039819
13 H6         -2.0930   -1.4150   -0.1850 ha        1 MOL      0.110328
14 C6          0.0300   -1.4710   -0.0390 ca        1 MOL     -0.415901
15 H7          0.0680   -2.5460   -0.0010 ha        1 MOL      0.200811
16 C7          1.2150   -0.7420    0.0100 ca        1 MOL      0.462337
17 N1          2.4400   -1.3810    0.1590 nh        1 MOL     -0.947570
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.393629
19 H9          2.4130   -2.3620   -0.0210 hn        1 MOL      0.393629

1 C1          0.8450   -0.4030   -0.1980 ca        1 MOL     -0.165459
2 C2          1.9750   -1.3720   -0.4250 c3        1 MOL      0.310672
3 H1          1.5770   -2.3780   -0.5140 h1        1 MOL      0.003116
4 O1          2.8910   -1.3020    0.6480 oh        1 MOL     -0.656692
5 H2          3.6190   -1.8840    0.4770 ho        1 MOL      0.411755
6 H3          2.4830   -1.1300   -1.3560 h1        1 MOL      0.003116
7 C3         -0.4610   -0.8520   -0.1310 ca        1 MOL     -0.109224
8 H4         -0.6550   -1.9080   -0.1960 ha        1 MOL      0.111286
9 C4         -1.5340    0.0210    0.0110 ca        1 MOL      0.065359
10 S1         -3.2040   -0.6010    0.1340 sh        1 MOL     -0.298425
11 H5         -3.5300   -0.5740   -1.1540 hs        1 MOL      0.161931
12 C5         -1.2700    1.3780    0.1110 ca        1 MOL     -0.038510
13 H6         -2.0810    2.0710    0.2370 ha        1 MOL      0.110611
14 C6          0.0290    1.8480    0.0610 ca        1 MOL     -0.416449
15 H7          0.2170    2.9040    0.1440 ha        1 MOL      0.199881
16 C7          1.1000    0.9730   -0.0910 ca        1 MOL      0.467910
17 N1          2.4020    1.4480   -0.1890 nh        1 MOL     -0.949833
18 H8          3.0790    0.8510    0.2360 hn        1 MOL      0.394478
19 H9          2.5130    2.3950    0.1020 hn        1 MOL      0.394478
``````

### RESP, ini structure(from gesp)

``````      1 C1          1.2080    0.6970    0.0000 ca        1 MOL     -0.257815
2 C2          2.4950    1.4410    0.0000 c3        1 MOL      0.340462
3 H1          2.2950    2.5120    0.0000 h1        1 MOL     -0.018413
4 O1          3.2380    1.0990   -1.1580 oh        1 MOL     -0.633372
5 H2          4.0610    1.5930   -1.1310 ho        1 MOL      0.395844
6 H3          3.0660    1.1780    0.8900 h1        1 MOL     -0.018413
7 C3          0.0000    1.3950    0.0000 ca        1 MOL     -0.107945
8 H4          0.0000    2.4820    0.0000 ha        1 MOL      0.218270
9 C4         -1.2080    0.6980    0.0000 ca        1 MOL     -0.031004
10 S1         -2.7370    1.5800    0.0000 sh        1 MOL     -0.403760
11 H5         -3.7520    0.7120    0.0000 hs        1 MOL      0.254530
12 C5         -1.2080   -0.6980    0.0000 ca        1 MOL     -0.031282
13 H6         -2.1490   -1.2410    0.0000 ha        1 MOL      0.136449
14 C6          0.0000   -1.3950    0.0000 ca        1 MOL     -0.453574
15 H7          0.0000   -2.4820    0.0000 ha        1 MOL      0.199507
16 C7          1.2080   -0.6980    0.0000 ca        1 MOL      0.631395
17 N1          2.4190   -1.3970    0.0000 nh        1 MOL     -1.110707
18 H8          3.2930   -0.8910    0.0000 hn        1 MOL      0.444914
19 H9          2.4190   -2.4070    0.0000 hn        1 MOL      0.444914

1 C1          0.8530   -0.4100   -0.2240 ca        1 MOL     -0.278993
2 C2          1.9900   -1.3480   -0.4080 c3        1 MOL      0.349096
3 H1          1.6110   -2.3670   -0.4880 h1        1 MOL     -0.019875
4 O1          2.8640   -1.2560    0.7050 oh        1 MOL     -0.632241
5 H2          3.5870   -1.8710    0.5590 ho        1 MOL      0.394561
6 H3          2.5320   -1.0890   -1.3170 h1        1 MOL     -0.019875
7 C3         -0.4530   -0.8990   -0.1790 ca        1 MOL     -0.095439
8 H4         -0.6380   -1.9650   -0.2790 ha        1 MOL      0.216303
9 C4         -1.5220   -0.0190   -0.0070 ca        1 MOL     -0.033467
10 S1         -3.1740   -0.6370    0.0490 sh        1 MOL     -0.403503
11 H5         -4.0240    0.3800    0.2210 hs        1 MOL      0.254100
12 C5         -1.2840    1.3500    0.1220 ca        1 MOL     -0.033933
13 H6         -2.1160    2.0360    0.2570 ha        1 MOL      0.139010
14 C6          0.0220    1.8380    0.0780 ca        1 MOL     -0.460436
15 H7          0.2070    2.9050    0.1780 ha        1 MOL      0.200621
16 C7          1.0900    0.9590   -0.0950 ca        1 MOL      0.647870
17 N1          2.3990    1.4480   -0.1390 nh        1 MOL     -1.113780
18 H8          3.1720    0.8100   -0.2650 hn        1 MOL      0.444990
19 H9          2.5710    2.4390   -0.0460 hn        1 MOL      0.444990
``````

### RESP, Fix.sh `%i`

``````      1 C1          1.1550    0.6590   -0.0420 ca        1 MOL     -0.161751
2 C2          2.4130    1.4800    0.0700 c3        1 MOL      0.333577
3 H1          2.1600    2.5350    0.0440 h1        1 MOL     -0.011242
4 O1          3.2910    1.1610   -0.9900 oh        1 MOL     -0.599179
5 H2          4.0950    1.6510   -0.8880 ho        1 MOL      0.371756
6 H3          2.8990    1.2760    1.0210 h1        1 MOL     -0.011242
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.148060
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.112207
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.128114
10 S1         -2.8330    1.3820   -0.3880 sh        1 MOL     -0.276854
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.142380
12 C5         -1.1910   -0.8320   -0.1410 ca        1 MOL     -0.105482
13 H6         -2.0930   -1.4150   -0.1850 ha        1 MOL      0.116274
14 C6          0.0300   -1.4710   -0.0390 ca        1 MOL     -0.315181
15 H7          0.0680   -2.5460   -0.0010 ha        1 MOL      0.159579
16 C7          1.2150   -0.7420    0.0100 ca        1 MOL      0.419384
17 N1          2.4400   -1.3810    0.1590 nh        1 MOL     -0.838418
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.342069
19 H9          2.4130   -2.3620   -0.0210 hn        1 MOL      0.342069

1 C1          0.8450   -0.4030   -0.1980 ca        1 MOL     -0.197801
2 C2          1.9750   -1.3720   -0.4250 c3        1 MOL      0.332137
3 H1          1.5770   -2.3780   -0.5140 h1        1 MOL     -0.009912
4 O1          2.8910   -1.3020    0.6480 oh        1 MOL     -0.598844
5 H2          3.6190   -1.8840    0.4770 ho        1 MOL      0.372008
6 H3          2.4830   -1.1300   -1.3560 h1        1 MOL     -0.009912
7 C3         -0.4610   -0.8520   -0.1310 ca        1 MOL     -0.093596
8 H4         -0.6550   -1.9080   -0.1960 ha        1 MOL      0.104183
9 C4         -1.5340    0.0210    0.0110 ca        1 MOL      0.057735
10 S1         -3.2040   -0.6010    0.1340 sh        1 MOL     -0.270447
11 H5         -3.5300   -0.5740   -1.1540 hs        1 MOL      0.150700
12 C5         -1.2700    1.3780    0.1110 ca        1 MOL     -0.061126
13 H6         -2.0810    2.0710    0.2370 ha        1 MOL      0.110181
14 C6          0.0290    1.8480    0.0610 ca        1 MOL     -0.339422
15 H7          0.2170    2.9040    0.1440 ha        1 MOL      0.159005
16 C7          1.1000    0.9730   -0.0910 ca        1 MOL      0.450824
17 N1          2.4020    1.4480   -0.1890 nh        1 MOL     -0.839315
18 H8          3.0790    0.8510    0.2360 hn        1 MOL      0.341801
19 H9          2.5130    2.3950    0.1020 hn        1 MOL      0.341801

1 C1          1.1550    0.6580   -0.0420 ca        1 MOL     -0.162088
2 C2          2.4130    1.4790    0.0710 c3        1 MOL      0.333112
3 H1          2.1600    2.5340    0.0440 h1        1 MOL     -0.011087
4 O1          3.2910    1.1600   -0.9890 oh        1 MOL     -0.599158
5 H2          4.0950    1.6510   -0.8870 ho        1 MOL      0.371796
6 H3          2.8990    1.2760    1.0220 h1        1 MOL     -0.011087
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.147541
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.111988
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.127889
10 S1         -2.8330    1.3830   -0.3870 sh        1 MOL     -0.276796
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.142292
12 C5         -1.1920   -0.8320   -0.1410 ca        1 MOL     -0.105290
13 H6         -2.0930   -1.4150   -0.1860 ha        1 MOL      0.116207
14 C6          0.0300   -1.4710   -0.0400 ca        1 MOL     -0.315288
15 H7          0.0680   -2.5460   -0.0020 ha        1 MOL      0.159498
16 C7          1.2150   -0.7420    0.0090 ca        1 MOL      0.419934
17 N1          2.4390   -1.3820    0.1590 nh        1 MOL     -0.838427
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.342023
19 H9          2.4130   -2.3620   -0.0220 hn        1 MOL      0.342023

1 C1         -0.8450   -0.4030    0.1980 ca        1 MOL     -0.002460
2 C2         -1.9750   -1.3730    0.4250 c3        1 MOL     -0.357360
3 H1         -1.5770   -2.3780    0.5140 h1        1 MOL      0.209924
4 O1         -2.8910   -1.3020   -0.6480 oh        1 MOL     -0.518857
5 H2         -3.6190   -1.8850   -0.4760 ho        1 MOL      0.379120
6 H3         -2.4830   -1.1310    1.3560 h1        1 MOL      0.209924
7 C3          0.4610   -0.8520    0.1310 ca        1 MOL      0.079757
8 H4          0.6550   -1.9080    0.1950 ha        1 MOL      0.006852
9 C4          1.5340    0.0210   -0.0110 ca        1 MOL     -0.259404
10 S1          3.2040   -0.6010   -0.1340 sh        1 MOL     -0.160999
11 H5          3.5290   -0.5730    1.1550 hs        1 MOL      0.146076
12 C5          1.2700    1.3780   -0.1110 ca        1 MOL      0.191306
13 H6          2.0810    2.0710   -0.2370 ha        1 MOL      0.034480
14 C6         -0.0290    1.8480   -0.0610 ca        1 MOL     -0.457445
15 H7         -0.2170    2.9040   -0.1440 ha        1 MOL      0.173420
16 C7         -1.1000    0.9730    0.0910 ca        1 MOL      0.469661
17 N1         -2.4020    1.4480    0.1890 nh        1 MOL     -0.879408
18 H8         -3.0790    0.8510   -0.2360 hn        1 MOL      0.367707
19 H9         -2.5130    2.3950   -0.1020 hn        1 MOL      0.367707
``````

### RESP, Fix.sh `%s`

``````      1 C1          1.1550    0.6590   -0.0420 ca        1 MOL     -0.130447
2 C2          2.4130    1.4800    0.0700 c3        1 MOL      0.333350
3 H1          2.1600    2.5350    0.0440 h1        1 MOL     -0.003752
4 O1          3.2910    1.1610   -0.9900 oh        1 MOL     -0.662975
5 H2          4.0950    1.6510   -0.8880 ho        1 MOL      0.411752
6 H3          2.8990    1.2760    1.0210 h1        1 MOL     -0.003752
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.180915
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.120993
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.157206
10 S1         -2.8330    1.3820   -0.3880 sh        1 MOL     -0.306690
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.154485
12 C5         -1.1910   -0.8320   -0.1410 ca        1 MOL     -0.101430
13 H6         -2.0930   -1.4150   -0.1850 ha        1 MOL      0.119336
14 C6          0.0300   -1.4710   -0.0390 ca        1 MOL     -0.379298
15 H7          0.0680   -2.5460   -0.0010 ha        1 MOL      0.196514
16 C7          1.2150   -0.7420    0.0100 ca        1 MOL      0.433842
17 N1          2.4400   -1.3810    0.1590 nh        1 MOL     -0.946080
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.393930
19 H9          2.4130   -2.3620   -0.0210 hn        1 MOL      0.393930

1 C1          0.8450   -0.4030   -0.1980 ca        1 MOL     -0.178577
2 C2          1.9750   -1.3720   -0.4250 c3        1 MOL      0.334818
3 H1          1.5770   -2.3780   -0.5140 h1        1 MOL     -0.003114
4 O1          2.8910   -1.3020    0.6480 oh        1 MOL     -0.661250
5 H2          3.6190   -1.8840    0.4770 ho        1 MOL      0.411407
6 H3          2.4830   -1.1300   -1.3560 h1        1 MOL     -0.003114
7 C3         -0.4610   -0.8520   -0.1310 ca        1 MOL     -0.111843
8 H4         -0.6550   -1.9080   -0.1960 ha        1 MOL      0.110821
9 C4         -1.5340    0.0210    0.0110 ca        1 MOL      0.074323
10 S1         -3.2040   -0.6010    0.1340 sh        1 MOL     -0.300720
11 H5         -3.5300   -0.5740   -1.1540 hs        1 MOL      0.163606
12 C5         -1.2700    1.3780    0.1110 ca        1 MOL     -0.043801
13 H6         -2.0810    2.0710    0.2370 ha        1 MOL      0.112099
14 C6          0.0290    1.8480    0.0610 ca        1 MOL     -0.423226
15 H7          0.2170    2.9040    0.1440 ha        1 MOL      0.201747
16 C7          1.1000    0.9730   -0.0910 ca        1 MOL      0.482003
17 N1          2.4020    1.4480   -0.1890 nh        1 MOL     -0.955682
18 H8          3.0790    0.8510    0.2360 hn        1 MOL      0.395252
19 H9          2.5130    2.3950    0.1020 hn        1 MOL      0.395252

1 C1          1.1550    0.6580   -0.0420 ca        1 MOL     -0.130494
2 C2          2.4130    1.4790    0.0710 c3        1 MOL      0.332882
3 H1          2.1600    2.5340    0.0440 h1        1 MOL     -0.003588
4 O1          3.2910    1.1600   -0.9890 oh        1 MOL     -0.662949
5 H2          4.0950    1.6510   -0.8870 ho        1 MOL      0.411785
6 H3          2.8990    1.2760    1.0220 h1        1 MOL     -0.003588
7 C3         -0.0770    1.2760   -0.1580 ca        1 MOL     -0.180700
8 H4         -0.1210    2.3480   -0.2150 ha        1 MOL      0.120902
9 C4         -1.2630    0.5510   -0.1970 ca        1 MOL      0.157062
10 S1         -2.8330    1.3830   -0.3870 sh        1 MOL     -0.306643
11 H5         -3.1360    1.5470    0.8960 hs        1 MOL      0.154428
12 C5         -1.1920   -0.8320   -0.1410 ca        1 MOL     -0.101181
13 H6         -2.0930   -1.4150   -0.1860 ha        1 MOL      0.119296
14 C6          0.0300   -1.4710   -0.0400 ca        1 MOL     -0.379635
15 H7          0.0680   -2.5460   -0.0020 ha        1 MOL      0.196603
16 C7          1.2150   -0.7420    0.0090 ca        1 MOL      0.434009
17 N1          2.4390   -1.3820    0.1590 nh        1 MOL     -0.946042
18 H8          3.1850   -0.9370   -0.3330 hn        1 MOL      0.393926
19 H9          2.4130   -2.3620   -0.0220 hn        1 MOL      0.393926

1 C1         -0.8450   -0.4030    0.1980 ca        1 MOL     -0.000094
2 C2         -1.9750   -1.3730    0.4250 c3        1 MOL     -0.344126
3 H1         -1.5770   -2.3780    0.5140 h1        1 MOL      0.211843
4 O1         -2.8910   -1.3020   -0.6480 oh        1 MOL     -0.567602
5 H2         -3.6190   -1.8850   -0.4760 ho        1 MOL      0.411686
6 H3         -2.4830   -1.1310    1.3560 h1        1 MOL      0.211843
7 C3          0.4610   -0.8520    0.1310 ca        1 MOL      0.086938
8 H4          0.6550   -1.9080    0.1950 ha        1 MOL      0.010714
9 C4          1.5340    0.0210   -0.0110 ca        1 MOL     -0.282636
10 S1          3.2040   -0.6010   -0.1340 sh        1 MOL     -0.179106
11 H5          3.5290   -0.5730    1.1550 hs        1 MOL      0.157669
12 C5          1.2700    1.3780   -0.1110 ca        1 MOL      0.243484
13 H6          2.0810    2.0710   -0.2370 ha        1 MOL      0.025461
14 C6         -0.0290    1.8480   -0.0610 ca        1 MOL     -0.544482
15 H7         -0.2170    2.9040   -0.1440 ha        1 MOL      0.215609
16 C7         -1.1000    0.9730    0.0910 ca        1 MOL      0.490765
17 N1         -2.4020    1.4480    0.1890 nh        1 MOL     -0.967096
18 H8         -3.0790    0.8510   -0.2360 hn        1 MOL      0.409565
19 H9         -2.5130    2.3950   -0.1020 hn        1 MOL      0.409565
``````

1. D0.1和B0.2的坐标结果基本一致, 电势结果也基本一致,偏差较小.
2. 使用nosymm和不适用nosymm结果有差异, 建议使用同一套方法处理即可. 要是使用B01,建议加入nosymm(加入后的值更准).
3. 由原子14起, fix工具使用`%i``%s`处理就开始出现较大偏差(约0.5,可能由于末尾`****`没加入fit所以不能正常被解析,导致最后一些坐标电荷缺失,从而影响结果.因此fix工具中使用`%s`将更好.这个问题可能作者都没有注意到..用fix工具出来的resp实际是优化后结构的,和正常途径获得的结果略有差异,但不大,主要可能是正常路径用于计算的charge center大概是6000多,而fix工具出来的约是12000左右,点更密集的缘故吧. Anyway, 还是推荐使用正常方法出来的结果.

## Reference

◆ 本文地址: http://platinhom.github.io/2015/09/17/AutoCalcRESP/, 转载请注明 ◆

/ 已阅读()
Source 类别: CompCB