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.

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.

Export of dissolved organic carbon from an upland peatland during storm events: implications for flux estimates

Most of the dissolved organic carbon (DOC) exported from catchments is transported during storm events. Accurate assessments of DOC fluxes are essential to understand long-term trends in the transport of DOC from terrestrial to aquatic systems, and also the loss of carbon from peatlands to determine changes in the source/sink status of peatland carbon stores. However, many long-term monitoring programmes collect water samples at a frequency (e.g. weekly/monthly) less than the time period of a typical storm event (typically <1–2 days). As widespread observations in catchments dominated by organo-mineral soils have shown that both concentration and flux of DOC increases during storm events, lower frequency monitoring could result in substantial underestimation of DOC flux as the most dynamic periods of transport are missed. However, our intensive monitoring study in a UK upland peatland catchment showed a contrasting response to these previous studies. Our results showed that (i) DOC concentrations decreased during autumn storm events and showed a poor relationship with flow during other seasons; and that (ii) this decrease in concentrations during autumn storms caused DOC flux estimates based on weekly monitoring data to be over-estimated, rather than under-estimated, because of over rather than under estimation of the flow-weighted mean concentration used in flux calculations. However, as DOC flux is ultimately controlled by discharge volume, and therefore rainfall, and the magnitude of change in discharge was greater than the magnitude of decline in concentrations, DOC flux increased during individual storm events. The implications for long-term DOC trends are therefore contradictory, as increased rainfall could increase flux but cause an overall decrease in DOC concentrations from peatland streams. Care needs to be taken when interpreting long-term trends in DOC flux rather than concentration; as flux is calculated from discharge estimates, and discharge is controlled by rainfall, DOC flux and rainfall/discharge will always be well correlated.

Influence of drought-induced acidification on the mobility of dissolved organic carbon in peat soils

A strong relationship between dissolved organic carbon (DOC) and sulphate (SO42−) dynamics under drought conditions has been revealed from analysis of a 10-year time series (1993–2002). Soil solution from a blanket peat at 10 cm depth and stream water were collected at biweekly and weekly intervals, respectively, by the Environmental Change Network at Moor House-Upper Teesdale National Nature Reserve in the North Pennine uplands of Britain. DOC concentrations in soil solution and stream water were closely coupled, displaying a strong seasonal cycle with lowest concentrations in early spring and highest in late summer/early autumn. Soil solution DOC correlated strongly with seasonal variations in soil temperature at the same depth 4-weeks prior to sampling. Deviation from this relationship was seen, however, in years with significant water table drawdown (>−25 cm), such that DOC concentrations were up to 60% lower than expected. Periods of drought also resulted in the release of SO42−, because of the oxidation of inorganic/organic sulphur stored in the peat, which was accompanied by a decrease in pH and increase in ionic strength. As both pH and ionic strength are known to control the solubility of DOC, inclusion of a function to account for DOC suppression because of drought-induced acidification accounted for more of the variability of DOC in soil solution (R2=0.81) than temperature alone (R2=0.58). This statistical model of peat soil solution DOC at 10 cm depth was extended to reproduce 74% of the variation in stream DOC over this period. Analysis of annual budgets showed that the soil was the main source of SO42− during droughts, while atmospheric deposition was the main source in other years. Mass balance calculations also showed that most of the DOC originated from the peat. The DOC flux was also lower in the drought years of 1994 and 1995, reflecting low DOC concentrations in soil and stream water. The analysis presented in this paper suggests that lower concentrations of DOC in both soil and stream waters during drought years can be explained in terms of drought-induced acidification. As future climate change scenarios suggest an increase in the magnitude and frequency of drought events, these results imply potential for a related increase in DOC suppression by episodic acidification.

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.

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.

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.