Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with nitric oxide

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with NO. 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 299-592 K. The second-order rate constants at 10 Torr fitted the Arrhenius equation log(k/cm(3) molecule(-1) s(-1)) = (- 11.66 +/- 0.01) + (6.20 +/- 0.10 kJ mol(-1))IRT In 10 The rate constants showed a variation with pressure of a factor of ca. 2 over the available range, almost independent of temperature. The data could not be fitted by RRKM calculations to a simple third body assisted association reaction alone. However, a mechanistic model with an additional (pressure independent) side channel gave a reasonable fit to the data. Ab initio calculations at the G3 level supported a mechanism in which the initial adduct, bent H2SiNO, can ring close to form cyclo-H2SiNO, which is partially collisionally stabilized. In addition, bent H2SiNO can undergo a low barrier isomerization reaction leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are NH2 + SiO. The rate controlling barrier for this latter pathway is only 16 kJ mol(-1) below the energy of SiH2 + NO. This is consistent with the kinetic findings. A particular outcome of this work is that, despite the pressure dependence and the effects of the secondary barrier (in the side reaction), the initial encounter of SiH2 with NO occurs at the collision rate. Thus, silylene can be as reactive with odd electron molecules as with many even electron species. Some comparisons are drawn with the reactions of CH2 + NO and SiCl2 + NO.

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.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 1. HCl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with HCL The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas, at five temperatures in the range of 296-611 K. The second-order rate constants fitted the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.51 +/- 0.06) + (1.92 +/- 0.47 kJ mol(-1))/RTIn10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range of 10-100 Torr, but showed small decreases in value of no more than 20% ( +/- 10%) at I Toff, at both the highest and lowest temperatures. The data are consistent with formation of an initial weakly bound donor-acceptor complex, which reacts by two parallel pathways. The first is by chlorine-to-silicon H-shift to make vibrationally excited chlorosilane, SiH3Cl*, which yields HSiCl by H-2 elimination from silicon. In the second pathway, the complex proceeds via H-2 elimination (4-center process) to make chlorosilylene, HSiCl, directly. This interpretation is supported by ab initio quantum calculations carried out at the G3 level which reveal the direct H-2 elimination route for the first time. RRKM modeling predicts the approximate magnitude of the pressure effect but is unable to determine the proportions of each pathway. The experimental data agree with the only previous measurements at room temperature. Comparisons with other reactions of SiH2 are also drawn.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 2. CH3Cl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of both silacyclopent-3-ene and phenylsilane, have been carried out to obtain second-order rate constants for its reaction with CH3Cl. The reaction was studied in the gas phase at six temperatures in the range 294-606 K. The second-order rate constants gave a curved Arrhenius plot with a minimum value at T approximate to 370 K. The reaction showed no pressure dependence in the presence of up to 100 Torr SF6. The rate constants, however, showed a weak dependence on laser pulse energy. This suggests an interpretation requiring more than one contributing reaction pathway to SiH2 removal. Apart from a direct reaction of SiH2 with CH3Cl, reaction of SiH2 with CH3 (formed by photodissociation of CH3Cl) seems probable, with contributions of up to 30% to the rates. Ab initio calculations (G3 level) show that the initial step of reaction of SiH2 with CH3Cl is formation of a zwitterionic complex (ylid), but a high-energy barrier rules out the subsequent insertion step. On the other hand, the Cl-abstraction reaction leading to CH3 + ClSiH2 has a low barrier, and therefore, this seems the most likely candidate for the main reaction pathway of SiH2 with CH3Cl. RRKM calculations on the abstraction pathway show that this process alone cannot account for the observed temperature dependence of the rate constants. The data are discussed in light of studies of other silylene reactions with haloalkanes.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

Experimental and theoretical evidence for homogeneous catalysis in the gas-phase reaction of SiH2 with H2O (and D2O): a combined kinetic and quantum chemical study

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 297 K and at 345 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied independently as a function of H2O (or D2O) and SF6 (bath gas) pressures. At a fixed pressure of SF6 (5 Torr), [SiH2] decay constants, k(obs), showed a quadratic dependence on [H2O] or [D2O]. At a fixed pressure of H2O or D2O, k(obs) Values were strongly dependent on [SF6]. The combined rate expression is consistent with a mechanism involving the reversible formation of a vibrationally excited zwitterionic donor-acceptor complex, H2Si...OH2 (or H2Si...OD2). This complex can then either be stabilized by SF6 or it reacts with a further molecule of H2O (or D2O) in the rate-determining step. Isotope effects are in the range 1.0-1.5 and are broadly consistent with this mechanism. The mechanism is further supported by RRKM theory, which shows the association reaction to be close to its third-order region of pressure (SF6) dependence. Ab initio quantum calculations, carried out at the G3 level, support the existence of a hydrated zwitterion H2Si...(OH2)(2), which can rearrange to hydrated silanol, with an energy barrier below the reaction energy threshold. This is the first example of a gas-phase-catalyzed silylene reaction.

