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.     

Equilibrium Acidities and Homolytic Bond Dissociation Enthalpies of the Acidic C-H Bonds in P-(Para-substituted benzyl)triphenylphosphonium Cations and Related Cations

Equilibrium acidities (pK(HA)) of six P-(para-substituted benzyl)triphenylphosphonium (p-GC(6)H(4)CH(2)PPh(3)(+)) cations, P-allyltriphenylphosphonium cation, P-cinnamyltriphenylphosphonium cation, and As-(p-cyanobenzyl)triphenylarsonium cation, together with the oxidation potentials [E(ox)(A(-))] of their conjugate anions (ylides) have been measured in dimethyl sulfoxide (DMSO) solution. The acidifying effects of the alpha-triphenylphosphonium groups on the acidic C-H bonds in toluene and propene were found to be ca 25 pK(HA) units (34 kcal/mol). Introduction of an electron-withdrawing group such as 4-NO(2), 4-CN, or 4-Br into the para position of the benzyl ring in p-GC(6)H(4)CH(2)PPh(3)(+) cations resulted in an additional acidity increase, but introduction of the 4-OEt electron-donating group decreases the acidity. The equilibrium acidities of p-GC(6)H(4)CH(2)PPh(3)(+) cations were nicely linearly correlated with the Hammett sigma(-) constants of the substituents (G) with a slope of 4.78 pK(HA) units (R(2) = 0.992) (Figure 1). Reversible oxidation potentials of the P-(para-substituted benzyl)triphenylphosphonium ylides were obtained by fast scan cyclic voltammetry. The homolytic bond dissociation enthalpies (BDEs) of the acidic C-H bonds in these cations, estimated by combining their equilibrium acidities with the oxidation potentials of their corresponding conjugate anions, showed that the alpha-Ph(3)P(+) groups have negligible stabilizing or destabilizing effects on the adjacent radicals. The equilibrium acidity of As-(p-cyanobenzyl)triphenylarsonium cation is 4 pK(HA) units weaker than that of P-(p-cyanobenzyl)triphenylphosphonium cation, but the BDE of the acidic C-H bond in As-(p-cyanobenzyl)triphenylarsonium cation is ca 2 kcal/mol higher than that in P-(p-cyanobenzyl)triphenylphosphonium cation.     

Calculation of the relative acidities and oxidation potentials of para-substituted phenols. A model for alpha-tocopherol in solution

Relative acidities (Delta pK(a)) of phenols and oxidation potentials (Delta E(ox)) of the phenoxide anions have been calculated for nine para-substituted phenols using density functional theory. Solvent effects were incorporated using the conductor-like polarisable continuum method. Using the calculated Delta pK(a) and Delta E(ox) values in a thermodynamic cycle, the DeltaBDE (bond dissociation enthalpy) of the phenols were also determined with all values calculated to within 1.5 kcal mol(-1) of experiment. The Delta pK(a) and Delta E(ox) values were calculated for 6-hydroxy-2,2,5,7,8-pentamethylchroman (HPMC), a model for alpha-tocopherol for which there are no known experimental values. The acidity of this compound is raised by 2.4 pK(a) units and lowered by -0.79 V relative to phenol with a calculated Delta BDE of -14.9 kcal mol(-1). There is a negative correlation (r(2) = 0.86) between the Delta pK(a) and the Delta BDE values. A stronger and positive correlation is found between the Delta E(ox) (r(2) = 0.98) and the Delta BDE values. Using these correlations it is uncovered that hydrogen abstraction of phenols, as measured by the Delta BDE, is driven by electron transfer rather than by proton transfer.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

Quantum-chemical predictions of redox potentials of organic anions in dimethyl sulfoxide and reevaluation of bond dissociation enthalpies measured by the electrochemical methods

A first-principle theoretical protocol was developed that could predict the absolute pK(a) values of over 250 structurally unrelated compounds in DMSO with a precision of 1.4 pK(a) units. On this basis we developed the first theoretical protocol that could predict the standard redox potentials of over 250 structurally unrelated organic anions in DMSO with a precision of 0.11 V. Using the two new protocols we systematically reevaluated the bond dissociation enthalpies (BDEs) measured previously by the electrochemical methods. It was confirmed that for most compounds the empirical equation (BDE = 1.37 pK(HA) + 23.1E(o) + constant) was valid. The constant in this equation was determined to be 74.0 kcal/mol, compared to 73.3 kcal/mol previously reported. Nevertheless, for a few compounds the empirical equation could not be used because the solvation energy changed dramatically during the bond cleavage, which resulted from the extraordinary change of dipole moment during the reaction. In addition, we found 40 compounds (mostly oximes and amides) for which the experimental values were questionable by over 5 kcal/mol. Further analyses revealed that all these questionable BDEs could be explained by one of the three following reasons: (1) the experimental pK(a) value is questionable; (2) the experimental redox potential is questionable; (3) the solvent effect cannot be neglected. Thus, by developing practical theoretical methods and utilizing them to solve realistic problems, we hope to demonstrate that ab initio theoretical methods can now be developed to make not only reliable, but also useful, predictions for solution-phase organic chemistry.     

Model studies of the histidine-tyrosine cross-link in cytochrome C oxidase reveal the flexible substituent effect of the imidazole moiety

[reaction: see text] Experimental and theoretical studies were carried out to interrogate the effect of an imidazole substituent in each of the ortho, meta, and para positions on the pK(a), E degrees , and O-H BDE of phenol. The results reveal that imidazole substitution lowers the pK(a) of phenol and increases the E degrees of phenoxide due to its sigma-electron withdrawing ability (sigma(p)(-) = +0.21, sigma(m)(-) = +0.45) but decreases the O-H BDE and E degrees of phenol due to its pi-electron-donating ability (sigma(p)(+) = -0.45).     

Kinetics and thermodynamics of H. transfer from (eta5-C5R5)Cr(CO)3H (R = Ph,Me, H) to methyl methacrylate and styrene

The rates of H/D exchange have been measured between (a) the activated olefins methyl methacrylate-d(5) and styrene-d(8), and (b) the Cr hydrides (eta(5)-C(5)Ph(5))Cr(CO)(3)H (2a), (eta(5)-C(5)Me(5))Cr(CO)(3)H (2b), and (eta(5)-C(5)H(5))Cr(CO)(3)H (2c). With a large excess of the deuterated olefin the first exchange goes to completion before subsequent exchanges begin, at a rate first order in olefin and in hydride. (Hydrogenation is insignificant except with styrene and CpCr(CO)(3)H; in most cases, the radicals arising from the first H. transfer are too hindered to abstract another H. .) Statistical corrections give the rate constants k(reinit) for H. transfer to the olefin from the hydride. With MMA, k(reinit) decreases substantially as the steric bulk of the hydride increases; with styrene, the steric bulk of the hydride has little effect. At longer times, the reaction of MMA or styrene with 2a gives the corresponding metalloradical 1a as termination depletes the concentration of the methyl isobutyryl radical 3 or the alpha-methylbenzyl radical 4; computer simulation of [1a] as f(t) gives an estimate of k(tr), the rate constant for H. transfer from 3 or 4 back to Cr. These rate constants imply a DeltaG (50 degrees C) of +11 kcal/mol for H. transfer from 2a to MMA, and a DeltaG (50 degrees C) of +10 kcal/mol for H. transfer from 2a to styrene. The CH(3)CN pK(a) of 2a, 11.7, implies a BDE for its Cr-H bond of 59.6 kcal/mol, and DFT calculations give 58.2 kcal/mol for the Cr-H bond in 2c. In combination the kinetic DeltaG values, the experimental BDE for 2a, and the calculated DeltaS values for H. transfer imply a C-H BDE of 45.6 kcal/mol for the methyl isobutyryl radical 3 (close to the DFT-calculated 49.5 kcal/mol), and a C-H BDE of 47.9 kcal/mol for the alpha-methylbenzyl radical 4 (close to the DFT-calculated 49.9 kcal/mol). A solvent cage model suggests 46.1 kcal/mol as the C-H BDE for the chain-carrying radical in MMA polymerization.     

Solvolysis of para-substituted cumyl chlorides. Brown and Okamoto's electrophilic substituent constants revisited using continuum solvent models

