The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.

Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE

We report rigorous calculations of rovibrational energies and dipole transition intensities for three molecules using a new version of the code MULTIMODE. The key features of this code which permit, for the first time, such calculations for moderately sized but otherwise general polyatomic molecules are briefly described. Calculations for the triatomic molecule BF(2) are done to validate the code. New calculations for H(2)CO and H(2)CS are reported; these make use of semiempirical potentials but ab initio dipole moment surfaces. The new dipole surface for H(2)CO is a full-dimensional fit to the dipole moment obtained with the coupled-cluster with single and double excitations and a perturbative treatment of triple excitations method with the augmented correlation consistent triple zeta basis set. Detailed comparisons are made with experimental results from a fit to relative data for H(2)CS and absolute intensities from the HITRAN database for H(2)CO.

Methylene spacer-regulated structural variation in cobalt(II/III) complexes with bridging acetate and Salen- or Salpn-type Schiff-base ligands

Two linear, trinuclear mixed-valence complexes, [Co-II{(mu-L-1)(mu-OAc)Co-III (OAc)}(2)] (1) and [Co-II(mu-L-2) (mu-OAc)Co-III(OAc)}(2)] (2) and two mononuclear Con' complexes [Co-III{L-3)(OAc)] (3), and [Co-III {L-4}(OAc)] (4) were prepared and the molecular structures of 1, 2 and 4 elucidated on the basis of X-ray crystallography [OAc = Acetate ion, H2L1 = H(2)Salen 1,6-bis(2-hydroxyphenyl)-2,5-diazahexa-1,5-diene, H2L2 H2Me2-Salen = 2,7-bis(2-hydroxyphenyl)-2,6-diazaocta-2,6-diene, H2L3 = H(2)Salpn = 1,7-bis(2-hydroxyphenyl)-2,6-diazahepta1,6-diene, H2L4 = H(2)Me(2)Salpn = 2,8-bis(2-hydroxyphenyl)3,7-diazanona-2,7-dienel. In complexes I and 2, the acetate groups show both monodentate and bridging bidentate coordination modes, whereas chelating bidentate acetate is present in 4. The terminal (CoN2O4)-N-III centres in 1 and 2 exhibit uniform facial arrangements of both non-bridged N2O and bridging O-3 donor sets and the Co-II centre is coordinated to six (four phenoxo and two acetato) oxygen atoms of the bridging ligands. The effective magnetic moment at room temperature corresponds to the presence of high-spin Coll in both 1 and 2. The complexes 1 and 2 are thus Co-III(S = 0)Co-II(S = 3/2)-Co-II(S = 0) trimers. Complexes 3 and 4 are monomeric and diamagnetic containing low-spin Co-III(S = 0) with chelating tetradentate Schiff base and bidentate acetate. Calculations based on DFT rationalise the formation of trinuclear or monomiclear complexes. (C) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2008).

Large vibrational variational calculations using 'multimode' and an iterative diagonalization technique

We have favoured the variational (secular equation) method for the determination of the (ro-) vibrational energy levels of polyatomic molecules. We use predominantly the Watson Hamiltonian in normal coordinates and an associated given potential in the variational code 'Multimode'. The dominant cost is the construction and diagonalization of matrices of ever-increasing size. Here we address this problem, using pertubation theory to select dominant expansion terms within the Davidson-Liu iterative diagonalization method. Our chosen example is the twelve-mode molecule methanol, for which we have an ab initio representation of the potential which includes the internal rotational motion of the OH group relative to CH3. Our new algorithm allows us to obtain converged energy levels for matrices of dimensions in excess of 100 000.

Surface geometry of Cu{531}

We present a combined quantitative low-energy electron diffraction (LEED) and density-functional theory (DFT) study of the chiral Cu{531} surface. The surface shows large inward relaxations with respect to the bulk interlayer distance of the first two layers and a large expansion of the distance between the fourth and fifth layers. (The latter is the first layer having the same coordination as the Cu atoms in the bulk.) Additional calculations have been performed to study the likelihood of faceting by comparing surface energies of possible facet terminations. No overall significant reduction in energy with respect to planar {531} could be found for any of the tested combinations of facets, which is in agreement with the experimental findings.

Cu(I) and Ag(I) complexes of chalcogenide derivatives of the organometallic ligand dppf and the dppa analogue

