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.

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.

Non-stoichiometric oxide and metal interfaces and reactions

We have employed a combination of experimental surface science techniques and density functional calculations to study the reduction of TiO2(110) surfaces through the doping with submonolayer transition metals. We concentrate on the role of Ti adatoms in self doping of rutile and contrast the behaviour to that of Cr. DFT+U calculations enable identification of probable adsorption structures and their spectroscopic characteristics. Adsorption of both metals leads to a broken symmetry and an asymmetric charge transfer localised around the defect site of a mixed localised/delocalised character. Charge transfer creates defect states with Ti 3d character in the band gap at similar to 1-eV binding energy. Cr adsorption, however, leads to a very large shift in the valence-band edge to higher binding energy and the creation of Cr 3d states at 2.8-eV binding energy. Low-temperature oxidation lifts the Ti-derived band-gap states and modifies the intensity of the Cr features, indicative of a change of oxidation state from Cr3+ to Cr4+. Higher temperature processing leads to a loss of Cr from the surface region, indicative of its substitution into the bulk.

Isomeric olefin tetracarbonyl complexes of tungsten(I): an infrared spectroelectrochemical study at low temperatures

In situ electrolysis within an optically transparent thin-layer electrochemical (OTTLE) cell was applied at 293-243 K in combination with FTIR spectroscopy to monitor spectral changes in the carbonyl stretching region accompanying oxidation of four tetracarbonyl olefin complexes of tungsten(0), viz., trans-[W(CO)(4)(eta(2)-ethene)(2)], trans-[W(CO)(4)(eta(2)-norbornene)(2)], [W(CO)(4)(eta(4)-cycloocta-1,5-diene)], and [W(CO)(4)(eta(4)-norbornadiene)]. In all cases, the one-electron-oxidized radical cations (17-electron complexes) have been identified by their characteristic nu(CO) patterns. For the bidentate diene ligands, the cis stereochemistry is essentially fixed in both the 18- and 17-electron complexes. The radical cation of the trans-bis(ethene) complex was observed only at 243 K, while at room temperature it isomerized rapidly to the corresponding cis-isomer. The thermal stability of the three studied radical cations in the cis configuration correlates with the relative strength of the W-CO bonds in the positions trans to the olefin ligand, which are more affected by the oxidation than the axial W-CO bonds. For the bulky norbornene ligands, their trans configuration in the bis(norbornene) complex remains preserved after the oxidation in the whole temperature range studied. The limited thermal stability of the radical cations of the trans-bis(alkene) complexes is ascribed to dissociation of the alkene ligands. The spectroelectrochemical results are in very good agreement with data obtained earlier by DFT (B3LYP) calculations.

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.