[AMBER] Hydrogen atoms with negative RESP charges from Antechamber

From: 库米尔 <418025.mail.dhu.edu.cn>
Date: Thu, 2 Sep 2021 23:23:24 +0800 (GMT+08:00)

Hello Amber Community
I have optimized ethanol molecule at HF/6-31g(d) level and calculated the esp charges at the same level with following guassian keywords Scf=tight pop=mk Iop(6/32=2, 6/41=10,6/42=17). Antechamber gives a resp charges below as shown below
.<TRIPOS>MOLECULE
ETH
    9 8 1 0 0
SMALL
resp




.<TRIPOS>ATOM
      1 C1 -0.0290 0.5450 0.0000 c3 1 ETH 0.566448
      2 C2 1.1650 -0.3870 0.0000 c3 1 ETH -0.216837
      3 O1 -1.1940 -0.2390 0.0000 oh 1 ETH -0.723897
      4 H1 0.0010 1.1890 0.8780 h1 1 ETH -0.099119
      5 H2 0.0010 1.1890 -0.8780 h1 1 ETH -0.099119
      6 H3 1.1470 -1.0220 0.8790 hc 1 ETH 0.052386
      7 H4 1.1470 -1.0220 -0.8790 hc 1 ETH 0.052386
      8 H5 2.0920 0.1780 0.0000 hc 1 ETH 0.052386
      9 H6 -1.9570 0.3220 0.0000 ho 1 ETH 0.415368
.<TRIPOS>BOND
     1 1 2 1
     2 1 3 1
     3 1 4 1
     4 1 5 1
     5 2 6 1
     6 2 7 1
     7 2 8 1
     8 3 9 1
.<TRIPOS>SUBSTRUCTURE
     1 ETH 1 TEMP 0 **** **** 0 ROOT

My concern is the negative charges on the hydrogen attached to the C1 makes the results appear quite unreasonable to me. Again, when i use the same method for 1-butyl-3-methyl-3-imidazolium chloride ionic liquid according united atom method employed by Liu et.al J.Phys.Chem.B.2010,114,4572,4572-4582, it produces similar negative charges on the hydrogen atom attached to the C2 carbon on the alkyl chain. I must add that no published paper have reported such negative charges in 1-butyl-3-methyl-3-imidazolium chloride ionic liquid. Could somebody help me understand why am getting these charges. Thanks in advance.
Alex




_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Sep 02 2021 - 08:30:02 PDT
Custom Search