Benchmarking Continuum Solvent Models for Keto-Enol Tautomerizations

Billy W. McCann, Stuart McFarland, Orlando Acevedo

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

Experimental free energies of tautomerization, ΔGT, were used to benchmark the gas-phase predictions of 17 different quantum mechanical methods and eight basis sets for seven keto-enol tautomer pairs dominated by their enolic form. The G4 method and M06/6-31+G(d,p) yielded the most accurate results, with mean absolute errors (MAE's) of 0.95 and 0.71 kcal/mol, respectively. Using these two theory levels, the solution-phase ΔGT values for 23 unique tautomer pairs composed of aliphatic ketones, β-dicarbonyls, and heterocycles were computed in multiple protic and aprotic solvents. The continuum solvation models, namely, polarizable continuum model (PCM), polarizable conductor calculation model (CPCM), and universal solvation model (SMD), gave relatively similar MAE's of ∼1.6-1.7 kcal/mol for G4 and ∼1.9-2.0 kcal/mol with M06/6-31+G(d,p). Partitioning the tautomer pairs into their respective molecular types, that is, aliphatic ketones, β-dicarbonyls, and heterocycles, and separating out the aqueous versus nonaqueous results finds G4/PCM utilizing the UA0 cavity to be the overall most accurate combination. Free energies of activation, ΔG, for the base-catalyzed keto-enol interconversion of 2-nitrocyclohexanone were also computed using six bases and five solvents. The M06/6-31+G(d,p) reproduced the ΔG with MAE's of 1.5 and 1.8 kcal/mol using CPCM and SMD, respectively, for all combinations of base and solvent. That specific enolization was previously proposed to proceed via a concerted mechanism in less polar solvents but shift to a stepwise mechanism in more polar solvents. However, the current calculations suggest that the stepwise mechanism operates in all solvents.

Original languageEnglish (US)
Pages (from-to)8724-8733
Number of pages10
JournalJournal of Physical Chemistry A
Volume119
Issue number32
DOIs
StatePublished - Jul 23 2015
Externally publishedYes

Fingerprint

Benchmarking
continuums
tautomers
Solvation
solvation
Ketones
ketones
Free energy
conductors
free energy
Pulse code modulation
Gases
Chemical activation
activation
vapor phases
cavities
shift
predictions

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Benchmarking Continuum Solvent Models for Keto-Enol Tautomerizations. / McCann, Billy W.; McFarland, Stuart; Acevedo, Orlando.

In: Journal of Physical Chemistry A, Vol. 119, No. 32, 23.07.2015, p. 8724-8733.

Research output: Contribution to journalArticle

McCann, Billy W. ; McFarland, Stuart ; Acevedo, Orlando. / Benchmarking Continuum Solvent Models for Keto-Enol Tautomerizations. In: Journal of Physical Chemistry A. 2015 ; Vol. 119, No. 32. pp. 8724-8733.
@article{f5eb177bcd264db49d7ae42e62b8e4f2,
title = "Benchmarking Continuum Solvent Models for Keto-Enol Tautomerizations",
abstract = "Experimental free energies of tautomerization, ΔGT, were used to benchmark the gas-phase predictions of 17 different quantum mechanical methods and eight basis sets for seven keto-enol tautomer pairs dominated by their enolic form. The G4 method and M06/6-31+G(d,p) yielded the most accurate results, with mean absolute errors (MAE's) of 0.95 and 0.71 kcal/mol, respectively. Using these two theory levels, the solution-phase ΔGT values for 23 unique tautomer pairs composed of aliphatic ketones, β-dicarbonyls, and heterocycles were computed in multiple protic and aprotic solvents. The continuum solvation models, namely, polarizable continuum model (PCM), polarizable conductor calculation model (CPCM), and universal solvation model (SMD), gave relatively similar MAE's of ∼1.6-1.7 kcal/mol for G4 and ∼1.9-2.0 kcal/mol with M06/6-31+G(d,p). Partitioning the tautomer pairs into their respective molecular types, that is, aliphatic ketones, β-dicarbonyls, and heterocycles, and separating out the aqueous versus nonaqueous results finds G4/PCM utilizing the UA0 cavity to be the overall most accurate combination. Free energies of activation, ΔG‡, for the base-catalyzed keto-enol interconversion of 2-nitrocyclohexanone were also computed using six bases and five solvents. The M06/6-31+G(d,p) reproduced the ΔG‡ with MAE's of 1.5 and 1.8 kcal/mol using CPCM and SMD, respectively, for all combinations of base and solvent. That specific enolization was previously proposed to proceed via a concerted mechanism in less polar solvents but shift to a stepwise mechanism in more polar solvents. However, the current calculations suggest that the stepwise mechanism operates in all solvents.",
author = "McCann, {Billy W.} and Stuart McFarland and Orlando Acevedo",
year = "2015",
month = "7",
day = "23",
doi = "10.1021/acs.jpca.5b04116",
language = "English (US)",
volume = "119",
pages = "8724--8733",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "32",

}

TY - JOUR

T1 - Benchmarking Continuum Solvent Models for Keto-Enol Tautomerizations

AU - McCann, Billy W.

AU - McFarland, Stuart

AU - Acevedo, Orlando

PY - 2015/7/23

Y1 - 2015/7/23

N2 - Experimental free energies of tautomerization, ΔGT, were used to benchmark the gas-phase predictions of 17 different quantum mechanical methods and eight basis sets for seven keto-enol tautomer pairs dominated by their enolic form. The G4 method and M06/6-31+G(d,p) yielded the most accurate results, with mean absolute errors (MAE's) of 0.95 and 0.71 kcal/mol, respectively. Using these two theory levels, the solution-phase ΔGT values for 23 unique tautomer pairs composed of aliphatic ketones, β-dicarbonyls, and heterocycles were computed in multiple protic and aprotic solvents. The continuum solvation models, namely, polarizable continuum model (PCM), polarizable conductor calculation model (CPCM), and universal solvation model (SMD), gave relatively similar MAE's of ∼1.6-1.7 kcal/mol for G4 and ∼1.9-2.0 kcal/mol with M06/6-31+G(d,p). Partitioning the tautomer pairs into their respective molecular types, that is, aliphatic ketones, β-dicarbonyls, and heterocycles, and separating out the aqueous versus nonaqueous results finds G4/PCM utilizing the UA0 cavity to be the overall most accurate combination. Free energies of activation, ΔG‡, for the base-catalyzed keto-enol interconversion of 2-nitrocyclohexanone were also computed using six bases and five solvents. The M06/6-31+G(d,p) reproduced the ΔG‡ with MAE's of 1.5 and 1.8 kcal/mol using CPCM and SMD, respectively, for all combinations of base and solvent. That specific enolization was previously proposed to proceed via a concerted mechanism in less polar solvents but shift to a stepwise mechanism in more polar solvents. However, the current calculations suggest that the stepwise mechanism operates in all solvents.

AB - Experimental free energies of tautomerization, ΔGT, were used to benchmark the gas-phase predictions of 17 different quantum mechanical methods and eight basis sets for seven keto-enol tautomer pairs dominated by their enolic form. The G4 method and M06/6-31+G(d,p) yielded the most accurate results, with mean absolute errors (MAE's) of 0.95 and 0.71 kcal/mol, respectively. Using these two theory levels, the solution-phase ΔGT values for 23 unique tautomer pairs composed of aliphatic ketones, β-dicarbonyls, and heterocycles were computed in multiple protic and aprotic solvents. The continuum solvation models, namely, polarizable continuum model (PCM), polarizable conductor calculation model (CPCM), and universal solvation model (SMD), gave relatively similar MAE's of ∼1.6-1.7 kcal/mol for G4 and ∼1.9-2.0 kcal/mol with M06/6-31+G(d,p). Partitioning the tautomer pairs into their respective molecular types, that is, aliphatic ketones, β-dicarbonyls, and heterocycles, and separating out the aqueous versus nonaqueous results finds G4/PCM utilizing the UA0 cavity to be the overall most accurate combination. Free energies of activation, ΔG‡, for the base-catalyzed keto-enol interconversion of 2-nitrocyclohexanone were also computed using six bases and five solvents. The M06/6-31+G(d,p) reproduced the ΔG‡ with MAE's of 1.5 and 1.8 kcal/mol using CPCM and SMD, respectively, for all combinations of base and solvent. That specific enolization was previously proposed to proceed via a concerted mechanism in less polar solvents but shift to a stepwise mechanism in more polar solvents. However, the current calculations suggest that the stepwise mechanism operates in all solvents.

UR - http://www.scopus.com/inward/record.url?scp=84939207944&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84939207944&partnerID=8YFLogxK

U2 - 10.1021/acs.jpca.5b04116

DO - 10.1021/acs.jpca.5b04116

M3 - Article

AN - SCOPUS:84939207944

VL - 119

SP - 8724

EP - 8733

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 32

ER -