Quantum mechanics studies of the tautomers of dlhydrodiacetylformoin, an important Maillard intermediate and the synthesis of furaneol

Conformational analyses have been carried out on the acyclic and cyclic forms of dihydrodiacetylformoin, an important Maillard intermediate and precursor for furaneol. For the acyclic forms, the 2,5-dicarbonyl isomers have the lowest energy, while for the cyclic forms, the 3-carbonyl are favoured over the 4-carbonyl isomers. The likely path for cyclisation is investigated and it is shown that the favoured path is dependent upon the relative chiralities of the carbon atoms and in particular that the reaction proceeds more readily if C2 and C3 have different chiralities. After cyclisation, the reaction path to produce furaneol proceeds via the loss of a water molecule. This reaction has been studied with a model including two water molecules and a hydroxide anion and shows relatively low-energy barriers. (C) 2008 Published by Elsevier B.V.

Quantum studies of the vibrations in H3O2- and D3O2-

The vibrations of H3O2- and D3O2- are investigated using diffusion Monte Carlo (DMC) and vibrational configuration-interaction approaches, as implemented in the program MULTIMODE. These studies use the potential surface recently developed by Huang [ J. Am. Chem. Soc. 126, 5042 (2004)]. The focus of this work is on the vibrational ground state and fundamentals which occur between 100 and 3700 cm(-1). In most cases, excellent agreement is obtained between the fundamental frequencies calculated by the two approaches. This serves to demonstrate the power of both methods for treating this very anharmonic system. Based on the results of the MULTIMODE and DMC treatments, the extent and nature of the couplings in H3O2- and D3O2- are investigated. (C) 2005 American Institute of Physics.

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.

Comparison of white metal and resin socket terminations for wire ropes

Previously the authors have presented both theoretical and experimental work discussing the operating mechanism of a wire rope held in a tapered socket by means of a cast resin cone. The work reported here extends the investigation to address the question of whether the same socket fabricated with white metal operates in the same manner. To date, previous investigations have compared the operational efficiency of resin and white metal in terms of both strength and/or fatigue endurance. Some other work has analysed the operation of resin sockets or specific cast metal terminations. This paper seeks to draw the results from this work together, and, in addition to a theoretical analysis, presents experimental data obtained from a direct comparison of the operation mechanism for the same sockets filled with resin or white metal. Results show that white metal terminations have a very different distribution of stresses along the length of the socket basket from resin terminations, and a smaller but still significant amount of socket draw. For both types of termination the socket draw develops high frictional gripping forces which can transfer the load from the rope to the socket. The different stress distributions mean that the consequences of termination fabrication defects may not be the same for resin and white metal terminations.

Modelling the torsional interaction of wire and polyester fibre ropes used for offshore moorings

Oil rig mooring lines have traditionally consisted of chain and wire rope. As production has moved into deeper water it has proved advantageous to incorporate sections of fibre rope into the mooring lines. However, this has highlighted torsional interaction problems that can occur when ropes of different types are joined together. This paper describes a method by which the torsional properties of ropes can be modelled and can then be used to calculate the rotation and torque for two ropes connected in series. The method uses numerical representations of the torsional characteristics of both the ropes, and equates the torque generated in each rope under load to determine the rotation at the connection point. Data from rope torsional characterization tests have been analysed to derive constants used in the numerical model. Constants are presented for: a six-strand wire rope; a torque-balanced fibre rope; and a fibre rope that has been designed to be torque-matched to stranded wire rope. The calculation method has been verified by comparing predicted rotations with measured test values. Worked examples are given for a six-strand wire rope connected, firstly, to a torque-balanced fibre rope that offers little rotational restraint, and, secondly, to a fibre rope whose torsional properties are matched to that of the wire rope.