DFT study of the effect of sigma-ligands on the structure of ester enolates in THF, as models of the active center in the anionic polymerization of methyl methacrylate.

A Density Functional Theory (DFT) study was carried out on structures of the lithium ester enolate of methyl isobutyrate (MIB-Li) in THF solution, in the presence of TMEDA, dimethoxyethane (DME), crown ether 12-crown-4, and cryptand-2,1,1, as electron donor ligands (sigma-ligands). Both specific solvation with THF and/or ligand molecules and nonspecific solvation by the solvent continuum were taken into account. The possibility of ligand-separated ion pair formation was analyzed for each of the ligands, including THF alone. In most cases peripherally solvated dimers are the most stable species. Only in the presence of cryptand-2,1,1 was a ligand-separated triple ion pair, (MIB-Li-MIB)(-)(THF)(2),Li(2,1,1)(1)(+), shown to be comparable in stability to the THF-solvated dimer, (MIB-Li)(2)(THF)(4). These results are in agreement with experimental NMR data on the structure of MIB-Li in the presence of DME, 12-crown-4, and cryptand-2,1,1. An upfield shift of the (13)C NMR signal of the alpha-carbon of MIB-Li observed in the presence of cryptand-2,1,1, originally attributed to a ligand-separated monomer, MIB(-),Li(2,1,1)(+), was well reproduced by Hartree--Fock calculated NMR shifts for the predicted ligand-separated triple ion pair.     

Investigation of noncovalent interactions in deprotonated peptides: structural and energetic competition between aggregation and hydration

Noncovalent peptide-peptide and peptide-water interactions in small model systems were examined using an electrospray mass spectrometer equipped with a high-pressure drift cell. The results of these aggregation and hydration experiments were interpreted with the aid of molecular mechanics (MM) and density functional theory (DFT) calculations. The systems investigated include bare deprotonated monomers and dimers [P(1,2)-H](-) and hydrated deprotonated monomers and dimers [P(1,2)-H](-).(H(2)O)(n)() for the peptides dialanine (P = AA) and diglycine (P = GG). Mass spectra indicated that both peptides AA and GG form exclusively dimer ions in the electrospray process. Monomeric ions were generated by high-energy injection of the dimers into the drift cell. Temperature-dependent hydration equilibrium experiments carried out in the drift cell yielded water binding energies ranging from 11.7 (first water molecule) to 7.1 kcal/mol (fourth water) for [AA-H](-) and 11.0 to 7.4 kcal/mol for [GG-H](-). The first water molecule adding to the dimer ions [AA-H](-).(AA) and [GG-H](-).(GG) is bound by 8.4 and 7.5 kcal/mol, respectively. The hydration mass spectra for the monomers and dimers provide a means to compare the ability of water and a neutral peptide to solvate a deprotonated peptide [P-H](-). The data indicate that a similar degree of solvation is achieved by four water molecules, [P-H](-).(H(2)O)(4), or one neutral peptide, [P-H](-).(P). Temperature-dependent kinetics experiments yielded activation energies for dissociation of the dimers [AA-H](-).(AA) and [GG-H](-).(GG) of 34.9 and 32.2 kcal/mol, respectively. MM and DFT calculations carried out for the dialanine system indicated that the dimer binding energy is 24.3 kcal/mol, when the [AA-H](-) and AA products are relaxed to their global minimum structures. However, a value of 38.9 kcal/mol is obtained if [AA-H](-) and AA dissociate but retain the structures of the moieties in the dimer, suggesting the transition state occurs early in the dissociation process. Similar results were found for the diglycine dimer.     

Absolute solvation free energy of Li+ and Na+ ions in dimethyl sulfoxide solution: a theoretical ab initio and cluster-continuum model study

The solvation of the lithium and sodium ions in dimethyl sulfoxide solution was theoretically investigated using ab initio calculations coupled with the hybrid cluster-continuum model, a quasichemical theory of solvation. We have investigated clusters of ions with up to five dimethyl sulfoxide (DMSO) molecules, and the bulk solvent was described by a dielectric continuum model. Our results show that the lithium and sodium ions have four and five DMSO molecules into the first coordination shell, and the calculated solvation free energies are -135.5 and -108.6 kcal mol(-1), respectively. These data suggest a solvation free energy value of -273.2 kcal mol(-1) for the proton in dimethyl sulfoxide solution, a value that is more negative than the present uncertain experimental value. This and previous studies on the solvation of ions in water solution indicate that the tetraphenylarsonium tetraphenylborate assumption is flawed and the absolute value of the free energy of transfer of ions from water to DMSO solution is higher than the present experimental values.     