Brown and Okamoto (J. Am. Chem. Soc. 1958, 80, 4979) derived their electrophilic substitutent constants, sigma(p)+, from the relative rates of solvolysis of ring-substituted cumyl chlorides in an acetone/water solvent mixture. Application of the Hammett equation to the rates for the meta-substituted cumyl chlorides, where there could be no resonance interaction with the developing carbocation, gave a slope, rho(+) = -4.54 ( identical with 6.2 kcal/mol free energy). Rates for the para-substituted chlorides were then used to obtain sigma(p)+ values. We have calculated gas-phase C-Cl heterolytic bond dissociation enthalpy differences, Delta BDE(het) (= BDE(het)(4-YC(6)H(4)CMe(2)Cl) - BDE(het)(C(6)H(5)CMe(2)Cl)), for 16 of the 4-Y substituents employed by Brown and Okamoto. The plot of Delta BDE(het) vs sigma(p)+ gave rho(+) (SD) = 16.3 (2.3) kcal/mol, i.e., a rho(+) value roughly 2.5 times greater than experiment. Inclusion of solvation (water) energies, calculated using three continuum solvent models, reduced rho(+) and SD. The computationally least expensive model used, SM5.42R (Li et al. Theor. Chem. Acc. 1999, 103, 9) gave the best agreement with experiment. This model yielded rho(+) (SD) = 7.7 (0.9) kcal/mol, i.e., a rho(+) value that is only 24% larger than experiment.     

Water effect on the o-h dissociation enthalpy of para-substituted phenols: a DFT study

The effect of water on the O-H bond dissociation enthalpy (BDE) of para-substituted phenols has been investigated by means of DFT calculations. It is shown that the experimental BDE values are fairly well-reproduced by simple B3LYP/6-31G* calculations carried out on the phenol/phenoxyl-water complexes taking into account only hydrogen-bonding (HB) interactions of water molecules with molecular sites (HB model). On the contrary, the BDE values computed with the polarizable continuum model (PCM/B3LYP/6-31G*)8 are overestimated by about 3-4 kcal/mol. Discrepancy between theory and experiment increases using the PCM method in addition to the HB model. Calculations show that, in general, the HB interaction with water molecules decreases the BDE of phenols bearing electron-releasing groups while increasing the BDE of phenols bearing electron-withdrawing substituents. This opposite effect is explained by considering the resonance structures with charge separation both in phenols and in phenoxyl radicals. With electron donors, the phenoxyl radical is preferentially stabilized by the HB acceptor interaction with two water molecules, while with electron acceptors the phenol is preferentially stabilized by the HB donor interaction with one water molecule.     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

Successive attachment of electrons to protonated Guanine: (G+H)* radicals and (G+H)- anions

The structures, energetics, and vibrational frequencies of nine hydrogenated 9H-keto-guanine radicals (G+H)(*) and closed-shell anions (G+H)(-) are predicted using the carefully calibrated (Chem. Rev. 2002, 102, 231) B3LYP density functional method in conjunction with a DZP++ basis set. These radical and anionic species come from consecutive electron attachment to the corresponding protonated (G+H)(+) cations in low pH environments. The (G+H)(+) cations are studied using the same level of theory. The proton affinity (PA) of guanine computed in this research (228.1 kcal/mol) is within 0.7 kcal/mol of the latest experiment value. The radicals range over 41 kcal/mol in relative energy, with radical r1, in which H is attached at the C8 site of guanine, having the lowest energy. The lowest energy anion is a2, derived by hydride ion attachment at the C2 site of guanine. No stable N2-site hydride should exist in the gas phase. Structure a9 was predicted to be dissociative in this research. The theoretical adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies were computed, with AEAs ranging from 0.07 to 3.12 eV for the nine radicals.     

Acidity of anilines: calculations vs experiment

The gas-phase acidities of ca. 60 monosubstituted anilines (with acidity span of almost 50 kcal mol(-1)) have been calculated using density functional theory (DFT) at the B3LYP/6-311+G** level. At this relatively simple level of theory the calculated (ΔG(calc)) and available experimental (ΔG(exp)) acidities are in reasonable quantitative correlation according to the following equation: ΔG(obs) = a + bΔG(calc), where a=20.79, b=0.942, n=27, R(2)=0.990, and s=0.78 kcal·mol(-1). The slope is not far from its ideal value. Substituent effects on the acidities were dissected separately into those operating in the neutral acid molecule and in its conjugated anion using the isodesmic homodesmotic reactions. All in all, both forms, neutral and anionic, are contributing in combination to make up the gross acidity of anilines. However, the contributions of the anions into the gross substituent effects are much larger than the substituent effects in the neutral anilines. Some of the systems were used in testing a relatively new theoretical model, COSMO-RS (conductor-like screening model for real solvents), using it for the prediction of pK(a) values in DMSO. The method proved to be rather accurate for showing pK(a) trends (R(2)=0.980 in DMSO). However, the predicted absolute pK(a) values were all somewhat lower (rmsd=2.49 kcal·mol(-1)) than the respective experimental values.     

