ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Gas-phase kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 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.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

The calculation of force constants and normal coordinates IV: XH4 and XH3 molecules

Normal coordinate calculations of XH4 and XH3 molecules are reviewed and discussed. It is shown that for most of these molecules the true values of the force constants in the most General Harmonic Force Field can be uniquely determined only by making use of vibration-rotation interaction constants. It is emphasized that without these extra data the GFF is not determined. The results are compared with various model force fields for these molecules.

Coriolis perturbations in the infrared spectrum of ethylene

Rotational structure has been resolved and analyzed in two of the infrared‐active perpendicular bands of C2H4 vapor: the Type b fundamental band, ν10, at 826 cm—1, and the Type c fundamental band, ν7, at 949 cm—1. Many of the individual PP and RR branch lines have been observed. The analysis has been confined to values of the quantum number K≥3, for which energy levels ethylene shows no detectable deviations from a symmetric‐top rotational structure. The analysis reveals a Coriolis interaction between ν7 and ν10, and between ν4 and ν10, and values of the Coriolis constants ζ7,10z and ζ4,10y are obtained; these are related to normal coordinate calculations for the appropriate symmetry species, and force constants are derived to fit the observed zeta constants. The band center of ν10 has been revised from the original figure of 810 cm—1 to the new value, 826 cm—1, and the inactive frequency ν4 is estimated to lie at 1023±3 cm—1, in good agreement with the previous estimate of 1027 cm—1. The change in the value of ν10 leads to a suggested change in the value of the Raman‐active fundamental ν6 from 1236 to 1222 cm—1. New combination bands have been observed at 2174 cm—1, assigned as ν3+ν10; and at 2252 cm—1, assigned as ν4+ν6; also rotational structure has been resolved and analyzed in the ν6+ν10 band at 2048 cm—1. The new data obtained for the C2H4 molecule are summarized in Table XII, with all of the other data presently available on the vibrational and rotational constants.

Mode-specific chemisorption of CH4 on Pt{110}-(1 x 2) explored by first-principles molecular dynamics

The chemisorption of CH4 on Pt{110}-(1 x 2) has been studied by vibrational analysis of the reaction pathway defined by the potential energy surface and, in time reversal, by first-principles molecular dynamics simulations of CH4 associative desorption, with the electronic structure treated explicitly using density functional theory. We find that the symmetric stretch vibration ν1 is strongly coupled to the reaction coordinate; our results therefore provide a firm theoretical basis for recently reported state-resolved reactivity measurements, which show that excitation of the ν1 normal mode is the most efficient way to enhance the reaction probability

Mode-specificity and transition state-specific energy redistribution in the chemisorption of CH4 on Ni{100}

We have investigated methane (CH4) dissociative chemisorption on the Ni{100} surface by first-principles molecular dynamics (MD) simulations. Our results show that this reaction is mode-specific, with the n1 state being the most strongly coupled to efficient energy flow into the reaction coordinate when the molecule reaches the transition state. By performing MD simulations for two different transition state (TS) structures we provide evidence of TS structure-specific energy redistribution in methane chemisorption. Our results are compared with recently reported state-resolved measurement of methane adsorption probability on nickel surfaces, and we find that a strong correlation exists between the highest vibrational efficacy measured on Ni{100} for the n1 state and the calculated highest fractional vibrational energy content in this mode.

Bond-selective energy redistribution in the chemisorption of CH3D and CD3H on Pt{1 1 0}-(1 2): a first-principles molecular dynamics study

We have investigated the chemisorption of CH3D and CD3H on Pt{11 0}-(1 2) by performing first-principles molecular dynamics simulations of the recombinative desorption of CH3D (from adsorbed methyl and deuterium) and of CD3H (from adsorbed trideuteromethyl and hydrogen). Vibrational analysis of the symmetry adapted internal coordinates of the desorbing molecules shows that excitation of the single C– D (C–H) bond in the parent molecule is strongly correlated with energy excess in the reaction coordinate. The results of the molecular dynamics simulations are consistent with observed mode- and bond-specific reactivity measurements for chemisorption of methane and its isotopomers on platinum and nickel surfaces.

Surprisingly slow reaction of dimethylsilylene with dimethylgermane: time-resolved kinetic studies and related quantum chemical calculations