Modeling protic to dipolar aprotic solvent rate acceleration and leaving group effects in S(N)2 reactions: A theoretical study of the reaction of acetate ion with ethyl halides in aqueous and dimethyl sulfoxide solutions

The S(N)2 reactions between acetate ions and ethyl chloride, ethyl bromide, and ethyl iodide in aqueous and dimethyl sulfoxide (DMSO) solutions were theoretically investigated at an ab initio second-order M?ller-Plesset perturbation level of theory for geometry optimizations and at a fourth-order M?ller-Plesset perturbation level for energy calculations. The solvent effect was included by the polarizable continuum model using the Pliego and Riveros parametrization for DMSO and the Luque et al. scale factor for the water solution. The calculated DeltaG() values of 24.9, 20.0, and 18.5 kcal mol(-1) in a DMSO solution for ethyl chloride, ethyl bromide, and ethyl iodide are in good agreement with the estimated experimental values of 22.3, 20.0, and 16.6 kcal mol(-1), respectively. In an aqueous solution, the theoretical Delta G++ barriers of 26.9, 23.1, and 22.1 kcal mol(-1) are also in good agreement with the estimated experimental values of 26.1, 25.2, and 24.7 kcal mol(-1), respectively. The present ab initio calculations are reliable to predict the absolute and relative reactivities of ethyl halides in a DMSO solution, but in the aqueous phase, the results are less accurate. The protic to dipolar aprotic solvent rate acceleration is theoretically predicted, although this effect is underestimated. We suggest that further improvement of the present results could be obtained by including liquid-phase optimization in both solvents and treating specific solvation by water molecules for the reaction in the aqueous phase.     

Role of solvation in the reduction of the uranyl dication by water: a density functional study

We have studied the solvation of uranyl, UO(2)(2+), and the reduced species UO(OH)(2+) and U(OH)(2)(2+) systematically using three levels of approximation: direct application of a continuum model (M1); explicit quantum-chemical treatment of the first hydration sphere (M2); a combined quantum-chemical/continuum model approach (M3). We have optimized complexes with varying numbers of aquo ligands (n = 4-6) and compared their free energies of solvation. Models M1 and M2 have been found to recover the solvation energy only partially, underestimating it by approximately 100 kcal/mol or more. With our best model M3, the calculated hydration free energy Delta(h)G degrees of UO(2)(2+) is about -420 kcal/mol, which shifts to about -370 kcal/mol when corrected for the expected error of the model. This value agrees well with the experimentally determined interval, -437 kcal/mol < Delta(h)G degrees < -318 kcal/mol. Complexes with 5 and 6 aquo ligands have been found to be about equally favored with models M2 and M3. The same solvation models have been applied to a two-step reduction of UO(2)(2+) by water, previously theoretically studied in the gas phase. Our results show that the solvation contribution to the reaction free energy, about 60 kcal/mol, dominates the endoergicity of the reduction.     

Thermochemistry of arylselanyl radicals and the pertinent ions in acetonitrile

Reduction and oxidation potentials of a series of parasubstituted phenylselanyl radicals, XC(6)H(4)Se(*), have been measured using photomodulated voltammetry in acetonitrile. The thermodynamic significance of these data was substantiated through a study of the oxidation process of the pertinent selenolates in linear sweep voltammetry. Both the reduction and the oxidation potentials correlate linearly with the Hammett substituent coefficients sigma and sigma(+) leading in the latter case to slopes, rho(+), of 2.5 and 3.8, respectively. Through comparison of these slopes with those published previously for the O- and S-centered analogues, it is revealed that the pi-interaction becomes progressively smaller as the size of the radical center increases in the order O, S, and Se. Solvation energies of the pertinent selenolates and selanylium ions have been extracted from thermochemical cycles incorporating the measured electrode potentials for XC(6)H(4)Se(*) as well as electron affinities and ionization potentials obtained from theoretical calculations at the B3LYP/6-31+G(d) level. The extracted data show the expected overall substituent dependency for both kinds of ions; that is, the absolute value of the solvation energy decreases as the charge becomes more delocalized. The data have also been compared with solvation energies computed using the polarizable continuum model (PCM). Interestingly, we find that, while the model seems to work well for selenolates, it underestimates the solvation of selanylium ions in acetonitrile by as much as 25 kcal mol(-)(1). These large deviations are ascribed to the fact that the PCM method does not take specific solvent effects into account as it treats the solvent as a continuum described solely by its dielectric constant. Gas-phase calculations show that the arylselanylium ions can coordinate covalently to one or two molecules of acetonitrile in strong Ritter-type adducts. When this strong interaction is included in the solvation energy calculations by means of a combined supermolecule and PCM approach, the experimental data are reproduced within a few kcal mol(-)(1). Although the energy difference of the singlet and triplet spin states of the arylselanylium ions is small for the gas-phase structures, the singlet cation is undoubtedly the dominating species in solution because the triplet cation lacks the ability to form covalent bonds.     

