CCL: Solvation Free Energy of Proton in THF

When I modeled protonated THF, I got more reasonable numbers. However, I observe a significant deviationÂÂfor the anionic system that makes me think about the anionic nature of the complex.

Thanks a lot for all suggestions.


On Wed, Aug 27, 2014 at 2:33 PM, Andreas Klamt klamt/./ <owner-chemistry(!)> wrote:

Sent to CCL by: Andreas Klamt [klamt^^]
Dear Busra,

obviously Mariusz is absolutely right. If you really should like to try to calculate absolute pKa, then you should get aware that the state of the proton in THF surely is not an isolated proton, but a protonated THF. As in water you have H3O+.

But I also would not expect that you can get an accurate prediction of the total free energy of the protonated H3O+ with a continuum model.

Best regards


Am 27.08.2014 11:38, schrieb Mariusz Radon mariusz.radon[A]

Sent to CCL by: Mariusz Radon []
On 08/25/2014 08:23 PM, BUSRA DERELI wrote:
Sent to CCL by: "BUSRAÂ DERELI" []

I am trying to calculate pKa values of two systems, one is anionic and the
other is cationic copper systems. I got so weird numbers not close to the
experimental values when I computed the free energy of proton in THF which is
the solvent of interest. I use M11-L local functional with SDD basis set on Cu
atom and 6-31G(d) on all other atoms in the system. Thus, I want to compare the
computed solvation free energy of proton in THF with the experimental value. In
some literature reports, I have seen absolute solvation free energy of proton
in a couple of solvents but not in THF.

Busra Dereli

Dear Busra:

Not very strange that you get weird number for solvation energy of
proton which (unlike other ions) is just a point charge! As far as I
know, it is extremely difficult to reproduce solvation energy of proton
in such calculations.

If your goal is to compare acidity of two substances, you are probably
more interested in the difference of their pKa values, than in the
individual pKa values. The point is that when you compute the
difference, the (huge) error on the proton solvation energy cancels out
exactly. Alternatively, you can compute relative pKa values of both
substances with respect to some other chemical, treated as reference,
whose actual pKa you take from experiment. That is much easier than
trying to reproduce absolute pKa values.

Best regards,
Mariusz Radon

Prof. Dr. Andreas Klamt
CEO / GeschÃftsfÃhrer
COSMOlogic GmbH & Co. KG
Imbacher Weg 46
D-51379 Leverkusen, Germany

phone Â+49-2171-731681
fax  Â+49-2171-731689
web  Â

[University address:Â Â Â Inst. of Physical and
Theoretical Chemistry, University of Regensburg]

HRA 20653 Amtsgericht Koeln, GF: Prof. Dr. Andreas Klamt
Komplementaer: COSMOlogic Verwaltungs GmbH
HRB 49501 Amtsgericht Koeln, GF: Prof. Dr. Andreas Klamt

-= This is automatically added to each message by the mailing script =-
E-mail to subscribers: CHEMISTRY(!) or use:

E-mail to administrators: CHEMISTRY-REQUEST(!) or use

Subscribe/Unsubscribe:Â Â Â

Before posting, check wait time at:

Job: Conferences:

Search Messages: