Syntheses, x-ray structures, photochemistry, redox properties, and DFT calculations of interconvertible fac- and mer-[Mn(SPS)(CO)3] isomers containing a flexible SPS-based pincer ligand

The lithium salt of the anionic SPS pincer ligand composed of a central hypervalent lambda(4)-phosphinine ring bearing two ortho-positioned diphenylphosphine sulfide side arms reacts with [Mn(CO)(5)Br] to give fac-[Mn(SPS)(CO)(3)], This isomer can be converted photochemicaily to mer-[Mn(SPS)(CO)(3)], with a very high quantum yield (0.80 +/- 0.05). The thermal backreaction is slow (taking ca. 8 h at room temperature), in contrast to rapid electrodecatalyzed mer-to-fac isomerization triggered by electrochemical reduction of mer-[Mn(SPS)(CO)(3)]. Both geometric isomers of [Mn(SPS)(CO)(3)] have been characterized by X-ray crystallography. Both isomers show luminescence from a low-lying (IL)-I-3 (SPS-based) excited state. The light emission of fac-[Mn(SPS)(CO)(3)] is largely quenched by the efficient photoisomerization occurring probably from a low-lying Mn-CO dissociative excited state. Density functional theory (DFT) and time-dependent DFT calculations describe the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of fac- and mer-[Mn(CO)(3)(SPS)] as ligand-centered orbitals, largely localized on the phosphinine ring of the SPS pincer ligand. In line with the ligand nature of its frontier orbitals, fac-[Mn(SPS)(CO)(3)] is electrochemically reversibly oxidized and reduced to the corresponding radical cation and anion, respectively. The spectroscopic (electron paramagnetic resonance, IR, and UV-vis) characterization of the radical species provides other evidence for the localization of the redox steps on the SIPS ligand. The smaller HOMO-LUMO energy difference in the case of mer-[Mn(CO)(3)(SPS)], reflected in the electronic absorption and emission spectra, corresponds with its lower oxidation potential compared to that of the fac isomer. The thermodynamic instability of mer-[Mn(CO)(3)(SPS)], confirmed by the DFT calculations, increases upon one-electron reduction and oxidation of the complex.

Excited states of nitro-polypyridine metal complexes and their ultrafast decay. Time-resolved IR absorption, spectroelectrochemistry, and TD-DFT calculations of fac-[Re(Cl)(CO)3(5-Nitro-1,10-phenanthroline)]

The lowest absorption band of fac-[Re(Cl)(CO)(3)(5-NO2-phen)] encompasses two close-lying MLCT transitions. The lower one is directed to LUMO, which is heavily localized on the NO2 group. The UV-vis absorption spectrum is well accounted for by TD-DFT (G03/PBEPBE1/CPCM), provided that the solvent, MeCN, is included in the calculations. Near-UV excitation of fac-[Re(Cl)(CO)(3)(5-NO2-phen)] populates a triplet metal to ligand charge-transfer excited state, (MLCT)-M-3, that was characterized by picosecond time-resolved IR spectroscopy. Large positive shifts of the v(CO) bands upon excitation (+70 cm(-1) for the A'(1) band) signify a very large charge separation between the Re(Cl)(CO)3 unit and the 5-NO2-phen ligand. Details of the excited-state character are revealed by TD-DFT calculated changes of electron density distribution. Experimental excited-state v(CO) wavenumbers agree well with those calculated by DFT. The (MLCT)-M-3 state decays with a ca. 10 ps lifetime (in MeCN) into another transient species, that was identified by TRIR and TD-DFT calculations as an intraligand (3)n pi* excited state, whereby the electron density is excited from the NO2 oxygen lone pairs to the pi* system of 5-NO2-phen. This state is short-lived, decaying to the ground state with a similar to 30 ps lifetime. The presence of an n pi* state seems to be the main factor responsible for the lack of emission and the very short lifetimes of 3 MLCT states seen in all d(6)-metal complexes of nitro-polypyridyl ligands. Localization of the excited electron density in the lowest (MLCT)-M-3 states parallels localization of the extra electron in the reduced state that is characterized by a very small negative shift of the v(CO) IR bands (-6 cm(-1) for A'(1)) but a large downward shift of the v(s)(NO2) IR band. The Re-Cl bond is unusually stable toward reduction, whereas the Cl ligand is readily substituted upon oxidation.

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 north Atlantic marine oundary layer experiment (NAMBLEX): Overview of the campaign held at Mace Head, Ireland, in summer 2002

The North Atlantic Marine Boundary Layer Experiment (NAMBLEX), involving over 50 scientists from 12 institutions, took place at Mace Head, Ireland (53.32° N, 9.90° W), between 23 July and 4 September 2002. A wide range of state-of-the-art instrumentation enabled detailed measurements of the boundary layer structure and atmospheric composition in the gas and aerosol phase to be made, providing one of the most comprehensive in situ studies of the marine boundary layer to date. This overview paper describes the aims of the NAMBLEX project in the context of previous field campaigns in the Marine Boundary Layer (MBL), the overall layout of the site, a summary of the instrumentation deployed, the temporal coverage of the measurement data, and the numerical models used to interpret the field data. Measurements of some trace species were made for the first time during the campaign, which was characterised by predominantly clean air of marine origin, but more polluted air with higher levels of NOx originating from continental regions was also experienced. This paper provides a summary of the meteorological measurements and Planetary Boundary Layer (PBL) structure measurements, presents time series of some of the longer-lived trace species (O3, CO, H2, DMS, CH4, NMHC, NOx, NOy, PAN) and summarises measurements of other species that are described in more detail in other papers within this special issue, namely oxygenated VOCs, HCHO, peroxides, organo-halogenated species, a range of shorter lived halogen species (I2, OIO, IO, BrO), NO3 radicals, photolysis frequencies, the free radicals OH, HO2 and (HO2+Σ RO2), as well as a summary of the aerosol measurements. NAMBLEX was supported by measurements made in the vicinity of Mace Head using the NERC Dornier-228 aircraft. Using ECMWF wind-fields, calculations were made of the air-mass trajectories arriving at Mace Head during NAMBLEX, and were analysed together with both meteorological and trace-gas measurements. In this paper a chemical climatology for the duration of the campaign is presented to interpret the distribution of air-mass origins and emission sources, and to provide a convenient framework of air-mass classification that is used by other papers in this issue for the interpretation of observed variability in levels of trace gases and aerosols.

Experimental measurement and numerical simulation of horizontal-coupled slinky ground source heat exchangers

Results from both experimental measurements and 3D numerical simulations of Ground Source Heat Pump systems (GSHP) at a UK climate are presented. Experimental measurements of a horizontal-coupled slinky GSHP were undertaken in Talbot Cottage at Drayton St Leonard site, Oxfordshire, UK. The measured thermophysical properties of in situ soil were used in the CFD model. The thermal performance of slinky heat exchangers for the horizontal-coupled GSHP system for different coil diameters and slinky interval distances was investigated using a validated 3D model. Results from a two month period of monitoring the performance of the GSHP system showed that the COP decreased with the running time. The average COP of the horizontal-coupled GSHP was 2.5. The numerical prediction showed that there was no significant difference in the specific heat extraction of the slinky heat exchanger at different coil diameters. However, the larger the diameter of coil, the higher the heat extraction per meter length of soil. The specific heat extraction also increased, but the heat extraction per meter length of soil decreased with the increase of coil central interval distance.

Gas-phase kinetics of chlorosilylene reactions. I. ClSiH + Me3SiH: absolute rate measurements and theoretical calculations for prototype Si−H insertion reactions

Time-resolved studies of chlorosilylene, CISiH, 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, Me3SiH, in the gas phase. The reaction was studied at total pressures up to 100 torr (with and without added SF6) over the temperature range 297-407 K. The rate constants were found to be pressure independent and gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-13.97 +/- 0.25) + (12.57 +/- 1.64) kJ mol(-1)/RT In 10. The Arrhenius parameters are consistent with a mechanism involving an intermediate complex, whose rearrangement is the rate-determining step. Quantum chemical calculations of the potential energy surface for this reaction and also the reactions of CISiH with SiH4 and the other methylsilanes support this conclusion. Comparisons of both experiment and theory with the analogous Si-H insertion processes of SiH2 and SiMe2 show that the main factor causing the lower reactivity of ClSiH is the secondary energy barrier. The calculations also show the existence of a novel intramolecular H-atom exchange process in the complex of ClSiH with MeSiH3.

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.

Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).

A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

Numerical modeling of inhaled charged aerosol deposition in human airways

A new numerical modeling of inhaled charge aerosol has been developed based on a modified Weibel's model. Both the velocity profiles (slug and parabolic flows) and the particle distributions (uniform and parabolic distributions) have been considered. Inhaled particles are modeled as a dilute dispersed phase flow in which the particle motion is controlled by fluid force and external forces acting on particles. This numerical study extends the previous numerical studies by considering both space- and image-charge forces. Because of the complex computation of interacting forces due to space-charge effect, the particle-mesh (PM) method is selected to calculate these forces. In the PM technique, the charges of all particles are assigned to the space-charge field mesh, for calculating charge density. The Poisson's equation of the electrostatic potential is then solved, and the electrostatic force acting on individual particle is interpolated. It is assumed that there is no effect of humidity on charged particles. The results show that many significant factors also affect the deposition, such as the volume of particle cloud, the velocity profile and the particle distribution. This study allows a better understanding of electrostatic mechanism of aerosol transport and deposition in human airways.

The distribution of solar wind speeds during solar minimum: calibration for numerical solar wind modeling constraints on the source of the slow solar wind

It took the solar polar passage of Ulysses in the early 1990s to establish the global structure of the solar wind speed during solar minimum. However, it remains unclear if the solar wind is composed of two distinct populations of solar wind from different sources (e.g., closed loops which open up to produce the slow solar wind) or if the fast and slow solar wind rely on the superradial expansion of the magnetic field to account for the observed solar wind speed variation. We investigate the solar wind in the inner corona using the Wang-Sheeley-Arge (WSA) coronal model incorporating a new empirical magnetic topology–velocity relationship calibrated for use at 0.1 AU. In this study the empirical solar wind speed relationship was determined by using Helios perihelion observations, along with results from Riley et al. (2003) and Schwadron et al. (2005) as constraints. The new relationship was tested by using it to drive the ENLIL 3-D MHD solar wind model and obtain solar wind parameters at Earth (1.0 AU) and Ulysses (1.4 AU). The improvements in speed, its variability, and the occurrence of high-speed enhancements provide confidence that the new velocity relationship better determines the solar wind speed in the outer corona (0.1 AU). An analysis of this improved velocity field within the WSA model suggests the existence of two distinct mechanisms of the solar wind generation, one for fast and one for slow solar wind, implying that a combination of present theories may be necessary to explain solar wind observations.