Methyl cation affinities of rare gases and nitrogen and the heat of formation of diazomethane

The methyl cation affinities of the rare gases have been calculated at 0 and 298 K by using coupled cluster theory including noniterative, quasiperturbative triple excitations with the new correlation-consistent basis sets for Xe up through aug-cc-pV5Z in some cases. To achieve near chemical accuracy (+/-1 kcal/mol) in the thermodynamic properties, we add to the estimated complete basis set valence binding energies, based on frozen core coupled cluster theory energies, two corrections: (1) a core/valence correction and (2) a scalar relativistic correction. Vibrational zero-point energies were computed at the coupled cluster level of theory at the CCSD(T)/aug-cc-pVDZ level. The calculated rare gas methyl cation affinities (MCA in kcal/mol) at 298 K are the following: MCA(He) = 1.7, MCA(Ne) = 2.5, MCA(Ar) = 16.9, MCA(Kr) = 25.5, and MCA(Xe) = 36.6. Because of the importance of the MCA(N(2)) in the experimental measurements of the MCA scale, we calculated a number of quantities associated with CH(3)N(2)(+) and CH(2)N(2). The calculated values for diazomethane at 298 K are: DeltaH(f)(CH(2)N(2)) = 65.3 kcal/mol, PA(CH(2)N(2)) = 211.9 kcal/mol, and MCA(N(2)) = 43.2 kcal/mol.     

Acidity and alkali metal adsorption on the SiO2-aqueous solution interface

The intrinsic deprotonation constant (pK(a(2))(int)) and the intrinsic ion exchange constants (pK(Me(+))(int)) of Li(+), Na(+), and K(+) on SiO(2) were uniquely determined at 30 degrees C by using the potentiometric titration data, the Gouy-Chapman-Stern-Grahame (GSCG) model for the structure of the electrical double-layer (edl) and the double-extrapolation method. The values of these constants were pK(a(2))(int) = 6.57, pK(Li(+))(int) = pK(Na(+))(int) = pK(K(+))(int) = 5.61. The chemical meaning of intrinsic equilibrium constants and the equality in the values of pK(Li(+))(int), pK(Na(+))(int) and pK(K(+))(int) were discussed.     

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.     

Thermochemical parameters of CHFO and CF2O