Probing the effects of heterogeneity on delocalized pi...pi interaction energies

Dimers composed of benzene (Bz), 1,3,5-triazine (Tz), cyanogen (Cy) and diacetylene (Di) are used to examine the effects of heterogeneity at the molecular level and at the cluster level on pi...pi stacking energies. The MP2 complete basis set (CBS) limits for the interaction energies (E(int)) of these model systems were determined with extrapolation techniques designed for correlation consistent basis sets. CCSD(T) calculations were used to correct for higher-order correlation effects (deltaE(CCSD)(T)(MP2)) which were as large as +2.81 kcal mol(-1). The introduction of nitrogen atoms into the parallel-slipped dimers of the aforementioned molecules causes significant changes to E(int). The CCSD(T)/CBS E(int) for Di-Cy is -2.47 kcal mol(-1) which is substantially larger than either Cy-Cy (-1.69 kcal mol(-1)) or Di-Di (-1.42 kcal mol(-1)). Similarly, the heteroaromatic Bz-Tz dimer has an E(int) of -3.75 kcal mol(-1) which is much larger than either Tz-Tz (-3.03 kcal mol(-1)) or Bz-Bz (-2.78 kcal mol(-1)). Symmetry-adapted perturbation theory calculations reveal a correlation between the electrostatic component of E(int) and the large increase in the interaction energy for the mixed dimers. However, all components (exchange, induction, dispersion) must be considered to rationalize the observed trend. Another significant conclusion of this work is that basis-set superposition error has a negligible impact on the popular deltaE(CCSD)(T)(MP2) correction, which indicates that counterpoise corrections are not necessary when computing higher-order correlation effects on E(int). Spin-component-scaled MP2 (SCS-MP2 and SCSN-MP2) calculations with a correlation-consistent triple-zeta basis set reproduce the trends in the interaction energies despite overestimating the CCSD(T)/CBS E(int) of Bz-Tz by 20-30%.     

Gas-Phase and Solution-Phase Homolytic Bond Dissociation Energies of H-N(+) Bonds in the Conjugate Acids of Nitrogen Bases

The oxidation potentials of 19 nitrogen bases (abbreviated as B: six primary amines, five secondary amines, two tertiary amines, three anilines, pyridine, quinuclidine, and 1,4-diazabicyclo[2,2,2]octane), i.e., E(ox)(B) values in dimethyl sulfoxide (DMSO) and/or acetonitrile (AN), have been measured. Combination of these E(ox)(B) values with the acidity values of the corresponding acids (pK(HB)(+)) in DMSO and/or AN using the equation: BDE(HB)(+) = 1.37pK(HB)(+) + 23.1 E(ox)(B) + C (C equals 59.5 kcal/mol in AN and 73.3 kcal/mol in DMSO) gave estimates of solution phase homolytic bond dissociation energies of H-B(+) bonds. Gas-phase BDE values of H-B(+) bonds were estimated from updated proton affinities (PA) and adiabatic ionization potentials (aIP) using the equation, BDE(HB(+))(g) = PA + aIP - 314 kcal/mol. The BDE(HB)(+) values estimated in AN were found to be 5-11 kcal/mol higher than the corresponding gas phase BDE(HB(+))(g) values. These bond-strengthening effects in solution are interpreted as being due to the greater solvation energy of the HB(+) cation than that of the B(+*) radical cation.     

Ab initio study of substituent effects in the interactions of dimethyl ether with aromatic rings

