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).

Antimalarial drugs based on artemisinin: DFT calculations on the principal reactions

Theoretical calculations have been carried out on the interactions of several endoperoxides which are potential antimalarials, including the clinically useful artemisinin, with two possible sources of iron in the parasite, namely the hexa-aquo ferrous ion [Fe(H2O)(6)](2+) and haeme. DFT calculations show that the reactions of all endoperoxides considered, with both sources of iron, initially generate a Fe-O bond followed by cleavage of the O-O bond to oxygen radical species. Subsequently, they can be transformed into carbon-centred radicals of greater stability. However, with [Fe(H2O)(6)](2+) as the iron source, the oxygen-centred radical species are more likely to react further akin to Fenton's reagent, whereby iron salts encourage hydrogen peroxide to act as an oxidizing agent, and that solvent plays a major role. In contrast, when reacting with haeme, the oxygen-centred radicals interconvert to more stable carbon-centred radicals, which can then alkylate haeme. Subsequent cleavage of the Fe-O bond leads to stable and inactive antimalarial products. These results indicate that the reactivity of the endoperoxides as antimalarials is greater with iron hexahydrates for radical-mediated damage as opposed to haeme, which leads to unreactive species. Since only nanomolar quantities of hydrated metal ions could catalyse the reactions leading to damage to the parasites, this could be an alternative or competitive reaction responsible for the antimalarial activity. (c) 2005 Elsevier B.V. All rights reserved.

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.

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.

Tuning of redox potential and visible absorption band of ruthenium(II) complexes of (benzimidazolyl) derivatives: synthesis, characterization, spectroscopic and redox properties, X-ray structures and DFT calculations

Six ruthenium(II) complexes have been prepared using the tridentate ligands 2,6-bis(benzimidazolyl) pyridine and bis(2-benzimidazolyl methyl) amine and having 2,2'-bipyridine, 2,2':6',2 ''-terpyridine, PPh3, MeCN and chloride as coligands. The crystal structures of three of the complexes trans-[Ru(bbpH(2))(PPh3)(2)(CH3CN)I(ClO4)(2) center dot 2H(2)O (2), [Ru(bbpH(2))(bpy)Cl]ClO4 (3) and [Ru(bbpH(2))(terpy)](ClO4)(2) (4) are also reported. The complexes show visible region absorption at 402-517 nm, indicating that it is possible to tune the visible region absorption by varying the ancillary ligand. Luminescence behavior of the complexes has been studied both at RT and at liquid nitrogen temperature (LNT). Luminescence of the complexes is found to be insensitive to the presence of dioxygen. Two of the complexes [Ru(bbpH(2))(bpy)Cl]ClO4 (3) and [Ru(bbpH(2))(terpy]ClO4)(2) (4) show RT emission in the NIR region, having lifetime, quantum yield and radiative constant values suitable for their application as NIR emitter in the solid state devices. The DFT calculations on these two complexes indicate that the metal t(2g) electrons are appreciably delocalized over the ligand backbone. (C) 2006 Elsevier B.V. All rights reserved.

Trinuclear cu(ll) complexes containing peripheral ketonic oxygen bridges and a mu(3)-OH core: Syntheses, crystal structures, spectroscopic and magnetic properties

Four new trinuclear copper(II) complexes, [(CuL1)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (1), [(CuL2)(3)(mu(3)-OH)](CIO4)(2) (2), [(CuL3)(3)-(mu(3)-OH)](ClO4)(4)center dot H2O (3), and [(CuL4)(3)(mu(3)-OH)](ClO4)(2)center dot H2O (4), where HL1 = 8-amino-4,7,7-trimethyl-5-azaoct-3-en-2-one, HL2 = 7-amino-4-methyl-5-azaoct-3-en-2-one, HL3 = 7(ethylamino)-4-methyl-5-azahept-3-en-2-one, and HL4 = 4-methyl-7-(methylamino)-5-azahept-3-en-2-one, have been derived from the four tridentate Schiff bases (HL1, HL2, HL3, and HL4) and structurally characterized by X-ray crystallography. For all compounds, the cationic part is trinuclear with a CU3OH core held by three carbonyl oxygen bridges between each pair of copper(II) atoms. The copper atoms are five-coordinate with a distorted square-pyramidal geometry; the equatorial plane consists of the bridging oxygen atom of the central OH group together with three atoms (N, N, O) from one ligand whereas an oxygen atom of a second ligand occupies the axial position. Magnetic measurements have been performed in the 2-300 K temperature range. The experimental data could be satisfactorily reproduced by using an isotropic exchange model, H = -J(S1S2+S2S3+S1S3) yielding as best-fit parameters: J = -66.7 and g = 2.19 for 1, J = -36.6 and g = 2.20 for 2, J = -24.5 and g = 2.20 for 3, and J = -14.9 and g = 2.05 for 4. EPR spectra at low temperature show the existence of spin frustration in complexes 3 and 4, but it has not been possible to carry out calculations of the antisymmetric exchange parameter, G, from magnetic data. In frozen methanolic solution, at 4 K, hyperfine splitting in all complexes and spin frustration in complex 4 seem to be confirmed. ((c) Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, Germany, 2005)

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.

Tests of MULTIMODE calculations of rovibrational energies of CH4

We report variational calculations of rovibrational energies of CH4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found. (c) 2006 Elsevier B.V. All rights reserved.

Verification of design calculations of a wind catcher/tower natural ventilation system with performance testing in real building

A wind catcher/tower natural ventilation system was installed in a seminar room in the building of the School of Construction Management and Engineering, the University of Reading in the UK . Performance was analysed by means of ventilation tracer gas measurements, indoor climate measurements (temperature, humidity, CO2) and occupant surveys. In addition, the potential of simple design tools was evaluated by comparing observed ventilation results with those predicted by an explicit ventilation model and the AIDA implicit ventilation model. To support this analysis, external climate parameters (wind speed and direction, solar radiation, external temperature and humidity) were also monitored. The results showed the chosen ventilation design provided a substantially greater ventilation rate than an equivalent area of openable window. Also air quality parameters stayed within accepted norms while occupants expressed general satisfaction with the system and with comfort conditions. Night cooling was maximised by using the system in combination with openable windows. Comparisons of calculations with ventilation rate measurements showed that while AIDA gave reasonably correlated results with the monitored performance results, the widely used industry explicit model was found to over estimate the monitored ventilation rate.

Fast and accurate SCFT calculations for periodic block-copolymer morphologies using the spectral method with Anderson mixing

We study the numerical efficiency of solving the self-consistent field theory (SCFT) for periodic block-copolymer morphologies by combining the spectral method with Anderson mixing. Using AB diblock-copolymer melts as an example, we demonstrate that this approach can be orders of magnitude faster than competing methods, permitting precise calculations with relatively little computational cost. Moreover, our results raise significant doubts that the gyroid (G) phase extends to infinite $\chi N$. With the increased precision, we are also able to resolve subtle free-energy differences, allowing us to investigate the layer stacking in the perforated-lamellar (PL) phase and the lattice arrangement of the close-packed spherical (S$_{cp}$) phase. Furthermore, our study sheds light on the existence of the newly discovered Fddd (O$^{70}$) morphology, showing that conformational asymmetry has a significant effect on its stability.