Computational characterization of a complete palladium-catalyzed cross-coupling process: the associative transmetalation in the Stille reaction

[reaction: see text] The species presumably involved in the associative ligand substitution mechanism for the Stille cross coupling of vinyl bromide and trimethylvinyl stannane with Pd(PMe(3))(2)/PMe(3) as catalysts in DMF (as ligand and solvent) have been structurally and energetically characterized. The cyclic four-coordinate transition state for the rate-determining transmetalation step explains the retention of configuration in the Stille coupling of chiral nonracemic alkyl stannanes.     

Theoretical and experimental analysis of ammonia ionic clusters produced by 252Cf fragment impact on an NH3 ice target

Positive and negatively charged ammonia clusters produced by the impact of (252)Cf fission fragments (FF) on an NH(3) ice target have been examined theoretical and experimentally. The ammonia clusters generated by (252)Cf FF show an exponential dependence of the cluster population on its mass, and the desorption yields for the positive (NH(3))(n)NH(4)(+) clusters are 1 order of magnitude higher than those for the negative (NH(3))(n)NH(2)(-) clusters. The experimental population analysis of (NH(3))(n)NH(4)(+) (n = 0-18) and (NH(3))(n)NH(2)(-) (n = 0-8) cluster series show a special stability at n = 4 and 16 and n = 2, 4, and 6, respectively. DFT/B3LYP calculations of the (NH(3))(0)(-)(8)NH(4)(+) clusters show that the structures of the more stable conformers follow a clear pattern: each additional NH(3) group makes a new hydrogen bond with one of the hydrogen atoms of an NH(3) unit already bound to the NH(4)(+) core. For the (NH(3))(0)(-)(8)NH(2)(-) clusters, the DFT/B3LYP calculations show that, within the calculation error, the more stable conformers follow a clear pattern for n = 1-6: each additional NH(3) group makes a new hydrogen bond to the NH(2)(-) core. For n = 7 and 8, the additional NH(3) groups bind to other NH(3) groups, probably because of the saturation of the NH(2)(-) core. Similar results were obtained at the MP2 level of calculation. A stability analysis was performed using the commonly defined stability function E(n)(-)(1) + E(n)(+1) - 2E(n), where E is the total energy of the cluster, including the zero point correction energy (E = E(t) + ZPE). The trend on the relative stability of the clusters presents an excellent agreement with the distribution of experimental cluster abundances. Moreover, the stability analysis predicts that the (NH(3))(4)NH(4)(+) and the even negative clusters [(NH(3))(n)NH(2)(-), n = 2, 4, and 6] should be the most stable ones, in perfect agreement with the experimental results.     

Rate coefficient for the reaction SiO + Si2O2 at T = 10-1000 K