Ab initio calculations have been used to investigate the interaction energies of dimers of dimethyl ether with benzene, hexafluorobenzene, and several monosubstituted benzenes. The potential energy curves were explored at the MP2/aug-cc-pVDZ level for two basic configurations of the dimers, one in which the oxygen atom of the dimethyl ether was pointed towards the aromatic ring and the other in which the oxygen atom was pointed away from the aromatic ring. Once the optimum intermolecular distances between the dimethyl and the aromatic ring had been determined for each of the dimers in both configurations at the MP2/aug-cc-pVDZ level, single point energy calculations were performed at the MP2/aug-cc-pVTZ level. A CCSD(T) correction term to the energy was determined and this was combined with the MP2/aug-cc-pVTZ energies to estimate the CCSD(T)/aug-cc-pVTZ interaction energies of the dimers. The estimated CCSD(T)/aug-cc-pVTZ interaction energies are predicted to be attractive for all of the dimers in both configurations and dispersion interactions are found to be a large component of the stabilization of the dimers. For the dimers with the dimethyl ether oxygen pointing towards the aromatic ring, the strengths of interaction energies are found to increase as the aromatic ring becomes more electron deficient, while for the dimers with the dimethyl ether oxygen pointing away from the aromatic ring, they increase as the aromatic ring becomes more electron rich. In both cases, the trends can be explained in terms of the electrostatic potentials of the dimethyl ether and the aromatic rings.     

Ab initio studies on the mechanism of the size-dependent hydrogen-loss reaction in Mg+(H2O)n

The mechanism of size-dependent intracluster hydrogen loss in the cluster ions Mg(+)(H(2)O)(n), which is switched on around n=6, and off around n=14, was studied by ab initio calculations at the MP2/6-31G* and MP2/6-31G** levels for n=1-6. The reaction proceeds by Mg(+)-assisted breaking of an H-O bond in one of the H(2)O molecules. The reaction barrier is dependent on both the cluster size and the solvation structure. As n increases from 1 to 6, there is a dramatic drop in the reaction barrier, from greater than 70 kcal mol(-1) for n=1 to less than 10 kcal mol(-1) for n=6. In the transition structures, the Mg atom is close to the oxidation state of +2, and H(2)O molecules in the first solvation shell are much more effective in stabilizing the transition structures and lowering the reaction barriers than H(2)O molecules in the other solvation shells. While the reaction barrier for trimer core structures with only three H(2)O molecules in the first shell is greater than 24 kcal mol(-1), even for Mg(+)(H(2)O)(6), it drops considerably for clusters with four-six H(2)O molecules in the first shell. The more highly coordinated complexes have comparable or slightly higher energy than the trimer core structures, and the presence of such high coordination number complexes is the underlying kinetic factor for the switching on of the hydrogen-loss reaction around n=6. For clusters with trimer core structures, the hydrogen loss reaction is much easier when it is preceded by an isomerization step that increases the coordination number around Mg(+). Delocalization of the electron on the singly occupied molecular orbital (SOMO) away from the Mg(+) ion is observed for the hexamer core structure, while at the same time this isomer is the most reactive for the hydrogen-loss reaction, with an energy barrier of only 2.7 kcal mol(-1) at the MP2/6-31G** level.     

Combining the polarizable Drude force field with a continuum electrostatic Poisson–Boltzmann implicit solvation model

In this work, we have combined the polarizable force field based on the classical Drude oscillator with a continuum Poisson–Boltzmann/solvent‐accessible surface area (PB/SASA) model. In practice, the positions of the Drude particles experiencing the solvent reaction field arising from the fixed charges and induced polarization of the solute must be optimized in a self‐consistent manner. Here, we parameterized the model to reproduce experimental solvation free energies of a set of small molecules. The model reproduces well‐experimental solvation free energies of 70 molecules, yielding a root mean square difference of 0.8 kcal/mol versus 2.5 kcal/mol for the CHARMM36 additive force field. The polarization work associated with the solute transfer from the gas‐phase to the polar solvent, a term neglected in the framework of additive force fields, was found to make a large contribution to the total solvation free energy, comparable to the polar solute–solvent solvation contribution. The Drude PB/SASA also reproduces well the electronic polarization from the explicit solvent simulations of a small protein, BPTI. Model validation was based on comparisons with the experimental relative binding free energies of 371 single alanine mutations. With the Drude PB/SASA model the root mean square deviation between the predicted and experimental relative binding free energies is 3.35 kcal/mol, lower than 5.11 kcal/mol computed with the CHARMM36 additive force field. Overall, the results indicate that the main limitation of the Drude PB/SASA model is the inability of the SASA term to accurately capture non‐polar solvation effects. © 2018 Wiley Periodicals, Inc.     