Thermochemical properties of CHFO and CF 2O and their derivatives were calculated by using coupled-cluster theory (U)CCSD(T) calculations with the aug-cc-pV nZ ( n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. The predicted properties include the following. Enthalpies of formation (298 K, kcal/mol): Delta H f (CF 2O) = -144.7, Delta H f(CHFO) = -91.1, Delta H f (CFO (*)) = -41.6. Bond dissociation energy (0 K, kcal/mol): BDE(CFO-F) = 120.7, BDE(CHO-F) = 119.1, BDE(CFO-H) = 100.2. Ionization potential (eV): IP 1(CF 2O) = 13.04, IP 2(CF 2O) = 14.09, IP 1(CHFO) = 12.41, IP 2(CHFO) = 13.99, IP 1(CFO (*)) = 9.34. Proton affinity (298 K, kcal/mol), PA O(CF 2O) = 148.8, PA O(CHFO) = 156.7, PA F(CHFO) = 154.5 kcal/mol. Electron affinity: EA(CFO (*)) = 2.38 eV. Triplet-singlet separation gap (eV): Delta E T1-S0(CF 2O) = 4.47, Delta E T1-S0(CHFO) = 4.36. Triplet-triplet transition energy (eV): Delta E T2-T1(CF 2O) = 0.44. The new calculated values contribute to solving some persistent discrepancies in the literature. The effects of F-atoms on thermochemical parameters are not linearly additive, and the changes are largely dominated by the first F-substitution. On the basis of the calculated proton affinities of CF 2O and CF 3OH, the nucleophilicities of the oxygen atoms are, within computational errors, the same in both compounds.     

Potential energy surface for activation of methane by Pt(+): a combined guided ion beam and DFT study

A guided-ion beam tandem mass spectrometer is used to study the reactions of Pt(+) with methane, PtCH(2)(+) with H(2) and D(2), and collision-induced dissociation of PtCH(4)(+) and PtCH(2)(+) with Xe. These studies experimentally probe the potential energy surface for the activation of methane by Pt(+). For the reaction of Pt(+) with methane, dehydrogenation to form PtCH(2)(+) + H(2) is exothermic, efficient, and the only process observed at low energies. PtH(+), formed in a simple C-H bond cleavage, dominates the product spectrum at high energies. The observation of a PtH(2)(+) product provides evidence that methane activation proceeds via a (H)(2)PtCH(2)(+) intermediate. Modeling of the endothermic reaction cross sections yields the 0 K bond dissociation energies in eV (kJ/mol) of D(0)(Pt(+)-H) = 2.81 +/- 0.05 (271 +/- 5), D(0)(Pt(+)-2H) = 6.00 +/- 0.12 (579 +/- 12), D(0)(Pt(+)-C) = 5.43 +/- 0.05 (524 +/- 5), D(0)(Pt(+)-CH) = 5.56 +/- 0.10 (536 +/- 10), and D(0)(Pt(+)-CH(3)) = 2.67 +/- 0.08 (258 +/- 8). D(0)(Pt(+)-CH(2)) = 4.80 +/- 0.03 eV (463 +/- 3 kJ/mol) is determined by measuring the forward and reverse reaction rates for Pt(+) + CH(4) right harpoon over left harpoon PtCH(2)(+) + H(2) at thermal energy. We find extensive hydrogen scrambling in the reaction of PtCH(2)(+) with D(2). Collision-induced dissociation (CID) of PtCH(4)(+), identified as the H-Pt(+)-CH(3) intermediate, with Xe reveals a bond energy of 1.77 +/- 0.08 eV (171 +/- 8 kJ/mol) relative to Pt(+) + CH(4). The experimental thermochemistry is favorably compared with density functional theory calculations (B3LYP using several basis sets), which also establish the electronic structures of these species and provide insight into the reaction mechanism. Results for the reaction of Pt(+) with methane are compared with those for the analogous palladium system and the differences in reactivity and mechanism are discussed.     

Arachno, nido, and closo aromatic isomers of the Li6B6H6 molecule

We analyzed chemical bonding in low-lying isomers of the recently computationally predicted B(6)H(6)Li(6) molecule. According to our calculations the benzene-like B(6)H(6)Li(6) (D(2h), (1)A(1g)) arachno structure with the planar aromatic B(6)H(6)(6-) anion is the most stable one. A nido isomer with two aromatic B(6)H(6)(4-) (pentagonal pyramid) and Li(3)(+) (triangular) moieties, which can be considered as derived from the global minimum structure through a two-electron intramolecular transfer from B(6)H(6)(6-) to three Li(+) cations, was found to be 10.7 kcal/mol higher in energy. A closo isomer with three aromatic moieties (octahedral B(6)H(6)(2-) and two Li(3)(+)) was found to be 31.3 kcal/mol higher in energy than the global minimum. Another isomer with three aromatic moieties (two B(3)H(3)(2-) and Li(3)(+)) was found to be substantially higher in energy (74.4 kcal/mol). Thus, the intramolecular electron transfers from the highly charged B(6)H(6)(6-) anion to cations are not favorable for the B(6)H(6)Li(6) molecule, even when a formation of three-dimensional aromatic B(6)H(6)(2-) anion and two sigma-aromatic Li(3)(+) cations occurs in the closo isomer.     

Theoretical proton affinity and fluoride affinity of nerve agent VX

Proton affinity and fluoride affinity of nerve agent VX at all of its possible sites were calculated at the RI-MP2/cc-pVTZ//B3LYP/6-31G* and RI-MP2/aug-cc-pVTZ//B3LYP/6-31+G* levels, respectively. The protonation leads to various unique structures, with H(+) attached to oxygen, nitrogen, and sulfur atoms; among which the nitrogen site possesses the highest proton affinity of -ΔE ～ 251 kcal/mol, suggesting that this is likely to be the major product. In addition some H(2), CH(4) dissociation as well as destruction channels have been found, among which the CH(4) + [Et-O-P(═O)(Me)-S-(CH(2))(2)-N(+)(iPr)═CHMe] product and the destruction product forming Et-O-P(═O)(Me)-SMe + CH(2)═N(+)(iPr)(2) are only 9 kcal/mol less stable than the most stable N-protonated product. For fluoridization, the S-P destruction channel to give Et-O-P(═O)(Me)(F) + [S-(CH(2))(2)-N-(iPr)(2)](-) is energetically the most favorable, with a fluoride affinity of -ΔE ～ 44 kcal. Various F(-) ion-molecule complexes are also found, with the one having F(-) interacting with two hydrogen atoms in different alkyl groups to be only 9 kcal/mol higher than the above destruction product. These results suggest VX behaves quite differently from surrogate systems.