The reaction paths for the formation of Si3O3 molecules have been investigated at high level ab initio quantum chemical calculations by using the QCISD method with the 6-311++G(d,p) basis set. The cis-Si2O2 isomer does not participate in the chemical mechanism for the formation of Si3O3 molecules. Although the SiO + cis-Si2O2 reaction is exothermic and spontaneous, it is not expected to explain the growth mechanism of Si3O3 in the interstellar silicate grains of circumstellar envelopes surrounding M-type giants. The reaction of SiO with cyclic Si2O2 molecules is exothermic, is spontaneous, and has a nonplanar transition state. The Gibbs free energy for the transition state formation, (DeltaG0#), is around 5.5 kcal mol-1 at 298 K. The bimolecular rate coefficient for this reaction, kT, is about 1 x 10-12 cm3 molecule-1 s-1 at 298 K and in the collision limit, 1.5 x 10-10 cm3 molecule-1 s-1, at 500 K. The activation energy, Ea, is about 8 kcal mol-1. The enthalpy of Si3O3 fragmentation is 53.9 kcal mol-1 at 298 K. The SiO + cyclic Si2O2 reaction is expected to be the most prominent reaction path for the Si3O3 formation in interstellar environment and fabrication of silicon nanowires.     

Parametrizing PCM to obtain solvation free energies from group contributions

A parametrization methodology for evaluating the solvation free energy, using the polarizable continuum model implemented in Gamess software, is presented in a formulation which makes use of a group contribution conception to construct the cavities. The systems studied include alkanes, alcohols, aldehydes and ketones embeded in a continuous medium simulating the water as the solvent. For each family, the CH2, OH, and C=O moieties of atoms are put together in single spheres forming a group. The cavities are constructed in two different ways, one for the electrostatic component and the other for nonelectrostatic contributions, i.e., the cavitation, dispersion, and repulsion components of free energy of solvation. A multivariate analysis is performed to obtain an assembly of variables, for each homologous series, able to give the results which are close to experiment. The analysis is addressed in order to (i) compare the theoretical free energy of solvation with the experimental trends of the solutes in aqueous media, when the chain is increased, (ii) compare the behavior of each component of free energy with the increasing CH2 number, (iii) investigate the influence of the oxygen atom on the components, and (iv) quantify the relative contribution of each component to the final free energy of solvation for some homologous series.     

Energetics of the thermal dimerization of acenaphthylene to heptacyclene

The energetics of the thermal dimerization of acenaphthylene to give Z- or E-heptacyclene was investigated. The standard molar enthalpy of the formation of monoclinic Z- and E-heptacyclene isomers at 298.15 K was determined as Delta(f)H(m)o (E-C24H16, cr) = 269.3 +/- 5.6 kJ x mol(-1) and Delta(f)H(m)o (Z-C24H16, cr) = 317.7 +/- 5.6 kJ x mol(-1), respectively, by microcombustion calorimetry. The corresponding enthalpies of sublimation, Delta(sub)H(m)o (E-C24H16) = (149.0 +/- 3.1) kJ x mol(-1) and Delta(sub)H(m)o (Z-C24H16) = (128.5 +/- 2.3) kJ x mol(-1) were also obtained by Knudsen effusion and Calvet-drop microcalorimetry methods, leading to Delta(f)H(m)o (E-C24H16, g) = (418.3 +/- 6.4) kJ x mol(-1) and Delta(f)H(m)o (Z-C24H16, g) = (446.2 +/- 6.1) kJ x mol(-1), respectively. These results, in conjunction with the reported enthalpies of formation of solid and gaseous acenaphthylene, and the entropies of acenaphthylene and both hepatcyclene isomers obtained by the B3LYP/6-31G(d,p) method led to the conclusion that at 298.15 K the thermal dimerization of acenaphthylene is considerably exothermic and exergonic in the solid and gaseous states (although more favorable when the E isomer is the product), suggesting that the nonobservation of the reaction under these conditions is of kinetic nature. A full determination of the molecular and crystal structure of the E dimer by X-ray diffraction is reported for the first time. Finally, molecular dynamics computer simulations on acenaphthylene and the heptacyclene solids were carried out and the results discussed in light of the corresponding structural and Delta(sub)H(m)o data experimentally obtained.     

Multilevel and density functional electronic structure calculations of proton affinities and gas-phase basicities involved in biological phosphoryl transfer

Five multilevel model chemistries (CBS-QB3, G3B3, G3MP2B3, MCG3/3, and MC-QCISD/3) and seven hybrid density functional methods (PBE0, B1B95, B3LYP, MPW1KCIS, PBE1KCIS, and MPW1B95) have been applied to the calculation of gas-phase basicity and proton affinity values for a series of 17 molecules relevant to the study of biological phosphoryl transfer. In addition, W1 calculations were performed on a subset of molecules. The accuracy of the methods was assessed and the nature of systematic errors was explored, leading to the introduction of a set of effective bond enthalpy and entropy correction terms. The multicoefficient correlation methods (MCG3/3 and MC-QCISD), with inclusion of specific zero-point scale factors, slightly outperform the other multilevel methods tested (CBS-QB3, G3B3, and G3MP2B3), with significantly less computational cost, and in the case of MC-QCISD, slightly less severe scaling. Four density functional methods, PBE1KCIS, MPW1B95, PBE0, and B1B95 perform nearly as well as the multilevel methods. These results provide an important set of benchmarks relevant to biological phosphoryl transfer reactions.     

Fluorous phase-transfer activation of catalysts: application of a new rate-enhancement strategy to alkene metathesis

Reactions of the bis(pyridine) complex (H2IMes)(Py)2(Cl)2Ru(=CHPh) and fluorous phosphines P(CH2CH2R(fn))3 (n = a, 6; b, 8; c, 10; R(fn) = (CF2)(n-1)CF3) give (H2IMes)(P(CH2CH2R(fn))3)(Cl)2Ru(=CHPh) (2a-c, 64-73%), which are analogs of Grubbs' second generation catalyst and effective alkene metathesis catalysts under organic monophasic and fluorous/organic biphasic conditions. The latter give rate accelerations, which are believed to arise from phase transfer of the dissociated fluorous phosphine.     

Sparkle/AM1 structure modeling of lanthanum (III) and lutetium (III) complexes

The sparkle/AM1 model for the quantum chemical prediction of coordination polyhedron crystallographic geometries from isolated lanthanide complex ion calculations, defined recently for Eu(III), Gd(III), and Tb(III) (Inorg. Chem. 2005, 44, 3299) is now extended to La(III) and Lu(III). Thus, for each of the metal ions we chose a training set of 15 complexes that possess various representative ligands of high crystallographic quality (R factor < 0.05 Angstroms) and oxygen and/or nitrogen as coordinating atoms. In the validation procedure we used a set of 60 more La(III) coordination compound structures, as well as 15 more Lu(III) coordination compound structures, all of high crystallographic quality. For both the 75 La(III) compounds and the 30 Lu(III) compounds, the Sparkle/AM1 unsigned mean error, for all interatomic distances between the metal ions and the ligand atoms of the first sphere of coordination, is 0.08 Angstroms, thus comparable to the accuracy normally achievable by present day ab initio/ECP calculations, while being hundreds of times faster.     

Computational study on the bond dissociation enthalpies in the enolic and ketonic forms of beta-diketones: their influence on metal-ligand bond enthalpies

A computational study on the thermodynamic properties of 13 beta-diketones is presented. The B3LYP//6-311+G(2d,2p)//B3LYP/6-31G(d) theoretical approach was employed to compute the O-H and C-H bond dissociation enthalpies and enthalpy of tautomerization and to estimate standard gas-phase enthalpies of formation for the radicals and for the parent molecules. The gas-phase enthalpies of formation for the neutral molecules are in excellent agreement with available experimental data, supporting the estimates made for the radicals. The latter are very important for the clarification of the thermochemistry of many beta-diketonato metal complexes previously reported in the literature. Importantly, when substituents R = -CHR' are attached to the beta-diketone's scaffold, C-H homolytic bond cleavage is always favored with respect to O-H bond scission.