Noncovalent complexes between dimethyl ether and formic acid--an ab initio and matrix isolation study.

The complexes formed by noncovalent interactions between formic acid and dimethyl ether are investigated by ab initio methods and characterized by matrix isolation spectroscopy. Six complexes with binding energies between -2.26 and -7.97 kcal mol(-1) (MP2/cc-pVTZ+zero point vibrational energy+basis set superposition erros) are identified. The two strongest bound complexes are, within a range of 0.3 kcal mol(-1), isoenergetic. The binding in these six dimers can be described in terms of OH...O, C=O...H, C-O...H and CH...O interactions. Matrix isolation spectroscopy allowed to characterize the two strongest bound complexes by their infrared spectra.     

Structure and stability of the [TCNE](2)(2-) dimers in dichloromethane solution: a computational study

The solution behavior of [TCNE](.-), which forms long-living pi-[TCNE]22- dimers, is computationally studied by B3LYP and MCQDPT/CASSCF(2,2) calculations (a multiconfigurational quasi-degenerate perturbative calculation using a CASSCF(2,2) wavefunction, which properly accounts for the dispersion interaction). B3LYP calculations indicate minimum-energy [TCNE](2)(2-)(dichloromethane)(4) aggregates, a solvent where pi-[TCNE](2)(2-) dimers are spectroscopically observed. Their existence is attributed to [TCNE](.-)...solvent interactions that exceed the [TCNE](.-)...[TCNE](.-) repulsion. The lowest energy minimum at the B3LYP level corresponds to an open-shell singlet electronic structure, a metastable minimum where the shortest interanion C...C distance is 5.23 A. A slightly less stable minimum is also found for the closed-shell singlet when double-occupancy of the orbitals is imposed, but it converts into the open-shell singlet minimum when the double occupancy is relaxed. At the MCQDPT/CASSCF(2,2) level, the only minimum is for the closed-shell singlet (24.0 kcal/mol (101 kJ/mol) more stable than the dissociation products), consistent with experimental enthalpy of dimerization of [TCNE](.-) in dichloromethane solutions. It has an interanion C...C distance of 2.75 A and is in accord with the UV-vis experimental properties of the [TCNE](.-) solutions.     

Synthesis,structure, and dynamic behavior of cyclopentadienyl-lithium, -sodium,and -potassium annelated with bicyclo[2.2.2]octene units: a systematic study on site exchange of alkali metals on a cyclopentadienyl ring in tetrahydrofuran

Novel cyclopentadienyl (Cp)-alkali metal complexes 1-M and 2-M (M = Li, Na, K), in which the Cp ring is annelated with two bicyclo[2.2.2]octene units and substituted with a phenyl group for 1 and a tert-butyl group for 2, were synthesized, and their structures and dynamic behaviors were investigated by means of X-ray crystallography, dynamic (13)C NMR, and DFT calculations. The X-ray crystallography results indicated that 1-Li, 1-Na, and 2-Na form monomeric contact ion pairs (CIP) with three THF molecules coordinated to the metal atom. Also, in THF-d(8), all of the 1-M and 2-M form monomeric CIP in the ground state. However, variable-temperature (13)C NMR measurements of 1-M and 2-M in THF-d(8) demonstrated dynamic behavior in which the metal ion exchanges positions between the upper and lower faces of the Cp ring. From a study of the concentration dependence of the dynamic behavior, the exchange was found to proceed principally as an intramolecular process at concentration ranges lower than 0.2 M. The experimentally observed deltaG values for the intramolecular exchange process for all the 1-M and 2-M (except for 2-Li, whose intramolecular process was too slow to observe) were found to be quite similar in THF-d(8) solution and to fall within the range of 12-14 kcal mol(-)(1). Within this range, a tendency was observed for the deltaG values to increase as the size of the metal decreased. Theoretical calculations (B3LYP/6-31G(d)) afforded considerably large values as the gas-phase dissociation energy for 1-M (162.7 kcal mol(-)(1) for M = Li; 131.6 kcal mol(-)(1) for M = Na; 110.9 kcal mol(-)(1) for M = K) and for 2-M (170.0 kcal mol(-)(1) for M = Li; 137.5 kcal mol(-)(1) for M = Na; 115.4 kcal mol(-)(1) for M = K). These values should be compensated for by a decrease in the solvation energies for the metal ions with increasing size, as exemplified by the calculated solvation energy for M(+)(Me(2)O)(4), which serves as a model for metal ions solvated with four molecules of THF (-122.9 kcal mol(-)(1) for M = Li; -94.7 kcal mol(-)(1) for M = Na; -67.7 kcal mol(-)(1) for M = K). This compensation results in a small difference in the overall energy for dissociation of 1-M or 2-M in ethereal solutions, thus supporting the similar deltaG values observed for the intramolecular metal exchange.     