Time-resolved studies of silylene, SiH2, and dimethylsilylene, SiMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to obtain rate constants for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas phase. SiMe2 + Me2GeH2 was studied at five temperatures in the range 299-555 K. Problems of substrate UV absorption at 193 nm at temperatures above 400 K meant that only three temperatures could be used reliably for rate constant measurement. These rate constants gave the Arrhenius parameters log(A/cm(3) molecule(-1) s(-1)) = -13.25 +/- 0.16 and E-a = -(5.01 +/- 1.01) kJ mol(-1). Only room temperature studies of SiH2 were carried out. These gave values of (4.05 +/- 0.06) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + Me2GeH2 at 295 K) and also (4.41 +/- 0.07) x 10(-10) cm(3) molecule(-1) s(-1) (SiH2 + MeGeH3 at 296 K). Rate constant comparisons show the surprising result that SiMe2 reacts 12.5 times slower with Me2GeH2 than with Me2SiH2. Quantum chemical calculations (G2(MP2,SVP)//B3LYP level) of the model Si-H and Ge-H insertion processes of SiMe2 with SiH4/MeSiH3 and GeH4/MeGeH3 support these findings and show that the lower reactivity of SiMe2 with Ge-H bonds is caused by a higher secondary barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper. Other, related, comparisons of silylene reactivity are also presented.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

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.

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.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Unusual isotope effect in the reaction of chlorosilylene with trimethylsilane-1-d: Absolute rate studies and quantum chemical and Rice–Ramsperger–Kassel–Marcus calculations provide strong evidence for the involvement of an intermediate complex

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1- silacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with trimethylsilane-1-d, Me3SiD, in the gas phase. The reaction was studied at total pressures up to 100 Torr (with and without added SF6) over the temperature range of 295−407 K. The rate constants were found to be pressure independent and gave the following Arrhenius equation: log[(k/(cm3 molecule−1 s−1)] = (−13.22 ± 0.15) + [(13.20 ± 1.00) kJ mol−1]/(RT ln 10). When compared with previously published kinetic data for the reaction of ClSiH with Me3SiH, kinetic isotope effects, kD/kH, in the range from 7.4 (297 K) to 6.4 (407 K) were obtained. These far exceed values of 0.4−0.5 estimated for a single-step insertion process. Quantum chemical calculations (G3MP2B3 level) confirm not only the involvement of an intermediate complex, but also the existence of a low-energy internal isomerization pathway which can scramble the D and H atom labels. By means of Rice−Ramsperger−Kassel−Marcus modeling and a necessary (but small) refinement of the energy surface, we have shown that this mechanism can reproduce closely the experimental isotope effects. These findings provide the first experimental evidence for the isomerization pathway and thereby offer the most concrete evidence to date for the existence of intermediate complexes in the insertion reactions of silylenes.

Homoleptic copper(II) and copper(I) complexes of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline. Six-coordinate copper(I) in solution

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Cu(ClO(4))(2)center dot 6H(2)O in methanol in 3:1 M ratio at room temperature yields light green [CuL(3)](ClO(4))(2)center dot H(2)O (1). The X-ray crystal structure of the hemi acetonitrile solvate [CuL(3)](ClO(4))(2)center dot 0.5CH(3)CN has been determined which shows Jahn-Teller distortion in the CuN(6) core present in the cation [CuL(3)](2+). Complex 1 gives an axial EPR spectrum in acetonitrile-toluene glass with g(parallel to) = 2.262 (A(parallel to) = 169 x 10 (4) cm (1)) and g(perpendicular to) = 2.069. The Cu(II/I) potential in 1 in CH(2)Cl(2) at a glassy carbon electrode is 0.32 V versus NHE. This potential does not change with the addition of extra L in the medium implicating generation of a six-coordinate copper(I) species [CuL(3)](+) in solution. B3LYP/LanL2DZ calculations show that the six Cu-N bond distances in [CuL(3)](+) are 2.33, 2.25, 2.32, 2.25, 2.28 and 2.25 angstrom while the ideal Cu(I)-N bond length in a symmetric Cu(I)N(6) moiety is estimated as 2.25 angstrom. Reaction of L with Cu(CH(3)CN)(4)ClO(4) in dehydrated methanol at room temperature even in 4:1 M proportion yields [CuL(2)]ClO(4) (2). Its (1)H NMR spectrum indicates that the metal in [CuL(2)](+) is tetrahedral. The Cu(II/I) potential in 2 is found to be 0.68 V versus NHE in CH(2)Cl(2) at a glassy carbon electrode. In presence of excess L, 2 yields the cyclic voltammogram of 1. From (1)H NMR titration, the free energy of binding of L to [CuL(2)](+) to produce [CuL(3)](+) in CD(2)Cl(2) at 298 K is estimated as -11.7 (+/-0.2) kJ mol (1).