New Cu(I) and Ag(I) complexes were prepared by reaction of [M(NCCH3)(4)][X] (M = Cu or Ag; X = BF4 or PF6) with the bidentate chalcogenide ligands Ph2P(E)NHP(E)Ph-2 (E = S, S(2)dppa; E = Se, Se(2)dppa), and dpspf (1, 1'-bis(diphenylselenophosphoryl)ferrocene). Copper and silver behaved differently. While three molecules of either S(2)dppa and Se(2)dppa bind to a distorted tetrahedral Cu-4 cluster, with deprotonation of the ligand, 1:2 complexes of the neutral ligands are formed with Ag(l), with a tetrahedral coordination of the metal. The [Cu-4{Ph2P(Se)NP(Se)Ph-2}(3)](+) clusters assemble as dimers, held together by weak Se...Se distances interactions. Another dimer was observed for the [Ag(dpspf)](+) cation, with two short Ag...Se distances. DFT and MP2 calculations indicated the presence of attracting interactions, reflected in positive Mayer indices (MI). The electrochemistry study of this species showed that both oxidation and reduction took place at silver. (C) 2004 Elsevier B.V. All rights reserved.

A DFT study of propylene polymerization using neutral salicyladiminato nickel(II) and palladium(II) as catalysts

Propylene polymerization using salicyladiminato metal catalalysts has been studied using density functional theory at the B3LYP/LANL2DZ level. In particular, the effects on the reaction mechanisms of changing the metal from Pd(II) to Ni(II) have been investigated. While the reaction mechanisms involving the salicyladiminato Ni(II) catalyst have been found to be similar to those established previously for the salicyladiminato Pd(II) catalyst, the nickel catalyst was found to differentiate the trans-O intermediate from the trans-.N intermediate with an energy difference of 46.63 U mol(-1) significantly more than the palladium catalyst for which the energy difference was calculated as 35.82 kJ mol(-1). The energy difference between the trans-O configuration and the trans-N configuration is decreased significantly when combining a molecule of propylene with the catalyst. For the Ni catalyst, the trans-O isomer is more stable than the trans-N isomer to a greater extent than for Pd, so that the insertion of propylene from 20 is relatively less favoured for Ni than for Pd. It is predicted that the mechanism of isomerization from 20 to 2N through a rotational transition state TS2O2N is more appropriate for the Ni catalyst system. The palladium system shows a larger preference for pi-coordination than its nickel counterpart, although the latter possesses a lower reaction barrier. It was found that the occupation of the trans-O position in the asymmetric salicyladiminato catalyst is also more favored by the alkene as it is by the alkyl so that insertion of the alkene may always start from a particular configuration so that specific products are obtained. (c) 2005 Elsevier B.V. All rights reserved.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Chalcogenide-halides of niobium (V). 1. Gas-phase structures of NbOBr3, NbSBr3, and NbSCl3. 2. Matrix infrared spectra and vibrational force fields of NbOBr3, NbSBr3, NbSCl3, and NbOCl3

The molecular structures of NbOBr3, NbSCl3, and NbSBr3 have been determined by gas-phase electron diffraction (GED) at nozzle-tip temperatures of 250 degreesC, taking into account the possible presence of NbOCl3 as a contaminant in the NbSCl3 sample and NbOBr3 in the NbSBr3 sample. The experimental data are consistent with trigonal-pyramidal molecules having C-3v symmetry. Infrared spectra of molecules trapped in argon or nitrogen matrices were recorded and exhibit the characteristic fundamental stretching modes for C-3v species. Well resolved isotopic fine structure (Cl-35 and Cl-37) was observed for NbSCl3, and for NbOCl3 which occurred as an impurity in the NbSCl3 spectra. Quantum mechanical calculations of the structures and vibrational frequencies of the four YNbX3 molecules (Y = O, S; X = Cl, Br) were carried out at several levels of theory, most importantly B3LYP DFT with either the Stuttgart RSC ECP or Hay-Wadt (n + 1) ECP VDZ basis set for Nb and the 6-311 G* basis set for the nonmetal atoms. Theoretical values for the bond lengths are 0.01-0.04 Angstrom longer than the experimental ones of type r(a), in accord with general experience, but the bond angles with theoretical minus experimental differences of only 1.0-1.5degrees are notably accurate. Symmetrized force fields were also calculated. The experimental bond lengths (r(g)/Angstrom) and angles (angle(alpha)/deg) with estimated 2sigma uncertainties from GED are as follows. NbOBr3: r(Nb=O) = 1.694(7), r(Nb-Br) = 2.429(2), angle(O=Nb-Br) = 107.3(5), angle(Br-Nb-Br) = 111.5(5). NbSBr3: r(Nb=S) = 2.134(10), r(Nb-Br) = 2.408(4), angle(S=Nb-Br) = 106.6(7), angle(Br-Nb-Br) = 112.2(6). NbSCl3: Nb=S) = 2.120(10), r(Nb-Cl) = 2.271(6), angle(S=Nb-Cl) = 107.8(12), angle(Cl-Nb-Cl) = 111.1(11).

A new tetrameric Cu-II cluster with square topology exhibiting ferro- and antiferromagnetic magnetic pathways: which is which?

[Cu4L2(bpy)(4)(H2O)(3)](ClO4)(4).2.5H(2)O, 1, a new tetranuclear Cu-II cluster showing square planar geometry, formed with aspartate bridging ligand (L) has been synthesized. The global magnetic coupling is ferromagnetic but theoretical DFT/B3LYP calculations are necessary to assign which Cu-L-Cu side is ferro or antiferromagnetically coupled.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.