Methylene transfer or carbometalation? A theoretical study to determine the mechanism of lithium carbenoid-promoted cyclopropanation reactions in aggregation and solvation States

Density functional theory calculations for the lithium carbenoid-promoted cyclopropanations in aggregation and solvation states are presented in order to investigate the controversy of the mechanistic dichotomy, that is, the methylene-transfer mechanism and the carbometalation mechanism. The methylene-transfer mechanism represents the reaction reality, whereas the carbometalation pathway does not appear to compete significantly with the methylene-transfer pathway and should be ruled out as a major factor. A simple model calculation for monomeric lithium carbenoid-promoted cyclopropanations with ethylene in the gas phase is not sufficient to reflect the reaction conditions accurately or to determine the reaction mechanism since its result is inconsistent with the experimental facts. The aggregated lithium carbenoids are the most probable reactive species in the reaction system. The calculated reaction barriers of the methylene-transfer pathways are 10.1 and 8.0 kcal/mol for the dimeric (LiCH2F)2 and tetrameric (LiCH2F)4 species, respectively, compared with the reaction barrier of 16.0 kcal/mol for the monomeric LiCH2F species. In contrast, the reaction barriers of the carbometalation pathways are 26.8 kcal/mol for the dimeric (LiCH2F)2 and 33.9 kcal/mol for the tetrameric (LiCH2F)4 species, compared with the reaction barrier of 12.5 kcal/mol for the monomeric LiCH2F species. The effects of solvation were investigated by explicit coordination of the solvent molecules to the lithium centers. This solvation effect is found to enhance the methylene-transfer pathway, while it is found to impede the carbometalation pathway instead. The combined effects of the aggregation and solvation lead to barriers to reaction in the range of 7.2-9.0 kcal/mol for lithium carbenoid-promoted cyclopropanation reactions along the methylene-transfer pathway. Our computational results are in good agreement with the experimental observations.     

Theoretical comparison of ketene dimerization in the gas and liquid phase

We present the first theoretical comparison between ketene dimerization in gas phase and ketene dimerization in solution. Density functional theory (DFT) calculations on the ketene dimerization were carried out considering the following product dimers: diketene (d-I), 1,3-cyclobutanedione (d-II), 2,4-dimethylene-1,3-dioxetane (d-III), and 2-methyleneoxetan-3-one (d-IV). All structures were optimized at the PW86x+PBEc/DZP level of theory. Based on these geometries, a total of 58 meta and hybrid functionals were used to evaluate the heat of dimerization. The MPW1K functional was found to fit the experimental data best and subsequently used in the final analyses for all energy calculations. It was found on both kinetic and thermodynamic grounds that only d-I and d-II are formed during ketene dimerization in gas phase and solution. In gas phase, d-I is favored over d-II by 2 kcal/mol. However, the dimerization barrier for d-I is 1 kcal/mol higher than for d-II. Solvation makes dimerization more favorable. On the enthalpic surface this is due to a favorable interaction between the dimer dipole moment and solvent molecules. The dimer is stabilized further on the Gibbs energy surface by an increase of the dimerization entropy in solution compared to gas phase. The species d-I remains the most stable dimer in solution by 1 kcal/mol. Kinetically, the dimerization barriers for the relevant species d-I and d-II are cut in half by solvation, due to both favorable dimer-dipole/solvent interactions (DeltaH++, DeltaG++) and an increase in the activation entropies (DeltaS++). While the dimerization barrier for d-II is lowest for the gas phase and toluene, the barrier for d-I formation becomes lowest for the more polar solvent acetone by 1 kcal/mol as d-I dimerization has the most polar transition state.     

