The Journal of Organic Chemistry
Article
orienting the methyl groups to the boundaries of the cavity where the
solvated complex is placed within the framework of the PCM model
used. The water molecule in this complex is coordinated by one of the
hydrogen atoms to the oxygen of the DMSO molecule, and the
second hydrogen atom is coordinated to the atom, which is bound to
EXPERIMENTAL SECTION
■
Computational Details. Estimation of the enthalpy of the aldol
reaction is associated with certain difficulties when using many
popular density functionals.31 In particular, for the aldol reaction of
two acetone molecules, the enthalpy change calculated on the basis of
the heats of formation is ΔH = −10.5 0.2 kcal/mol; the CBS-QB3
method predicts an exothermicity of the reaction ΔH = −9.2 kcal/
mol. At the same time, the B3LYP/6-31+G (d, p) approach estimates
this reaction as endothermic with ΔH = +1.0 kcal/mol. Extension of
the basis set to the triple-zeta further deteriorates the estimatefor
B3LYP/6-311+G (2df, 2p), ΔH increases to +2.4 kcal/mol. The same
trend is typical for such functionals as B3PW91, B1B95,
MPW1PW91, PBE1PBE, and only M06-2X gives a sufficiently
adequate value. The latter is most likely due to the fact that this
reaction was included in the training set when parametrizing this
functional.31
the potassium cation in Nu·K+ (see e.g.,, the transition states TS2′→3
,
TS6″→7, and TS7→8 in Figures 1 and 2). Thus, when searching for
minima and transition states on the PES, we considered the
conformations that concerned only the Nu·K+ fragment. The number
of localized conformations at the stationary points on the PES usually
did not exceed three. The search for transition states, characterized by
a variety of possible mutual orientations of the reagents (trans-
formations 2′ → 3, 6″ → 7, 8′ → 9b, and 6c′ → 9c), was carried out
in a special way. In the cases 2′ → 3, 8′ → 9b, and 6c′ → 9c, scanning
of the dihedral angle involving the CO bond of the acetone
molecule 1 and the CC double bond of the anions 2, 8, and 6c was
carried out, while for the search along the path 6″ → 7, the dihedral
angle involving the CC double bonds of croton 6 and enolate ion 2
was scanned. In all the cases, the interatomic distances between the
bonded atoms were additionally frozen at 2.05 Å, which corresponds
to a somewhat larger distance that is characteristic of the transition
state geometry. The minima obtained for a particular PES cross
section were selected for the optimization of the transition states
followed by energy refinement using the adopted B2PLYP//B3LYP
approach.
The B2PLYP functional,40 which provides “chemical” accuracy in
various tests,41−43 was not among the functionals studied in ref 31. In
our recent studies, the B2PLYP/6-311+G**//B3LYP/6-31+G*
approach has been successfully used to describe various reactions of
ketones with acetylenes carried out in superbasic media, including the
steps of vinylation and ethynylation, isomerization of the multiple
bond position, reactions of nucleophilic addition to a double bond, in
which unsaturated ketones act as the Michael acceptors, and so
forth.38 For the aldol reaction of acetone in the gas phase, this
approach with the inclusion of the D2 Grimme dispersion
correction44 gives the estimate ΔH = −8.2 kcal/mol, which is in
good agreement with the data of the reference calculation CBS-QB3
(ΔH = −9.2 kcal/mol). For this reason, the B2PLYP(D2)/6-
311+G**//B3LYP/6-31+G* approach was chosen to describe the
energy profiles of the aldol−croton condensation.
Optimization of the structural parameters of the molecular species
involved in the reaction was carried out, accounting for solvent effects
at the level of the polarization continuum model IEFPCM45 using
density functional theory at the B3LYP46,47 level of theory with the 6-
31+G* basis set. The vibrational corrections to enthalpies and Gibbs
free energies were calculated at the same level of theory (B3LYP/6-
31+G*) for a standard temperature 298.15 K.
For all stationary points, the number of negative eigenvalues of the
Hessian matrix was analyzed; the connection of the transition states
with the corresponding potential energy surface (PES) minima was
proved by the reaction coordinate following, using the local quadratic
approximation algorithm.48
The energies of the stationary points were further refined using the
double-hybrid functional B2PLYP with Grimme’s dispersion
correction D249 in combination with the 6-311+G** basis set. In
addition, the solvation energies in DMSO were calculated within the
framework of the polarizable continuum model IEFPCM,45 using the
B3LYP/6-31+G* approach. Modeling of the cavity was performed
within the approach GEPOL50 using the universal force field radii and
universal scaling factor for the solvent α = 1.1.51 For the calculation of
the solvation energy, the dimethyl sulfoxide dielectric constant ε =
46.8 was used.
All the kinetic curves were modeled at T = 353 K to the Arrhenius
equation k(T) = kBT/h × exp(−ΔG⧧/RT), using the KINET
program.39 All quantum chemistry calculations were performed with
Gaussian 09.54
Reaction of Acetone with Acetylene in the KOH/DMSO
System. A mixture of acetone (5.81 g, 0.1 mol) and KOH·0.5H2O
(6.50 g, 0.1 mol) in DMSO (50 mL) was placed into a 0.25 dm3 steel
rotating autoclave. The autoclave was fed with acetylene up to 14 atm
and then decompressed to atmospheric pressure to remove air. The
autoclave was fed with acetylene again (initial pressure at ambient
temperature was 14 atm) and heated to 80 °C for 1 h. The pressure
reached its maximum (18−20 atm) as the temperature was raised to
80 °C and then dropped upon acetylene consumption in the course of
the reaction. The final pressure at ambient temperature was 2 atm.
The reaction mixture, after cooling to room temperature, was diluted
with cold (5−10 °C) water (100 mL) and extracted with Et2O (5 ×
25 mL). The combined organic extracts were washed with H2O (3 ×
20 mL) and dried over K2CO3. Et2O was evaporated, and the crude
product (4.96 g, viscous, light brown oil) was analyzed by GC MS
(Shimadzu GCMSQP5050A spectrometer).
ASSOCIATED CONTENT
■
sı
* Supporting Information
The Supporting Information is available free of charge at
Description of Wertz52 entropy correction for dimethyl
sulfoxide, complete ref 54, details of modeling the
kinetic curves of the interaction of acetone with
acetylene in KOH/DMSO superbase, calculated imag-
inary frequencies of transition-state species, and tables of
the Cartesian coordinates and electronic energies of all
the stationary points (PDF)
To estimate the activation free energy in solution, we used an
approach that is based on results proposed by Wertz52 in ref 53. When
applied to dimethyl sulfoxide, this approach suggests that the entropy
in the DMSO solution Ssol can be obtained from the entropy for ideal
gas in the harmonic approximation (Sharm) as follows
Ssol = 0.74Sharm − 3.21 cal·mol−1·K−1
(1)
AUTHOR INFORMATION
The superbasicity was taken into account at the level of the
monosolvate model, with the KOH·DMSO complex being explicitly
included in the calculation (MONOPCM model).38 It should be noted
that in ref 38 we determined the most preferable mutual orientation of
the molecules in the monosolvate complexes of the type Nu·K+·
DMSO·H2O (where Nu stands for nucleophile), as well as in the
transition states connecting these complexes. The DMSO molecule is
coordinated by the oxygen atom to the potassium cation of Nu·K+,
■
Corresponding Author
Nadezhda M. Vitkovskaya − Laboratory of Quantum-
Chemical Modeling of Molecular Systems, Irkutsk State
University, 664003 Irkutsk, Russian Federation;
7447
J. Org. Chem. 2021, 86, 7439−7449