Journal
PHYSICAL CHEMISTRY CHEMICAL PHYSICS
Volume 10, Issue 34, Pages 5238-5249Publisher
ROYAL SOC CHEMISTRY
DOI: 10.1039/b802376j
Keywords
-
Funding
- EPSRC
Ask authors/readers for more resources
The question of how to compute acidity constants (pK(a)) treating solvent and solute at the same level of theory remains of some interest, for example in the case of high or low pH conditions. We have developed a density functional theory based molecular dynamics implementation of such a method. The method is based on a half reaction scheme computing free energies of dissociation from the vertical energy gaps for insertion or removal of protons. Finite system size effects are important, but largely cancel when half reactions are combined to full reactions. We verified the method by investigating a series of organic and inorganic acids and bases spanning a wide range of pK(a) values (20 units). We find that the response of the aqueous solvent to vertical protonation/deprotonation is almost always asymmetric and correlated with the strength of the hydrogen bonding of the deprotonated base. We interpret these observations in analogy with the picture of solvent response to electronic ionization.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available