Accurate pK(a) calculations for carboxylic acids using complete basis set and Gaussian-n models combined with CPCM continuum solvation methods

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pK(a) values for six carboxylic acids. An experimental value of -264.61 kcal/mol for the free energy of solvation of H(+), DeltaG(s)(H(+)), was combined with a value for G(gas)(H(+)) of -6.28 kcal/mol to calculate pK(a) values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pK(a) values accurate to less than half a pK(a) unit.     

Amine-chelated aryllithium reagents--structure and dynamics

Multinuclear NMR studies of five-membered-ring amine chelated aryllithium reagents 2-lithio-N,N-dimethylbenzylamine (1), the diethylamine and diisopropylamino analogues (2, 3), and the o-methoxy analogue (4), isotopically enriched in (6)Li and (15)N, have provided a detailed picture of the solution structures in ethereal solvents (usually in mixtures of THF and dimethyl ether, ether, and 2,5-dimethyltetrahydrofuran). The effect of cosolvents such as TMEDA, PMDTA, and HMPA has also been determined. All compounds are strongly chelated, and the chelation is not disrupted by these cosolvents. Reagents 1, 2, and 3 are dimeric in solvents containing a large fraction of THF. Below -120 degrees C, three chelation isomers of the dimers are detectable by NMR spectroscopy: one (A) with both nitrogens coordinated to one lithium of the dimer, and two (B and C) in which each lithium bears one chelating group. Dynamic NMR studies have provided rates and activation energies for the interconversion of the 1-A, 1-B, and 1-C isomers. They interconvert either by simple ring rotation, which interconverts B and C, or by amine decoordination (probably associative, DeltaG(++)(-93) = 8.5 kcal/mol), which can interconvert all of the isomers. The dimers of 1 are thermodynamically more stable than those of model systems such as phenyllithium, o-tolyllithium, or 2-isoamylphenyllithium (5, DeltaDeltaG > or = 3.3 kcal/mol). They are not detectably deaggregated by TMEDA or PMDTA, although HMPA causes partial deaggregation. The dimers are also more robust kinetically with rates of interaggregate exchange, measured by DNMR line shape analysis of the C-Li signal, orders of magnitude smaller than those of models (DeltaDeltaG(++) > or = 4.4 kcal/mol). Similarly, the mixed dimer of 1 and phenyllithium, 13, is kinetically more stable than the phenyllithium dimer by >2.2 kcal/mol. X-ray crystal structures of the TMEDA solvate of 1-A and the THF solvate of 3-B showed them to be dimeric and chelated in the solid state as well. Compound 4, which has a methoxy group ortho to the C-Li group, differs from the others in being only partially dimeric in THF, presumably for steric reasons. This compound is fully deaggregated by 1 equiv of HMPA. Excess HMPA leads to the formation of ca. 15% of a triple ion (4-T) in which both nitrogens appear to be chelated to the central lithium.     

Equilibrium acidities and homolytic bond dissociation energies of acidic C-H bonds in alpha-arylacetophenones and related compounds

The equilibrium acidities (pK(AH)s) and the oxidation potentials of the congugate anions [E(ox)(A(-))s] were determined in dimethyl sulfoxide (DMSO) for eight ketones of the structure GCOCH(3) and 20 of the structure RCOCH(2)G, (where R = alkyl, phenyl and G = alkyl, aryl). The homolytic bond dissociation energies (BDEs) for the acidic C-H bonds of the ketones were estimated using the equation BDE(AH) = 1.37pK(AH) + 23.1E(ox)(A(-)) + 73.3. While the equilibrium acidities of GCOCH(3) were found to be dependent on the remote substituent G, the BDE values for the C-H bonds remained essentially invariant (93.5 +/- 0.5 kcal/mol). A linear correlation between pK(AH) values and [E(ox)(A(-))s] was found for the ketones. For RCOCH(2)G ketones, both pK(AH) and BDE values for the adjacent C-H bonds are sensitive to the nature of the substituent G. However, the steric bulk of the aryl group tends to exert a leveling effect on BDEs. The BDE of alpha-9-anthracenylacetophenone is higher than that of alpha-2-anthracenylacetophenone by 3 kcal/mol, reflecting significant steric inhibition of resonance in the 9-substituted system. A range of 80.7-84.4 kcal/mol is observed for RCOCH(2)G ketones. The results are discussed in terms of solvation, steric, and resonance effects. Ab initio density functional theory (DFT) calculations are employed to illustrate the effect of steric interactions on radical and anion geometries. The DFT results parallel the trends in the experimental BDEs of alpha-arylacetophenones.